Changeset 7252 for anuga_work/production/new_south_wales
- Timestamp:
- Jun 23, 2009, 2:53:07 PM (16 years ago)
- Location:
- anuga_work/production/new_south_wales
- Files:
-
- 8 edited
-
batemans_bay/export_results_max.py (modified) (3 diffs)
-
batemans_bay/get_timeseries.py (modified) (1 diff)
-
batemans_bay/project.py (modified) (3 diffs)
-
gosford/Arc_extract_by_mask.py (modified) (4 diffs)
-
gosford/export_results_max.py (modified) (3 diffs)
-
gosford/get_timeseries.py (modified) (1 diff)
-
gosford/project.py (modified) (1 diff)
-
gosford/run_model.py (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/new_south_wales/batemans_bay/export_results_max.py
r7225 r7252 22 22 directory = project.output_folder 23 23 24 ##time_dir1 = '20081217_101310_run_final_0_27283_250m_none_dt_kvanputt' 25 ##time_dir2 = '20081217_115336_run_final_0_27283_250m_none_dp_kvanputt' 26 ###time_dir1 = '20080924_123626_run_final_0_27283_250m_all_kvanputt' # This uses the 250m bathymetry with all interior polygons 27 ###time_dir2 = '20080912_154716_run_final_0_27283_alpha0.1_kvanputt' # This uses original bathyemetry data 28 ###time_dir3 = '20080924_123601_run_final_0_27283_250m_none_kvanputt' # This uses the 250m bathymetry without any interior polygons 29 ##time_dirs = [time_dir1, time_dir2]#, time_dir3] 24 #Specify output directories 25 time_dir1 = '20090601_172248_run_final_0.0_51204_jgriffin' 26 time_dir2 = '20090529_143458_run_final_0.0_58346_jgriffin' 30 27 31 time_dir1 = '20090529_143527_run_final_0.0_51424_jgriffin' 32 ##time_dir2 = '20090529_143458_run_final_0.0_51292_jgriffin' 33 ##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt' 34 ##time_dir4 = '20081211_162433_run_final_0.6_68693_alpha0.1_kvanputt' 35 ##time_dir5 = '20081211_162656_run_final_0.6_27283_alpha0.1_kvanputt' 36 ##time_dir6 = '20081211_162744_run_final_0_27283_alpha0.1_kvanputt' 37 ## 38 time_dirs = [time_dir1]#,time_dir2] #5, time_dir4]# , time_dir4, time_dir5, time_dir6] 39 28 time_dirs = [time_dir1, time_dir2] 40 29 41 cellsize = 20 42 ##cellsize = 5 30 cellsize = 20 #dependent on data resolution in area of interest. 43 31 44 32 timestep = None # None means no timestep! … … 49 37 ###### 50 38 51 #area = [' Bunbury', 'Busselton'] # strings must match keys in var_equations below39 #area = ['', ''] # strings must match keys in var_equations below 52 40 area = ['All'] # 'All' means no special areas - the whole thing 53 41 … … 94 82 name1 = join(directory, time_dir, project.scenario_name) 95 83 name2 = join(directory, time_dir, project.scenario_name)+'_time_37860_0' 96 ## name3 = join(directory, time_dir, project.scenario_name)+'_time_12620_0'97 ## name4 = join(directory, time_dir, project.scenario_name)+'_time_18930_0'98 ## name5 = join(directory, time_dir, project.scenario_name)+'_time_25240_0'99 ## name6 = join(directory, time_dir, project.scenario_name)+'_time_31550_0'100 ## name7 = join(directory, time_dir, project.scenario_name)+'_time_37860_0'101 ## name8 = join(directory, time_dir, project.scenario_name)+'_time_44170_0'102 ## name9 = join(directory, time_dir, project.scenario_name)+'_time_50480_0'103 ## name10 = join(directory, time_dir, project.scenario_name)+'_time_56790_0'104 84 105 85 106 86 107 names = [name1, name2] #, name3,name4,name5,name6,name7,name8,name9,name10]87 names = [name1, name2] 108 88 109 89 asc_name = [] -
anuga_work/production/new_south_wales/batemans_bay/get_timeseries.py
r7225 r7252 20 20 time_dir1 = '20090529_143419_run_final_0.0_51077_jgriffin' 21 21 time_dir2 = '20090529_143442_run_final_0.0_51347_jgriffin' 22 time_dir3 = '20090529_143458_run_final_0.0_51292_jgriffin'23 time_dir4 = '20090529_143527_run_final_0.0_51424_jgriffin'24 time_dir5 = '20090529_143332_run_final_0.0_51378_jgriffin'25 22 26 #time_dir1 = '20090528_143200_run_final_1.0_51204_jgriffin' 27 #time_dir2 = '20090528_143323_run_final_1.0_58346_jgriffin' 28 #time_dir3 = '20090528_153025_run_final_0.0_58346_jgriffin' 29 #time_dir4 = '20090528_153109_run_final_0.0_51204_jgriffin' 30 ##time_dir5 = '20090520_145510_run_final_0.0_58129_jgriffin' 31 ##time_dir6 = '20090520_145608_run_final_0.0_58226_jgriffin' 32 ##time_dir7 = '20090520_145616_run_final_0.0_58284_jgriffin' 33 ##time_dir8 = '20090520_145700_run_final_0.0_58286_jgriffin' 34 ##time_dir9 = '20090525_093040_run_final_1.0_51445_jgriffin' 35 36 37 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5]#, time_dir6, time_dir7, time_dir8, time_dir9] 23 time_dirs = [time_dir1, time_dir2] 38 24 39 25 -
anuga_work/production/new_south_wales/batemans_bay/project.py
r7225 r7252 29 29 30 30 # the event number or the mux file name 31 ##event_number = 58129 #1 in 200 yr Puyesgur32 ##event_number = 58115 #1 in 500 yr Puysegur33 ##event_number = 58226 #1 in 1000 yr Puysegur34 ##event_number = 58284 #1 in 2000 yr Puysegur35 ##event_number = 58286 #1 in 5000 yr Puysegur36 ##event_number = 58346 #1 in 10000 yr Puysegur 37 38 ##event_number = 51077 #1 in 200 yr New Hebrides39 ##event_number = 51378 #1 in 500 yr New Hebrides40 ##event_number = 51347 #1 in 1000 yr New Hebrides41 ##event_number = 51292 #1 in 2000 yr New Hebrides42 event_number = 51424#1 in 5000 yr New Hebrides43 ##event_number = 51204 #1 in 10000 yr New Hebrides 44 45 ######event_number = 58368 #1 in 100 000 yr Puysegur46 ######event_number = 51436 #1 in 100 000 yr New Hebrides31 ##event_number = 58129 #1 in 200 yr Puyesgur 32 ##event_number = 58115 #1 in 500 yr Puysegur 33 ##event_number = 58226 #1 in 1000 yr Puysegur 34 ##event_number = 58284 #1 in 2000 yr Puysegur (Event 3) 35 ##event_number = 58286 #1 in 5000 yr Puysegur 36 event_number = 58346 #1 in 10000 yr Puysegur (Event 1) 37 38 ##event_number = 51077 #1 in 200 yr New Hebrides 39 ##event_number = 51378 #1 in 500 yr New Hebrides 40 ##event_number = 51347 #1 in 1000 yr New Hebrides 41 ##event_number = 51292 #1 in 2000 yr New Hebrides 42 ##event_number = 51424 #1 in 5000 yr New Hebrides 43 ##event_number = 51204 #1 in 10000 yr New Hebrides (Event 2) 44 45 ######event_number = 58368 #1 in 100 000 yr Puysegur 46 ######event_number = 51436 #1 in 100 000 yr New Hebrides 47 47 ######event_number = 58272 #1 in ~15000 yr Puysegur 48 ######event_number = 51445 #1 in ~15000 yr New Hebrides48 ######event_number = 51445 #1 in ~15000 yr New Hebrides 49 49 50 50 alpha = 0.1 # smoothing parameter for mesh … … 78 78 # Used in build_elevation.py 79 79 # Format for ascii grids, as produced in ArcGIS + a projection file 80 ascii_grid_filenames = ['1a',# Topographic data 81 '1b', 82 '2b', 83 '2a_3', 84 '3b', 85 '3a', 86 '4a_2', 87 '4b', 88 'off1', 89 'off2', 90 'off3', 91 'bbhd', 92 'sd100031996_p', 93 'sd100031996_p2', 94 'sd100031996_p3', 95 'sd100031996_p4'] 80 ascii_grid_filenames = [] 96 81 97 82 # Format for point is x,y,elevation (with header) 98 point_filenames = ['SD100031996_clip.csv', 99 'tomaga_offshore_AHD_MGA_1997.csv', 100 'Batemans_BBHD_MGA_1995.csv', 101 'moruya_AHD_MGA_2000.csv'] 83 point_filenames = [] 102 84 103 85 … … 135 117 # Used in run_building_inundation.py 136 118 # Format latitude,longitude etc (geographic) 137 ##building_exposure_filename = ' busselton_res_clip.csv' # from NEXIS119 ##building_exposure_filename = '' # from NEXIS 138 120 139 121 # AREA OF IMAGES - Extent of each image to find out highest runup 140 122 # Header - easting,northing,id,value 141 123 # Used in get_runup.py 142 images_filename = ' images.csv'124 images_filename = '' 143 125 144 126 # BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively -
anuga_work/production/new_south_wales/gosford/Arc_extract_by_mask.py
r7225 r7252 16 16 # Load required toolboxes... 17 17 gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Spatial Analyst Tools.tbx") 18 gp.overwriteoutput = 1 18 19 19 20 scenario_dir="\\\\nas2\\gemd\\georisk_models\\inundation\\data\\new_south_wales\\gosford_tsunami_scenario_2009\\" … … 21 22 22 23 23 variable = 'depth' 24 time_dir1 = '20090525_233002_run_final_1.0_58242_jgriffin' 24 variable = 'speed' 25 time_dir1 = '20090615_160527_run_final_1.0_51436_jgriffin' 26 time_dir2 = '20090615_160610_run_final_0.0_51436_jgriffin' 27 time_dir3 = '20090615_160900_run_final_0.0_58242_jgriffin' 28 time_dir4 = '20090615_160939_run_final_1.0_58242_jgriffin' 29 time_dir5 = '20090615_161039_run_final_0.0_58284_jgriffin' 30 time_dir6 = '20090615_161145_run_final_0.0_58349_jgriffin' 31 time_dir7 = '20090615_161318_run_final_0.0_58187_jgriffin' 32 time_dir8 = '20090615_161339_run_final_0.0_58113_jgriffin' 33 time_dir9 = '20090615_161420_run_final_0.0_58025_jgriffin' 25 34 26 time_dirs = [time_dir1] 27 35 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6, time_dir7, time_dir8, time_dir9] 28 36 29 37 for time_dir in time_dirs: 30 38 39 print time_dir 40 41 #get event number 42 time_dir_split = time_dir.split('_') 43 event_number = time_dir_split[5] 44 print 'event_number', event_number 45 31 46 # Local variables... 32 47 folder = scenario_dir + output_dir + time_dir + '\\' 33 raster_g bd= folder + 'raster.gdb'48 raster_gdb = folder + 'raster.gdb' 34 49 35 gp.Workspace = raster_g bd50 gp.Workspace = raster_gdb 36 51 37 52 print gp.Workspace 38 53 39 54 #replication dictionary 40 replicate = ((' gosford', ''),55 replicate = (('batemans_bay', ''), 41 56 ('_', ''),('Geordie', 'Geo'),('Sorrento', 'Sor'), ('max','M_'), 42 57 ('Fremantle', 'Fre'),('Rockingham', 'Roc'),('depth','_dep_'), … … 44 59 45 60 generate_filename = [] 46 input_raster = glob.glob(raster_gbd + 'All_dep*) 47 print time_dir 61 input_raster = raster_gdb + '\\All_spe_M_' 62 #input_raster = glob.glob(raster_gbd + 'All_dep*') 63 print 'input raster', input_raster 48 64 49 for infile in input_raster:50 output_raster = os.path.basename(infile) + 'clip'65 output_raster = raster_gdb + '\speed_clip' 66 print 'output_raster', output_raster 51 67 52 68 # Local variables... … … 55 71 #All_dep_M_ = "N:\\georisk_models\\inundation\\data\\new_south_wales\\batemans_bay_tsunami_scenario_2009\\anuga\\outputs\\20090528_153025_run_final_0.0_58346_jgriffin\\raster.gdb\\All_dep_M_" 56 72 #All_spe_M_ = "N:\\georisk_models\\inundation\\data\\new_south_wales\\batemans_bay_tsunami_scenario_2009\\anuga\\outputs\\20090528_153025_run_final_0.0_58346_jgriffin\\raster.gdb\\All_spe_M_" 57 aoi_umina_land_Union = "N:\\georisk_models\\inundation\\data\\new_south_wales\\gosford_tsunami_scenario_2009\\anuga\\polygons\\polygons.gdb\\aoi_umina_land_Union"58 #aoi_terrigal_umina_allarea_Union= "N:\\georisk_models\\inundation\\data\\new_south_wales\\gosford_tsunami_scenario_2009\\anuga\\polygons\\polygons.gdb\\aoi_terrigal_umina_allarea_Union"73 ## aoi_land = "N:\\georisk_models\\inundation\\data\\new_south_wales\\gosford_tsunami_scenario_2009\\anuga\\polygons\\polygons.gdb\\aoi_umina_land_union" 74 aoi_all = "N:\\georisk_models\\inundation\\data\\new_south_wales\\gosford_tsunami_scenario_2009\\anuga\\polygons\\polygons.gdb\\aoi_terrigal_umina_allarea_Union" 59 75 60 76 # Process: Extract by Mask... 61 gp.ExtractByMask_sa(input_raster, aoi_umina_land_Union, output_raster) 77 print 'extracting by mask' 78 gp.ExtractByMask_sa(input_raster, aoi_all, output_raster) 62 79 63 80 # Process: Extract by Mask (2)... 64 #gp.ExtractByMask_sa(All_spe_M_, aoi_ terrigal_umina_allarea_Union, speed_58346_msl)81 #gp.ExtractByMask_sa(All_spe_M_, aoi_all, speed_58346_msl) 65 82 -
anuga_work/production/new_south_wales/gosford/export_results_max.py
r7225 r7252 22 22 directory = project.output_folder 23 23 24 ##time_dir1 = '20081217_101310_run_final_0_27283_250m_none_dt_kvanputt' 25 ##time_dir2 = '20081217_115336_run_final_0_27283_250m_none_dp_kvanputt' 26 ###time_dir1 = '20080924_123626_run_final_0_27283_250m_all_kvanputt' # This uses the 250m bathymetry with all interior polygons 27 ###time_dir2 = '20080912_154716_run_final_0_27283_alpha0.1_kvanputt' # This uses original bathyemetry data 28 ###time_dir3 = '20080924_123601_run_final_0_27283_250m_none_kvanputt' # This uses the 250m bathymetry without any interior polygons 29 ##time_dirs = [time_dir1, time_dir2]#, time_dir3] 24 #Specify output directories 25 time_dir1 = '20090615_160610_run_final_0.0_51436_jgriffin' 26 time_dir2 = '20090615_161318_run_final_0.0_58187_jgriffin' 30 27 31 time_dir1 = '20090615_160610_run_final_0.0_51436_jgriffin' 32 ##time_dir2 = '20090615_161318_run_final_0.0_58187_jgriffin' 33 ##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt' 34 ##time_dir4 = '20081211_162433_run_final_0.6_68693_alpha0.1_kvanputt' 35 ##time_dir5 = '20081211_162656_run_final_0.6_27283_alpha0.1_kvanputt' 36 ##time_dir6 = '20081211_162744_run_final_0_27283_alpha0.1_kvanputt' 37 ## 38 time_dirs = [time_dir1]#, time_dir2] #5, time_dir4]# , time_dir4, time_dir5, time_dir6] 39 28 time_dirs = [time_dir1, time_dir2] 29 40 30 41 cellsize = 20 42 ##cellsize = 5 31 cellsize = 20 #dependent on data resolution in area of interest. 43 32 44 33 timestep = None # None means no timestep! … … 49 38 ###### 50 39 51 #area = [' Bunbury', 'Busselton'] # strings must match keys in var_equations below40 #area = ['', ''] # strings must match keys in var_equations below 52 41 area = ['All'] # 'All' means no special areas - the whole thing 53 42 … … 93 82 94 83 name1 = join(directory, time_dir, project.scenario_name) 95 ## name2 = join(directory, time_dir, project.scenario_name)+'_time_6310_0' 96 ## name3 = join(directory, time_dir, project.scenario_name)+'_time_12620_0' 97 ## name4 = join(directory, time_dir, project.scenario_name)+'_time_18930_0' 98 ## name5 = join(directory, time_dir, project.scenario_name)+'_time_25240_0' 99 ## name6 = join(directory, time_dir, project.scenario_name)+'_time_31550_0' 100 ## name7 = join(directory, time_dir, project.scenario_name)+'_time_37860_0' 101 ## name8 = join(directory, time_dir, project.scenario_name)+'_time_44170_0' 102 ## name9 = join(directory, time_dir, project.scenario_name)+'_time_50480_0' 103 ## name10 = join(directory, time_dir, project.scenario_name)+'_time_56790_0' 104 105 106 107 names = [name1]#, name2, name3,name4,name5,name6,name7,name8,name9,name10] 84 85 names = [name1] 108 86 109 87 asc_name = [] -
anuga_work/production/new_south_wales/gosford/get_timeseries.py
r7225 r7252 18 18 directory = project.output_folder 19 19 20 ##time_dir1 = '20090615_160527_run_final_1.0_51436_jgriffin' 21 ##time_dir2 = '20090615_160610_run_final_0.0_51436_jgriffin' 22 ##time_dir3 = '20090615_160900_run_final_0.0_58242_jgriffin' 23 ##time_dir4 = '20090615_160939_run_final_1.0_58242_jgriffin' 24 ##time_dir5 = '20090615_161039_run_final_0.0_58284_jgriffin' 25 ##time_dir6 = '20090615_161145_run_final_0.0_58349_jgriffin' 26 ##time_dir7 = '20090615_161318_run_final_0.0_58187_jgriffin' 27 ##time_dir8 = '20090615_161339_run_final_0.0_58113_jgriffin' 28 time_dir9 = '20090615_161420_run_final_0.0_58025_jgriffin' 20 time_dir1 = '20090615_160527_run_final_1.0_51436_jgriffin' 21 time_dir2 = '20090615_160610_run_final_0.0_51436_jgriffin' 29 22 30 31 time_dirs = [time_dir9] 32 ##time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6, time_dir7, time_dir8, time_dir9] 33 ##time_dirs = [time_dir7, time_dir8, time_dir9, time_dir10, time_dir11, time_dir12, time_dir13] 23 time_dirs = [time_dir1, time_dir2] 34 24 35 25 -
anuga_work/production/new_south_wales/gosford/project.py
r7251 r7252 107 107 # Header - easting,northing,id,value 108 108 # Used in get_runup.py 109 images_filename = ' images.csv'109 images_filename = '' 110 110 111 111 # BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively -
anuga_work/production/new_south_wales/gosford/run_model.py
r7225 r7252 158 158 159 159 t0 = time.time() 160 ##for t in domain.evolve(yieldstep=project.yieldstep,161 ## finaltime=1440,162 ## skip_initial_step=False):163 ## print domain.timestepping_statistics()164 ## print domain.boundary_statistics(tags='ocean')165 ##166 ##for t in domain.evolve(yieldstep=5,167 ## finaltime=6000,168 ## skip_initial_step=True):169 ## print domain.timestepping_statistics()170 ## print domain.boundary_statistics(tags='ocean')171 172 160 173 161 for t in domain.evolve(yieldstep=project.yieldstep,
Note: See TracChangeset
for help on using the changeset viewer.
