- Timestamp:
- Jun 13, 2007, 4:36:08 PM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/pt_hedland_2006/run_pt_hedland_urs.py
r4467 r4542 30 30 from anuga.shallow_water import Field_boundary 31 31 from Numeric import allclose 32 from anuga.shallow_water.data_manager import export_grid 32 33 33 34 from anuga.pmesh.mesh_interface import create_mesh_from_regions 34 from anuga. abstract_2d_finite_volumes.utilimport start_screen_catcher, copy_code_files,store_parameters35 from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files,store_parameters 35 36 from anuga_parallel.parallel_api import distribute, numprocs, myid, barrier 36 37 from anuga_parallel.parallel_abstraction import get_processor_name 37 38 from anuga.caching import myhash 39 from anuga.damage_modelling.inundation_damage import add_depth_and_momentum2csv, inundation_damage 38 40 # Application specific imports 39 41 import project_urs # Definition of file names and polygons … … 65 67 kwargs['bathy_file']=project_urs.combined_dir_name + '.txt' 66 68 # kwargs['bathy_file']=project_urs.combined_small_dir_name + '.pts' 67 kwargs['boundary_file']=project_urs.boundaries_ dir_name + '.sww'69 kwargs['boundary_file']=project_urs.boundaries_in_dir_name + '.sww' 68 70 # kwargs['Completed']='' 71 72 start_screen_catcher(kwargs['output_dir'], myid, numprocs) 69 73 70 74 print 'output_dir',kwargs['output_dir'] … … 76 80 77 81 barrier() 78 79 start_screen_catcher(kwargs['output_dir'], myid, numprocs)80 82 81 83 print "Processor Name:",get_processor_name() … … 224 226 if myid == 0: 225 227 store_parameters(**kwargs) 226 barrier() 227 228 229 swwfile = kwargs['output_dir']+kwargs['aa_scenario_name'] 230 export_grid(swwfile, extra_name_out = 'town', 231 quantities = ['elevation','depth','stage','speed'], # '(xmomentum**2 + ymomentum**2)**0.5' defaults to elevation 232 timestep = None, 233 reduction = max, 234 cellsize = 25, 235 NODATA_value = -9999, 236 easting_min = project_urs.eastingmin, 237 easting_max = project_urs.eastingmax, 238 northing_min = project_urs.northingmin, 239 northing_max = project_urs.northingmax, 240 verbose = True, 241 origin = None, 242 datum = 'WGS84', 243 format = 'asc') 244 245 buildings_filename = project_urs.buildings_filename 246 buildings_filename_out = project_urs.buildings_filename_out 247 248 inundation_damage(swwfile+'.sww', buildings_filename, buildings_filename_out) 249 print '\n Augmented building file written to %s \n' %buildings_filename_out 250 251 barrier() 228 252 #------------------------------------------------------------- 229 253 if __name__ == "__main__":
Note: See TracChangeset
for help on using the changeset viewer.