source: anuga_work/production/pt_hedland_2006/build_pt_hedland.py @ 4631

Last change on this file since 4631 was 4587, checked in by nick, 17 years ago

update model scripts

File size: 4.9 KB
Line 
1"""Script for running tsunami inundation scenario for Dampier, WA, Australia.
2
3Source data such as elevation and boundary data is assumed to be available in
4directories specified by project_urs.py
5The output sww file is stored in project_urs.output_time_dir
6
7The scenario is defined by a triangular mesh created from project_urs.polygon,
8the elevation data and a simulated submarine landslide.
9
10Ole Nielsen and Duncan Gray, GA - 2005 and Jane Sexton, Nick Bartzis, GA - 2006
11"""
12
13#------------------------------------------------------------------------------
14# Import necessary modules
15#------------------------------------------------------------------------------
16
17# Standard modules
18from os import sep
19from os.path import dirname, basename
20from os import mkdir, access, F_OK
21from shutil import copy
22import time
23import sys
24
25# Related major packages
26from anuga.shallow_water import Domain
27from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts,urs_ungridded2sww
28from anuga.geospatial_data.geospatial_data import *
29from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files
30from anuga_parallel.parallel_abstraction import get_processor_name
31
32# Application specific imports
33import project_urs   # Definition of file names and polygons
34
35#------------------------------------------------------------------------------
36# Copy scripts to time stamped output directory and capture screen
37# output to file
38#------------------------------------------------------------------------------
39
40copy_code_files(project_urs.output_build_time_dir,__file__, 
41               dirname(project_urs.__file__)+sep+ project_urs.__name__+'.py' )
42
43start_screen_catcher(project_urs.output_build_time_dir)
44
45print "Processor Name:",get_processor_name()
46print 'USER: ', project_urs.user
47
48#-------------------------------------------------------------------------------
49# Preparation of topographic data
50#
51# Convert ASC 2 DEM 2 PTS using source data and store result in source data
52# Do for coarse and fine data
53# Fine pts file to be clipped to area of interest
54#-------------------------------------------------------------------------------
55print"project_urs.combined_dir_name",project_urs.combined_dir_name
56
57# topography directory filenames
58onshore_in_dir_name = project_urs.onshore_in_dir_name
59coast_in_dir_name = project_urs.coast_in_dir_name
60offshore_in_dir_name = project_urs.offshore_in_dir_name
61offshore1_in_dir_name = project_urs.offshore1_in_dir_name
62offshore2_in_dir_name = project_urs.offshore2_in_dir_name
63
64onshore_dir_name = project_urs.onshore_dir_name
65coast_dir_name = project_urs.coast_dir_name
66offshore_dir_name = project_urs.offshore_dir_name
67offshore1_dir_name = project_urs.offshore1_dir_name
68offshore2_dir_name = project_urs.offshore2_dir_name
69'''
70## creates DEM from asc data
71#print "creates DEMs from asc data"
72#convert_dem_from_ascii2netcdf(onshore_in_dir_name,
73#                              basename_out=onshore_dir_name,
74#                              use_cache=True, verbose=True)
75##creates pts file for onshore DEM
76#print "creates pts file for onshore DEM"
77#dem2pts(onshore_dir_name, use_cache=True, verbose=True)
78
79print'create Geospatial onshore objects from topographies'
80G = Geospatial_data(file_name = onshore_dir_name + '.pts') +\
81    Geospatial_data(file_name = coast_in_dir_name + '.txt') +\
82    Geospatial_data(file_name = offshore_in_dir_name + '.txt')
83     
84
85print'clip combined geospatial object by bounding polygon'
86G_clipped = G.clip(project_urs.poly_all)
87
88print'export combined DEM file'
89if access(project_urs.topographies_dir,F_OK) == 0:
90    mkdir (project_urs.topographies_dir)
91G.export_points_file(project_urs.combined_dir_name + '.txt')
92
93print'split combined data set'
94G_small, G_other = G_clipped.split(0.1, True)
95
96print'export split DEM file'
97G_small.export_points_file(project_urs.combined_dir_name + '_small' + '.txt')
98#G_other.export_points_file(project_urs.combined_dir_name + '_other' + '.pts')
99
100'''
101print 'start reading:',project_urs.combined_dir_name + '.txt'
102G = Geospatial_data(file_name = (project_urs.combined_dir_name + '.txt'))
103print 'start writing',project_urs.combined_dir_name + '.pts'
104G.export_points_file(project_urs.combined_dir_name + '.pts')
105'''
106G_clipped = G.clip(project_urs.poly_topo_clip)
107
108print 'start split'
109#G_smallest, G_other = G.split(0.1,True)
110
111print 'start export',project_urs.combined_smallest_dir_name + '.txt'
112#G.export_points_file(project_urs.combined_smaller_dir_name + '.txt')
113G_clipped.export_points_file(project_urs.combined_smallest_dir_name + '.txt')
114
115#-------------------------------------------------------------------------
116# Convert URS to SWW file for boundary conditions
117#-------------------------------------------------------------------------
118print 'starting to create boundary conditions'
119     
120print 'boundaries_in_dir_name',project_urs.boundaries_in_dir_name
121
122urs_ungridded2sww(project_urs.boundaries_in_dir_name, project_urs.boundaries_in_dir_name,
123                  verbose=True, mint=5000, maxt=35000, zscale=1)
124
125'''
126
127
128
Note: See TracBrowser for help on using the repository browser.