Changeset 7082
- Timestamp:
- May 26, 2009, 10:46:11 AM (16 years ago)
- Location:
- anuga_work/production/hobart_2009
- Files:
-
- 5 added
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/hobart_2009/export_results_max.py
r6984 r7082 21 21 22 22 directory = project.output_folder 23 time_dir1 = '20090420_153539_run_final_0_58274_kvanputt' 24 time_dir2 = '20090420_153654_run_final_0_58279_kvanputt' 25 time_dir3 = '20090420_153739_run_final_0_58281_kvanputt' 26 time_dir4 = '20090420_153826_run_final_0_58078_kvanputt' 27 time_dir5 = '20090330_093828_run_final_0_58292_kvanputt' 28 time_dir6 = '20090504_123724_run_final_0_58282_small_kvanputt' 29 time_dir7 = '20090504_124158_run_final_0_58068_small_kvanputt' 23 time_dir1 = '20090505_150430_run_final_0.8_58292_None_kvanputt' 24 time_dir2 = '20090505_150517_run_final_0_58292_None_kvanputt' 25 time_dir3 = '20090505_150711_run_final_0_58280_None_kvanputt' 26 time_dir4 = '20090505_150805_run_final_0.8_58280_None_kvanputt' 27 time_dir5 = '20090505_151322_run_final_0.8_64477_None_kvanputt' 28 time_dir6 = '20090505_151447_run_final_0_64477_None_kvanputt' 30 29 31 time_dirs = [time_dir6, time_dir7] 32 #time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5] 30 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6] 33 31 34 32 … … 43 41 ###### 44 42 45 area = [' SArms', 'Bruny'] # strings must match keys in var_equations below43 area = ['Hobart', 'NW', 'South'] # strings must match keys in var_equations below 46 44 #area = ['All'] # 'All' means no special areas - the whole thing 47 45 … … 57 55 58 56 # one or more key strings from var_equations above 59 var = [' depth', 'speed']57 var = ['speed'] 60 58 61 59 ###### … … 86 84 for time_dir in time_dirs: 87 85 name1 = join(directory, time_dir, project.scenario_name) 88 name2 = join(directory, time_dir, project.scenario_name)+'_time_ 22680_0'89 name3 = join(directory, time_dir, project.scenario_name)+'_time_ 45360_0'90 name4 = join(directory, time_dir, project.scenario_name)+'_time_ 68040_0'86 name2 = join(directory, time_dir, project.scenario_name)+'_time_16980_0' 87 name3 = join(directory, time_dir, project.scenario_name)+'_time_33960_0' 88 name4 = join(directory, time_dir, project.scenario_name)+'_time_50940_0' 91 89 92 names = [name1 ] #, name2, name3, name4]90 names = [name1, name2, name3, name4] 93 91 94 ##asc_name = []92 asc_name = [] 95 93 96 94 for name in names: … … 113 111 format = 'asc') 114 112 115 ##asc_name.append(outname + '.asc')116 ##117 ##maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'118 ## 119 ##print 'max asc outname ', maxasc_outname120 ##print 'asc_name ', str(asc_name)121 ##122 ##MaxAsc(maxasc_outname, asc_name)113 asc_name.append(outname + '.asc') 114 115 maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc' 116 117 print 'max asc outname ', maxasc_outname 118 print 'asc_name ', str(asc_name) 119 120 MaxAsc(maxasc_outname, asc_name) -
anuga_work/production/hobart_2009/project.py
r6984 r7082 128 128 # Used in get_timeseries.py. 129 129 # Format easting,northing,name,elevation (with header) 130 gauges_filename = 't sunamipointsMGA.csv'130 gauges_filename = 'time_of_arrival.csv' #'TideGaugesPoints.csv' #'tsunamipointsMGA.csv' 131 131 132 132 # BUILDINGS EXPOSURE - for identifying inundated houses … … 134 134 # Format latitude,longitude etc (geographic) 135 135 building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS 136 137 # AREA OF IMAGES - Extent of each image to find out highest runup 138 # Used in get_runup.py 139 images_filename = 'images.csv' 136 140 137 141 # BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively … … 173 177 yminSArms = 5237028 174 178 ymaxSArms = 5238465 179 180 # ASCII export grid for Hobart Large 181 xminHobart = 520000 182 xmaxHobart = 545000 183 yminHobart = 5230000 184 ymaxHobart = 5260000 185 186 # ASCII export grid for North West Communities 187 xminNW = 548000 188 xmaxNW = 561000 189 yminNW = 5250000 190 ymaxNW = 5258000 191 192 # ASCII export grid for South Communities 193 xminSouth = 523000 194 xmaxSouth = 533300 195 yminSouth = 5197000 196 ymaxSouth = 5214000 175 197 176 198 … … 273 295 building_exposure = join(gauges_folder, building_exposure_filename) 274 296 297 # The absolute pathname for the image file 298 # Used for get_runup.py 299 if images_filename: 300 images = join(polygons_folder, images_filename) 301 302 275 303 # full path to where MUX files (or meta-files) live 276 304 mux_input = join(event_folder, mux_input_filename) -
anuga_work/production/hobart_2009/setup_model.py
r6984 r7082 154 154 project.interior_regions = [] 155 155 156 if not project.PriorityArea_filename ==None:156 if project.PriorityArea_filename is not None: 157 157 polygons, maxareas = csv2polygons(project.PriorityAreas) 158 158 print 'Creating %d internal polygons' % len(polygons)
Note: See TracChangeset
for help on using the changeset viewer.