Changeset 4546
- Timestamp:
- Jun 19, 2007, 10:56:41 AM (17 years ago)
- Location:
- anuga_work/production
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/broome_2006/project_urs.py
r4544 r4546 32 32 33 33 #tide = -5.3 34 tide = 035 #tide = 4.934 #tide = 0 35 tide = 4.9 36 36 37 37 #Maybe will try to make project a class to allow these parameters to be passed in. … … 41 41 finaltime=25000 42 42 setup='final' 43 source=' dampier'43 source='onslow' 44 44 45 45 if setup =='trial': -
anuga_work/production/broome_2006/run_broome_urs.py
r4544 r4546 38 38 from anuga.caching import myhash 39 39 from anuga.damage_modelling.inundation_damage import add_depth_and_momentum2csv, inundation_damage 40 from anuga.fit_interpolate.benchmark_least_squares import mem_usage 40 41 41 42 # Application specific imports … … 44 45 def run_model(**kwargs): 45 46 46 # tide = kwargs['tide'] 47 # alpha = kwargs['alpha'] 48 # friction = kwargs['friction'] 49 # time_thinning = kwargs['time_thinning'] 50 # scenario_name = kwargs['aa_scenario_name'] 51 52 53 47 54 48 #------------------------------------------------------------------------------ 55 49 # Copy scripts to time stamped output directory and capture screen … … 71 65 # kwargs['bathy_file']=project_urs.combined_small_dir_name + '.pts' 72 66 kwargs['boundary_file']=project_urs.boundaries_in_dir_name + '.sww' 73 start_screen_catcher(kwargs['output_dir'], myid, numprocs)74 67 75 68 print 'output_dir',kwargs['output_dir'] … … 82 75 barrier() 83 76 84 #start_screen_catcher(kwargs['output_dir'], myid, numprocs)77 start_screen_catcher(kwargs['output_dir'], myid, numprocs) 85 78 86 79 print "Processor Name:",get_processor_name() … … 125 118 #above don't work 126 119 domain = Domain(project_urs.meshes_dir_name+'.msh', use_cache=False, verbose=True) 120 print 'memory usage before del domain',mem_usage() 127 121 128 122 print domain.statistics() … … 234 228 if myid == 0: 235 229 store_parameters(**kwargs) 236 barrier() 237 238 swwfile = kwargs['output_dir']+kwargs['aa_scenario_name'] 239 export_grid(swwfile, extra_name_out = 'town', 230 231 print 'memory usage before del domain',mem_usage() 232 del domain 233 print 'memory usage after del domain',mem_usage() 234 235 swwfile = kwargs['output_dir']+kwargs['aa_scenario_name'] 236 export_grid(swwfile, extra_name_out = 'town', 240 237 quantities = ['elevation','stage','speed'], # '(xmomentum**2 + ymomentum**2)**0.5' defaults to elevation 241 238 timestep = None, … … 252 249 format = 'asc') 253 250 254 buildings_filename = project_urs.buildings_filename255 buildings_filename_out = project_urs.buildings_filename_out251 buildings_filename = project_urs.buildings_filename 252 buildings_filename_out = project_urs.buildings_filename_out 256 253 257 inundation_damage(swwfile+'.sww', buildings_filename, buildings_filename_out)258 print '\n Augmented building file written to %s \n' %buildings_filename_out259 254 inundation_damage(swwfile+'.sww', buildings_filename, buildings_filename_out) 255 print '\n Augmented building file written to %s \n' %buildings_filename_out 256 barrier() 260 257 261 258 #------------------------------------------------------------- -
anuga_work/production/dampier_2006/project_urs.py
r4544 r4546 25 25 26 26 #tide = -2.5 27 tide = 0.028 #tide = 2.427 #tide = 0.0 28 tide = 2.4 29 29 30 30 #Maybe will try to make project a class to allow these parameters to be passed in. -
anuga_work/production/exmouth_2006/project.py
r4544 r4546 30 30 31 31 #tide = -1.4 32 tide = 033 #tide = 1.432 #tide = 0 33 tide = 1.4 34 34 35 35 #Maybe will try to make project a class to allow these parameters to be passed in. … … 39 39 finaltime=25000 40 40 setup='final' 41 extra='dampier' 41 source='onslow' 42 #source='dampier' 43 #source='exmouth' 42 44 43 45 if setup =='trial': … … 63 65 64 66 65 dir_comment='_'+setup+'_'+str(tide)+'_'+str( extra)+'_'+str(user)67 dir_comment='_'+setup+'_'+str(tide)+'_'+str(source)+'_'+str(user) 66 68 67 69 onshore_name = 'DLI' … … 111 113 tide_dir = anuga_dir+'tide_data'+sep 112 114 115 if source =='dampier': 116 boundaries_name = 'exmouth_3854_17042007' #Dampier gun 117 boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'dampier'+sep+'1_10000'+sep 118 119 if source=='onslow': 120 boundaries_name = 'exmouth_3103_18052007' #onslow_hedland_broome gun 121 boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'onslow_hedland_broome'+sep+'1_10000'+sep 122 123 if source=='exmouth': 124 boundaries_name = 'exmouth_3859_16052007' #exmouth gun 125 boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep 126 113 127 #boundaries locations 114 boundaries_name = 'exmouth_3854_17042007' #dampier 115 #boundaries_name = 'exmouth_3103_18052007' #exmouth 116 #boundaries_name = 'exmouth_3859_16052007' #onslow 117 118 boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'dampier'+sep+'1_10000'+sep 119 #boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep 120 #boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'onslow_hedland_broome'+sep+'1_10000'+sep 128 boundaries_source = '' 129 #boundaries_in_dir = anuga_dir+'boundaries'+sep 121 130 boundaries_in_dir_name = boundaries_in_dir + boundaries_name 122 boundaries_dir = anuga_dir+'boundaries'+sep123 boundaries_dir_name = boundaries_dir + boundaries_name124 131 125 132 #output locations … … 161 168 res_coast = 10000*res_factor 162 169 163 poly_coast1 = read_polygon(polygons_dir+'pts_coast1.csv')170 #poly_coast1 = read_polygon(polygons_dir+'pts_coast1.csv') 164 171 #res_0 = 20000*res_factor 165 172 … … 197 204 res_airport = 5000*res_factor 198 205 199 interior_regions = [[poly_coast,res_coast],[poly_coast 1,res_coast],[poly_coast2,res_coast]206 interior_regions = [[poly_coast,res_coast],[poly_coast2,res_coast] 200 207 ,[poly_coast3,res_coast],[poly_middle,res_middle],[poly_exmouth,res_exmouth] 201 208 ,[poly_exmouth_south,res_exmouth],[poly_carpark1,res_exmouth],[poly_carpark2,res_exmouth] -
anuga_work/production/exmouth_2006/run_exmouth.py
r4542 r4546 38 38 from anuga.caching import myhash 39 39 from anuga.damage_modelling.inundation_damage import add_depth_and_momentum2csv, inundation_damage 40 from anuga.fit_interpolate.benchmark_least_squares import mem_usage 40 41 41 42 # Application specific imports … … 44 45 def run_model(**kwargs): 45 46 46 scenario_name = kwargs['aa_scenario_name']47 #scenario_name = kwargs['aa_scenario_name'] 47 48 48 49 #------------------------------------------------------------------------------ … … 107 108 use_cache=False, 108 109 verbose=True) 110 copy(project.meshes_dir_name+'.msh',project.output_dir+project.scenario_name+'.msh') 109 111 barrier() 110 112 … … 226 228 if myid == 0: 227 229 store_parameters(**kwargs) 228 229 swwfile = kwargs['output_dir']+kwargs['aa_scenario_name'] 230 export_grid(swwfile, extra_name_out = 'town', 230 231 print 'memory usage before del domain',mem_usage() 232 del domain 233 print 'memory usage after del domain',mem_usage() 234 235 swwfile = kwargs['output_dir']+kwargs['aa_scenario_name'] 236 export_grid(swwfile, extra_name_out = 'town', 231 237 # quantities = ['depth'], # '(xmomentum**2 + ymomentum**2)**0.5' defaults to elevation 232 238 quantities = ['elevation','depth','stage','speed'], # '(xmomentum**2 + ymomentum**2)**0.5' defaults to elevation … … 244 250 format = 'asc') 245 251 246 buildings_filename = project.buildings_filename247 buildings_filename_out = project.buildings_filename_out252 buildings_filename = project.buildings_filename 253 buildings_filename_out = project.buildings_filename_out 248 254 249 inundation_damage(swwfile+'.sww', buildings_filename, buildings_filename_out)250 print '\n Augmented building file written to %s \n' %buildings_filename_out255 inundation_damage(swwfile+'.sww', buildings_filename, buildings_filename_out) 256 print '\n Augmented building file written to %s \n' %buildings_filename_out 251 257 252 258 barrier()
Note: See TracChangeset
for help on using the changeset viewer.