Changeset 6263
- Timestamp:
- Feb 3, 2009, 2:49:34 PM (16 years ago)
- Location:
- anuga_work/production/busselton/standardised_version
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/busselton/standardised_version/build_busselton.py
r6262 r6263 1 """ 2 Script for building the elevation data to run a tsunami inundation scenario 3 for busselton, WA, Australia. 1 """Build the elevation data to run a tsunami inundation scenario 2 for Busselton, WA, Australia. 4 3 5 4 Input: elevation data from project.py … … 17 16 18 17 # Related major packages 19 from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts 18 from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf 19 from anuga.shallow_water.data_manager import dem2pts 20 from anuga.shallow_water.data_manager import start_screen_catcher 21 from anuga.shallow_water.data_manager import copy_code_files 20 22 from anuga.geospatial_data.geospatial_data import Geospatial_data 21 from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files22 23 23 24 # Application specific imports 24 25 import project # Definition of file names and polygons 26 25 27 26 28 #------------------------------------------------------------------------------ … … 28 30 # output to file 29 31 #------------------------------------------------------------------------------ 30 31 32 copy_code_files(project.output_build_time_dir,__file__, 32 os.path.dirname(project.__file__)+os.sep+ project.__name__+'.py' ) 33 os.path.dirname(project.__file__)+os.sep+\ 34 project.__name__+'.py' ) 33 35 start_screen_catcher(project.output_build_time_dir) 34 36 35 37 36 #------------------------------------------------------------------------------ -38 #------------------------------------------------------------------------------ 37 39 # Preparation of topographic data 38 40 # … … 40 42 # Do for coarse and fine data 41 43 # Fine pts file to be clipped to area of interest 42 #------------------------------------------------------------------------------ -44 #------------------------------------------------------------------------------ 43 45 44 46 print 'project.poly_all', project.poly_all … … 83 85 G_clipped.export_points_file(project.combined_dir_name + '.pts') 84 86 print 'Do txt version too' 85 G_clipped.export_points_file(project.combined_dir_name + '.txt') # Use for comparision in ARC 87 # Use for comparision in ARC 88 G_clipped.export_points_file(project.combined_dir_name + '.txt') 86 89 -
anuga_work/production/busselton/standardised_version/run_busselton.py
r6262 r6263 1 """ Script for running a tsunami inundation scenario for busselton, WA, Australia.1 """Run a tsunami inundation scenario for Busselton, WA, Australia. 2 2 3 The scenario is defined by a triangular mesh created from project.polygon, 4 the elevation data is compiled into a pts file through build_busselton.py 5 and asimulated tsunami is generated through an sts file from build_boundary.py.3 The scenario is defined by a triangular mesh created from project.polygon, the 4 elevation data is compiled into a pts file through build_busselton.py and a 5 simulated tsunami is generated through an sts file from build_boundary.py. 6 6 7 7 Input: sts file (build_boundary.py for respective event) … … 42 42 #----------------------------------------------------------------------- 43 43 # Copy scripts to time stamped output directory and capture screen 44 # output to file 44 # output to file. Copy script must be before screen_catcher 45 45 #----------------------------------------------------------------------- 46 47 copy script must be before screen_catcher48 46 copy_code_files(project.output_run_time_dir, __file__, 49 os.path.dirname(project.__file__)+os.sep+ project.__name__+'.py' ) 50 start_screen_catcher(project.output_run_time_dir, myid, numprocs) 47 os.path.dirname(project.__file__)+os.sep+\ 48 project.__name__+'.py' ) 49 start_screen_catcher(project.output_run_time_dir, 0, 1) 51 50 52 51 … … 89 88 domain.set_datadir(project.output_run_time_dir) 90 89 domain.set_minimum_storable_height(0.01) # Don't store depth less than 1cm 90 91 91 92 92 #------------------------------------------------------------------------- … … 140 140 141 141 print 'Simulation took %.2f seconds' %(time.time()-t0) 142 143 144 145 146 147
Note: See TracChangeset
for help on using the changeset viewer.