Changeset 7225 for anuga_work/production
- Timestamp:
- Jun 19, 2009, 12:38:12 PM (16 years ago)
- Location:
- anuga_work/production/new_south_wales
- Files:
-
- 5 added
- 10 edited
-
batemans_bay/Arc_asc2raster_GDA94z56.py (modified) (1 diff)
-
batemans_bay/Arc_txt2feature_class_runup.py (added)
-
batemans_bay/export_results_max.py (modified) (2 diffs)
-
batemans_bay/get_runup.py (added)
-
batemans_bay/get_timeseries.py (modified) (1 diff)
-
batemans_bay/project.py (modified) (4 diffs)
-
gosford/Arc_asc2raster_GDA94z56.py (modified) (1 diff)
-
gosford/Arc_extract_by_mask.py (added)
-
gosford/build_elevation.py (modified) (1 diff)
-
gosford/export_results_max.py (modified) (3 diffs)
-
gosford/get_runup.py (added)
-
gosford/get_timeseries.py (modified) (2 diffs)
-
gosford/project.py (modified) (4 diffs)
-
gosford/run_model.py (modified) (1 diff)
-
gosford/run_model_later.py (added)
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/new_south_wales/batemans_bay/Arc_asc2raster_GDA94z56.py
r7076 r7225 25 25 output_dir = "anuga\\outputs\\" 26 26 27 ##time_dir1 = '20081209_155431_run_final_0_27255_250m_none_lfountai' 28 ##time_dir2 = '20081210_100528_run_final_0_68693_250m_none_kvanputt' 29 ##time_dirs = [time_dir1, time_dir2] 30 31 time_dir1 = '20090520_144937_run_final_0.0_58272_jgriffin' 32 time_dir2 = '20090520_145453_run_final_0.0_58115_jgriffin' 33 time_dir3 = '20090520_145156_run_final_1.0_58272_jgriffin' 34 time_dir4 = '20090520_145208_run_final_0.0_51445_jgriffin' 35 time_dir5 = '20090520_145510_run_final_0.0_58129_jgriffin' 36 time_dir6 = '20090520_145608_run_final_0.0_58226_jgriffin' 37 time_dir7 = '20090520_145616_run_final_0.0_58284_jgriffin' 38 time_dir8 = '20090520_145700_run_final_0.0_58286_jgriffin' 39 time_dir9 = '20090525_093040_run_final_1.0_51445_jgriffin' 27 time_dir5 = '20090529_143332_run_final_0.0_51378_jgriffin' 28 time_dir1 = '20090529_143419_run_final_0.0_51077_jgriffin' 29 time_dir2 = '20090529_143442_run_final_0.0_51347_jgriffin' 30 time_dir3 = '20090529_143458_run_final_0.0_51292_jgriffin' 31 time_dir4 = '20090529_143527_run_final_0.0_51424_jgriffin' 40 32 41 33 42 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5 , time_dir6, time_dir7, time_dir8, time_dir9]34 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5]#, time_dir6, time_dir7, time_dir8, time_dir9] 43 35 44 36 for time_dir in time_dirs: -
anuga_work/production/new_south_wales/batemans_bay/export_results_max.py
r7139 r7225 29 29 ##time_dirs = [time_dir1, time_dir2]#, time_dir3] 30 30 31 time_dir1 = '2009052 8_153109_run_final_0.0_51204_jgriffin'32 ##time_dir2 = '200 81211_162311_run_final_0_27255_alpha0.1_kvanputt'31 time_dir1 = '20090529_143527_run_final_0.0_51424_jgriffin' 32 ##time_dir2 = '20090529_143458_run_final_0.0_51292_jgriffin' 33 33 ##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt' 34 34 ##time_dir4 = '20081211_162433_run_final_0.6_68693_alpha0.1_kvanputt' … … 36 36 ##time_dir6 = '20081211_162744_run_final_0_27283_alpha0.1_kvanputt' 37 37 ## 38 time_dirs = [time_dir1] #5, time_dir4]# , time_dir4, time_dir5, time_dir6]38 time_dirs = [time_dir1]#,time_dir2] #5, time_dir4]# , time_dir4, time_dir5, time_dir6] 39 39 40 40 -
anuga_work/production/new_south_wales/batemans_bay/get_timeseries.py
r7132 r7225 18 18 directory = project.output_folder 19 19 20 time_dir1 = '20090520_144937_run_final_0.0_58272_jgriffin' 21 time_dir2 = '20090520_145453_run_final_0.0_58115_jgriffin' 22 time_dir3 = '20090520_145156_run_final_1.0_58272_jgriffin' 23 time_dir4 = '20090520_145208_run_final_0.0_51445_jgriffin' 24 time_dir5 = '20090520_145510_run_final_0.0_58129_jgriffin' 25 time_dir6 = '20090520_145608_run_final_0.0_58226_jgriffin' 26 time_dir7 = '20090520_145616_run_final_0.0_58284_jgriffin' 27 time_dir8 = '20090520_145700_run_final_0.0_58286_jgriffin' 28 time_dir9 = '20090525_093040_run_final_1.0_51445_jgriffin' 20 time_dir1 = '20090529_143419_run_final_0.0_51077_jgriffin' 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 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' 29 35 30 36 31 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5 , time_dir6, time_dir7, time_dir8, time_dir9]37 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5]#, time_dir6, time_dir7, time_dir8, time_dir9] 32 38 33 39 -
anuga_work/production/new_south_wales/batemans_bay/project.py
r7139 r7225 43 43 ##event_number = 51204 #1 in 10000 yr New Hebrides 44 44 45 ######event_number = 58368 #1 in 100 000 yr Puysegur !!!46 ######event_number = 51436 #1 in 100 000 yr New Hebrides !!!47 ######event_number = 58272 #1 in 10000 yr Puysegur????48 ######event_number = 51445 #1 in 10000 yr New Hebrides????45 ######event_number = 58368 #1 in 100 000 yr Puysegur 46 ######event_number = 51436 #1 in 100 000 yr New Hebrides 47 ######event_number = 58272 #1 in ~15000 yr Puysegur 48 ######event_number = 51445 #1 in ~15000 yr New Hebrides 49 49 50 50 alpha = 0.1 # smoothing parameter for mesh … … 96 96 97 97 # Format for point is x,y,elevation (with header) 98 point_filenames = ['SD100031996_ jgriffin_clip.csv',98 point_filenames = ['SD100031996_clip.csv', 99 99 'tomaga_offshore_AHD_MGA_1997.csv', 100 100 'Batemans_BBHD_MGA_1995.csv', … … 136 136 # Format latitude,longitude etc (geographic) 137 137 ##building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS 138 139 # AREA OF IMAGES - Extent of each image to find out highest runup 140 # Header - easting,northing,id,value 141 # Used in get_runup.py 142 images_filename = 'images.csv' 138 143 139 144 # BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively … … 272 277 ##building_exposure = join(gauges_folder, building_exposure_filename) 273 278 279 # The absolute pathname for the image file 280 # Used for get_runup.py 281 if images_filename: 282 images = join(polygons_folder, images_filename) 283 274 284 # full path to where MUX files (or meta-files) live 275 285 mux_input = join(event_folder, mux_input_filename) -
anuga_work/production/new_south_wales/gosford/Arc_asc2raster_GDA94z56.py
r7069 r7225 29 29 ##time_dirs = [time_dir1, time_dir2] 30 30 31 time_dir1 = '20090521_114544_run_final_0.0_58242_jgriffin' 32 ##time_dir2 = '20090323_161308_run_final_0_58115_jgriffin' 33 ##time_dir3 = '20090323_161328_run_final_0_58226_jgriffin' 34 ##time_dir4 = '20090323_164424_run_final_0_58284_jgriffin' 35 ##time_dir5 = '20090323_164755_run_final_0_58286_jgriffin' 36 ##time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt' 31 time_dir1 = '20090615_160527_run_final_1.0_51436_jgriffin' 32 time_dir2 = '20090615_160610_run_final_0.0_51436_jgriffin' 33 time_dir3 = '20090615_160900_run_final_0.0_58242_jgriffin' 34 time_dir4 = '20090615_160939_run_final_1.0_58242_jgriffin' 35 time_dir5 = '20090615_161039_run_final_0.0_58284_jgriffin' 36 time_dir6 = '20090615_161145_run_final_0.0_58349_jgriffin' 37 time_dir7 = '20090615_161318_run_final_0.0_58187_jgriffin' 38 time_dir8 = '20090615_161339_run_final_0.0_58113_jgriffin' 39 time_dir9 = '20090615_161420_run_final_0.0_58025_jgriffin' 40 ##time_dir10 = '20090527_103636_run_final_0.0_58242_jgriffin' 41 ##time_dir11 = '20090527_103718_run_final_1.0_58242_jgriffin' 42 ##time_dir12 = '20090527_103748_run_final_1.0_51436_jgriffin' 43 ##time_dir13 = '20090527_103817_run_final_0.0_51436_jgriffin' 37 44 38 39 time_dirs = [time_dir1]#, time_dir2, time_dir3, time_dir4, time_dir5]#, time_dir6] 45 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6, time_dir7, time_dir8, time_dir9]#, time_dir10, time_dir11, time_dir12, time_dir13] 40 46 41 47 for time_dir in time_dirs: -
anuga_work/production/new_south_wales/gosford/build_elevation.py
r7069 r7225 55 55 convert_dem_from_ascii2netcdf(absolute_filename, 56 56 basename_out=absolute_filename, 57 use_cache= False,57 use_cache=True, 58 58 verbose=True) 59 59 dem2pts(absolute_filename, use_cache=True, verbose=True) -
anuga_work/production/new_south_wales/gosford/export_results_max.py
r7132 r7225 29 29 ##time_dirs = [time_dir1, time_dir2]#, time_dir3] 30 30 31 time_dir1 = '20090 526_091618_run_final_0.0_58025_jgriffin'32 ##time_dir2 = '20090 526_091528_run_final_0.0_58113_jgriffin'31 time_dir1 = '20090615_160610_run_final_0.0_51436_jgriffin' 32 ##time_dir2 = '20090615_161318_run_final_0.0_58187_jgriffin' 33 33 ##time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt' 34 34 ##time_dir4 = '20081211_162433_run_final_0.6_68693_alpha0.1_kvanputt' … … 36 36 ##time_dir6 = '20081211_162744_run_final_0_27283_alpha0.1_kvanputt' 37 37 ## 38 time_dirs = [time_dir1]# time_dir2] #5, time_dir4]# , time_dir4, time_dir5, time_dir6]38 time_dirs = [time_dir1]#, time_dir2] #5, time_dir4]# , time_dir4, time_dir5, time_dir6] 39 39 40 40 … … 63 63 64 64 # one or more key strings from var_equations above 65 var = ['depth', 'speed','stage']66 65 #var = ['depth', 'speed','stage'] 66 var = ['depth', 'speed', 'stage'] 67 67 ###### 68 68 # Start running the various conversions we require. -
anuga_work/production/new_south_wales/gosford/get_timeseries.py
r7069 r7225 18 18 directory = project.output_folder 19 19 20 ##time_dir1 = '20081009_184721_run_final_0.6_68693_alpha0.1_kvanputt'+sep+'gauge' 21 ##time_dir2 = '20081031_120316_run_final_0_68693_alpha0.1_kvanputt'+sep+'gauge' 22 ##time_dir3 = '20081031_120515_run_final_0_27255_alpha0.1_kvanputt'+sep+'gauge' 23 ##time_dir4 = '20081031_120642_run_final_0_27283_alpha0.1_kvanputt'+sep+'gauge' 24 ##time_dir5 = '20081031_120753_run_final_0.6_27283_alpha0.1_kvanputt'+sep+'gauge' 25 ##time_dir6 = '20081111_113832_run_final_0.6_27255_alpha0.1_kvanputt'+sep+'gauge' 26 ## 27 ##time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6] 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' 28 29 29 ##time_dir1 = '20081217_115336_run_final_0_27283_250m_none_dp_kvanputt'30 time_dir6 = '20090306_165253_run_final_0_58368_jgriffin'31 30 32 time_dirs = [time_dir6] 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] 34 33 35 34 36 for time_dir in time_dirs: … … 36 38 gauge = project.gauges 37 39 sww2csv_gauges(name+'.sww',gauge, 38 quantities = ['stage','speed','depth' ,'elevation'],40 quantities = ['stage','speed','depth'], 39 41 verbose=True) 40 42 -
anuga_work/production/new_south_wales/gosford/project.py
r7132 r7225 9 9 from os.path import join, exists 10 10 import csv 11 #import run_model_later 11 12 12 13 … … 25 26 #------------------------------------------------------------------------------- 26 27 27 28 # Model specific parameters. 29 # One or all can be changed each time the run_model script is executed 30 tide = 0.0 # difference between MSL and HAT in metres 28 tide = 0.0 29 # difference between MSL and HAT in metres 31 30 zone = 56 # specify zone of model 32 31 #event_number = 58242 # Puysegur 1 in 10 000 # the event number or the mux file name 33 event_number = 51436 # New Hebrides 1 in 10 00032 #event_number = 51436 # New Hebrides 1 in 10 000 34 33 #event_number = 58349 # Puysegur 1 in 5000 35 34 #event_number = 58284 # Puysegur 1 in 2000 36 35 #event_number = 58187 # Puysegur 1 in 1000 37 36 #event_number = 58113 # Puysegur 1 in 500 38 #event_number = 58025 # Puysegur 1 in 200 39 37 event_number = 58025 # Puysegur 1 in 200 40 38 41 39 alpha = 0.1 # smoothing parameter for mesh … … 144 142 ##building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS 145 143 144 # AREA OF IMAGES - Extent of each image to find out highest runup 145 # Header - easting,northing,id,value 146 # Used in get_runup.py 147 images_filename = 'images.csv' 148 146 149 # BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively 147 150 # NOTE: when files are put together the points must be in sequence … … 279 282 ##building_exposure = join(gauges_folder, building_exposure_filename) 280 283 284 # The absolute pathname for the image file 285 # Used for get_runup.py 286 if images_filename: 287 images = join(polygons_folder, images_filename) 288 281 289 # full path to where MUX files (or meta-files) live 282 290 mux_input = join(event_folder, mux_input_filename) -
anuga_work/production/new_south_wales/gosford/run_model.py
r7069 r7225 176 176 print domain.timestepping_statistics() 177 177 print domain.boundary_statistics(tags='ocean') 178 print domain.volumetric_balance_statistics() 178 179 179 180 print 'Simulation took %.2f seconds' % (time.time()-t0)
Note: See TracChangeset
for help on using the changeset viewer.
