Changeset 7364
- Timestamp:
- Aug 14, 2009, 10:30:16 AM (15 years ago)
- Location:
- anuga_work/production/gold_coast_2009
- Files:
-
- 1 added
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/gold_coast_2009/Arc_asc2raster_GDA94z56.py
r7089 r7364 29 29 time_dir3 = '20090511_161526_run_final_1.1_51469_kvanputt' 30 30 time_dir4 = '20090518_154710_run_final_0_50994_lfountai' 31 ##time_dir5 = '20090519_160510_run_final_1.1_50994_lfountai'32 ##time_dir6 = '20090521_220101_run_final_1.1_50863_kvanputt'33 ##time_dir7 = '20090522_164526_run_final_0_51392_lfountai'34 ##time_dir8 = '20090522_164640_run_final_1.1_51392_lfountai'35 ##time_dir9 = '20090522_164948_run_final_0_51423_lfountai'36 ##time_dir10 = '20090522_165600_run_final_1.1_51423_lfountai'31 time_dir5 = '20090519_160510_run_final_1.1_50994_lfountai' 32 time_dir6 = '20090521_220101_run_final_1.1_50863_kvanputt' 33 time_dir7 = '20090522_164526_run_final_0_51392_lfountai' 34 time_dir8 = '20090522_164640_run_final_1.1_51392_lfountai' 35 time_dir9 = '20090522_164948_run_final_0_51423_lfountai' 36 time_dir10 = '20090522_165600_run_final_1.1_51423_lfountai' 37 37 38 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4] #, time_dir5, time_dir6] 38 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, 39 time_dir6, time_dir7, time_dir8, time_dir9, time_dir10] 39 40 40 41 for time_dir in time_dirs: … … 44 45 raster_gbd = folder + 'raster.gdb' 45 46 #contour = raster_gbd + '\\contour_dep' 46 ## land = scenario_dir + "ArcGIS\\gold_coast.gdb\\Coast_250"47 land = scenario_dir + "ArcGIS\\gold_coast.gdb\\Land_initial_conditions" 47 48 ## ocean = scenario_dir + "map_work\\Perth.gdb\\Outlines\\initial_condition_ocean" 48 49 … … 56 57 #replication dictionary 57 58 replicate = (('gold_coast', ''), ('_time_29220_0', 'b'), ('_time_58440_0', 'c'), 59 ('_time_28860_0', 'b'), ('_time_57720_0', 'c'), 58 60 ('_', ''), ('max','M_'), ('depth','_dep_'), 59 61 ('speed', '_spe_'), ('elevation', '_ele_'), ('stage','_stage')) … … 84 86 85 87 86 ##output_extract = output_DEM + 'E'87 ##print 'Output Extract ',output_extract88 ##89 ##print 'Process: Extract by Mask'90 ##gp.ExtractByMask_sa(output_DEM, land, output_extract)91 ## 88 output_extract = output_DEM + 'E' 89 print 'Output Extract ',output_extract 90 91 print 'Process: Extract by Mask' 92 gp.ExtractByMask_sa(output_DEM, land, output_extract) 93 92 94 -
anuga_work/production/gold_coast_2009/For_DVD/project.py
r7306 r7364 91 91 interior_regions_data = [['area_of_interest.csv', 500], 92 92 ['intermediate.csv', 25000]] 93 PriorityArea_filename = 'PriorityAreas.csv' 93 # If there are several priority areas of interest these can be defined 94 # in a single .csv file with headers 'easting, northing, id, value' where 95 # id defines each polygon, and the value is the mesh resolution 96 PriorityArea_filename = None 94 97 95 98 # LAND - used to set the initial stage/water to be offcoast only … … 137 140 # Environment variable names. 138 141 # The inundation directory, not the data directory. 139 ENV_INUNDATIONHOME = ' INUNDATIONHOME'142 ENV_INUNDATIONHOME = 'ANUGADATA' 140 143 141 144 #------------------------------------------------------------------------------- … … 159 162 # create paths generated from environment variables. 160 163 home = join(os.getenv(ENV_INUNDATIONHOME), 'data') # Absolute path for data folder 161 muxhome = os.getenv(ENV_MUXHOME)162 164 163 165 # check various directories/files that must exist -
anuga_work/production/gold_coast_2009/For_DVD/run_model.py
r7306 r7364 1 """Run a tsunami inundation scenario for Busselton, WA, Australia.1 """Run a tsunami inundation scenario for Gold Coast, Queensland, Australia. 2 2 3 3 The scenario is defined by a triangular mesh created from project.polygon, the -
anuga_work/production/gold_coast_2009/For_DVD/setup_model.py
r7306 r7364 11 11 from os.path import join, exists 12 12 from anuga.utilities.polygon import read_polygon, number_mesh_triangles 13 13 from anuga.shallow_water.data_manager import csv2polygons 14 14 import project 15 15 -
anuga_work/production/gold_coast_2009/export_results_max.py
r7089 r7364 22 22 directory = project.output_folder 23 23 24 time_dir1 = '20090521_220101_run_final_1.1_50863_kvanputt' 25 time_dir2 = '20090522_164526_run_final_0_51392_lfountai' 26 time_dir3 = '20090522_164640_run_final_1.1_51392_lfountai' 27 time_dir4 = '20090522_164948_run_final_0_51423_lfountai' 28 time_dir5 = '20090522_165600_run_final_1.1_51423_lfountai 24 time_dir1 = '20090605_120328_run_final_0_250m_51469_lfountai' 25 ##time_dir1 = '20090522_164640_run_final_1.1_51392_lfountai' 29 26 30 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]31 27 32 ##cellsize = 250 33 cellsize = 20 28 time_dirs = [time_dir1] #, time_dir2] #, time_dir3, time_dir4, time_dir5, time_dir6] 29 30 cellsize = 250 31 ##cellsize = 20 34 32 35 33 #timestep = None # None means no timestep! … … 84 82 85 83 name1 = join(directory, time_dir, project.scenario_name) 86 name2 = join(directory, time_dir, project.scenario_name)+'_time_29220_0'87 name3 = join(directory, time_dir, project.scenario_name)+'_time_58440_0'84 ## name2 = join(directory, time_dir, project.scenario_name)+'_time_29220_0' 85 ## name3 = join(directory, time_dir, project.scenario_name)+'_time_58440_0' 88 86 89 names = [name1 , name2, name3]87 names = [name1] #, name2, name3] 90 88 91 89 asc_name = [] … … 112 110 asc_name.append(outname + '.asc') 113 111 114 maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'115 116 print 'max asc outname ', maxasc_outname117 print 'asc_name ', str(asc_name)118 119 MaxAsc(maxasc_outname, asc_name)112 ## maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc' 113 ## 114 ## print 'max asc outname ', maxasc_outname 115 ## print 'asc_name ', str(asc_name) 116 ## 117 ## MaxAsc(maxasc_outname, asc_name) -
anuga_work/production/gold_coast_2009/get_timeseries.py
r7004 r7364 19 19 directory = project.output_folder 20 20 21 time_dir1 = '20090327_135943_run_final_0_51469_lfountai' 22 ##time_dir2 = '20090327_153206_run_final_0_51348_lfountai' 23 ##time_dir3 = '20090327_153231_run_final_0_51253_lfountai' 24 ##time_dir4 = '20090327_153304_run_final_0_50863_lfountai' 25 time_dir5 = '20090327_153403_run_final_0_51423_lfountai' 26 time_dir6 = '20090327_153602_run_final_0_62833_lfountai' 27 time_dirs = [time_dir1, time_dir5, time_dir6] #time_dir2, time_dir3, time_dir4, time_dir5, time_dir6] 21 time_dir1 = '20090508_150215_run_final_0_51469_lfountai' 22 time_dir2 = '20090511_161526_run_final_1.1_51469_kvanputt' 23 time_dir3 = '20090511_165539_run_final_0_50863_lfountai' 24 time_dir4 = '20090521_220101_run_final_1.1_50863_kvanputt' 25 time_dir5 = '20090522_164526_run_final_0_51392_lfountai' 26 time_dir6 = '20090522_164640_run_final_1.1_51392_lfountai' 27 time_dir7 = '20090518_154710_run_final_0_50994_lfountai' 28 time_dir8 = '20090519_160510_run_final_1.1_50994_lfountai' 29 time_dir9 = '20090522_164948_run_final_0_51423_lfountai' 30 time_dir10 = '20090522_165600_run_final_1.1_51423_lfountai' 31 32 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, 33 time_dir6, time_dir7, time_dir8, time_dir9, time_dir10] 28 34 29 35 for time_dir in time_dirs: -
anuga_work/production/gold_coast_2009/project.py
r7247 r7364 34 34 event_number = int(sys.argv[1]) 35 35 else: 36 event_number = 51 469# the event number or the mux file name36 event_number = 51253 # the event number or the mux file name 37 37 38 38 event_number_list = [51469, 51392, 50863] # To piggy back multiple events
Note: See TracChangeset
for help on using the changeset viewer.