Changeset 6672
- Timestamp:
- Mar 31, 2009, 5:39:21 PM (16 years ago)
- Location:
- anuga_work/production/australia_ph2
- Files:
-
- 12 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/australia_ph2/adelaide/project.py
r6617 r6672 29 29 tide = 0 # difference between MSL and HAT 30 30 event_number = 64322 # the event number 31 # events: 7867 (altiplano), 27346 (java), 46591 (midAmerica), 32 # 64322 (sandwich), 65368 (shetland), 33 # other events: 51466 (newHebrides), 57537 (peru), 58334 (puysegur), 68791 (sumatra) 34 # 11100 (andaman), 31 # events: 27346 (java 2.1e-6), 58334 (puysegur 8.8e-7), 64322 (sandwich 8.2e-6), 32 # other events: 51466 (newHebrides 7.2e-7), 57537 (peru 4.7e-7), 46591 (midAmerica 1.3e-6), 33 # 68791 (sumatra 7.5e-7), 11100 (andaman 4.9e-7), 65368 (shetland 2.5e-6), 7867 (altiplano 1.0e-6) 35 34 alpha = 0.1 # smoothing parameter for mesh 36 35 friction=0.01 # manning's friction coefficient … … 42 41 wave = 'Tb' # Bf (sts wave) Tb (index wave) 43 42 44 setup = ' trial' # This can be one of three values43 setup = 'final' # This can be one of three values 45 44 # trial - coarsest mesh, fast 46 45 # basic - coarse mesh -
anuga_work/production/australia_ph2/albany/project.py
r6617 r6672 27 27 # One or all can be changed each time the run_model script is executed 28 28 tide = 0 # difference between MSL and HAT 29 central_meridian = 117.0 # central meridian for projection (optional) 29 30 event_number = 64344 # the event number or the mux file name 30 31 # selected events: 27319 (java), 64344 (sandwich), 46697 (midamerica), … … 47 48 # final - fine mesh, slowest 48 49 49 internal_polygon = True50 internal_polygon = False 50 51 51 52 #------------------------------------------------------------------------------- -
anuga_work/production/australia_ph2/broome/project.py
r6617 r6672 36 36 # index is only used when wave = Tb 37 37 index = 1693 # index from the PTHA - Y2000 5.4259m 38 wave = ' Bf' # Bf (sts wave) Tb (index wave)38 wave = 'Tb' # Bf (sts wave) Tb (index wave) 39 39 40 40 -
anuga_work/production/australia_ph2/carnarvon/export_results_max.py
r6412 r6672 22 22 directory = project.output_folder 23 23 24 time_dir1 = '20090 224_102606_run_trial_0_58362_kvanputt'25 ##time_dir2 = ''24 time_dir1 = '20090324_150124_run_final_0_27285_109_Bf_mhingee' 25 time_dir2 = '20090324_161452_run_final_0_27285_109_Tb_mhingee' 26 26 ##time_dir3 = '' 27 27 ##time_dir4 = '' … … 29 29 ##time_dir6 = '' 30 30 ## 31 time_dirs = [time_dir1 ] #time_dir2]# , time_dir4, time_dir5, time_dir6]31 time_dirs = [time_dir1, time_dir2] # , time_dir4, time_dir5, time_dir6] 32 32 33 33 … … 35 35 ##cellsize = 5 36 36 37 #timestep = None # None means no timestep!38 timestep = 037 timestep = None # None means no timestep! 38 ##timestep = 0 39 39 40 40 ###### … … 56 56 57 57 # one or more key strings from var_equations above 58 var = [' elevation']58 var = ['stage','speed'] 59 59 60 60 ###### … … 86 86 87 87 name1 = join(directory, time_dir, project.scenario_name) 88 #name2 = join(directory, time_dir, project.scenario_name)+'_time_39600_0' 89 #name3 = join(directory, time_dir, project.scenario_name)+'_time_79200_0' 90 91 names = [name1] #, name2]#, name3] 88 name2 = join(directory, time_dir, project.scenario_name)+'_time_8760_0' 89 name3 = join(directory, time_dir, project.scenario_name)+'_time_17520_0' 90 name4 = join(directory, time_dir, project.scenario_name)+'_time_26280_0' 91 name5 = join(directory, time_dir, project.scenario_name)+'_time_35040_0' 92 name6 = join(directory, time_dir, project.scenario_name)+'_time_43800_0' 93 name7 = join(directory, time_dir, project.scenario_name)+'_time_52560_0' 94 name8 = join(directory, time_dir, project.scenario_name)+'_time_61320_0' 95 name9 = join(directory, time_dir, project.scenario_name)+'_time_70080_0' 96 name10 = join(directory, time_dir, project.scenario_name)+'_time_78840_0' 97 98 names = [name1, name2, name3, name4, name5, name6, name7, name8, name9, name10] 92 99 93 #asc_name = []100 asc_name = [] 94 101 95 102 for name in names: … … 112 119 format = 'asc') 113 120 114 ##asc_name.append(outname + '.asc')115 ##116 ##maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'117 ## 118 ##print 'max asc outname ', maxasc_outname119 ##print 'asc_name ', str(asc_name)120 ##121 ##MaxAsc(maxasc_outname, asc_name)121 asc_name.append(outname + '.asc') 122 123 maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc' 124 125 print 'max asc outname ', maxasc_outname 126 print 'asc_name ', str(asc_name) 127 128 MaxAsc(maxasc_outname, asc_name) -
anuga_work/production/australia_ph2/carnarvon/project.py
r6617 r6672 44 44 # final - fine mesh, slowest 45 45 46 internal_polygon = False46 internal_polygon = True 47 47 48 48 #------------------------------------------------------------------------------- -
anuga_work/production/australia_ph2/ceduna/project.py
r6617 r6672 38 38 friction=0.01 # manning's friction coefficient 39 39 starttime=0 # start time for simulation 40 finaltime= 1000 # final time for simulation40 finaltime=80000 # final time for simulation 41 41 42 42 # index is only used when wave = Tb -
anuga_work/production/australia_ph2/darwin/export_results_max.py
r6532 r6672 22 22 directory = project.output_folder 23 23 24 time_dir1 = '200903 16_120226_run_final_0_17898_1833_Bf_mhingee'25 ##time_dir2 = ' 20090316_120344_run_final_0_64977_1833_Bf_mhingee'26 ##time_dir3 = ' 20090316_120457_run_final_0_70906_1833_Bf_mhingee'27 ##time_dir4 = ' 20090316_120558_run_final_0_71261_1833_Bf_mhingee'24 time_dir1 = '20090327_103428_run_final_0_70906_1833_Bf_None_mhingee' 25 ##time_dir2 = '' 26 ##time_dir3 = '' 27 ##time_dir4 = '' 28 28 ##time_dir5 = '' 29 29 ##time_dir6 = '' -
anuga_work/production/australia_ph2/darwin/project.py
r6617 r6672 28 28 tide = 0 # difference between MSL and HAT 29 29 central_meridian = 129.0 # for projection (necessary for nonstandard) 30 event_number = 70906# the event number or the mux file name30 event_number = 17898 # the event number or the mux file name 31 31 # contributing events: 17898 (Flores), 64977 (Seram), 70906 (Sumba_normal), 32 32 # 71261 (Timor_OTF) -
anuga_work/production/australia_ph2/hobart/project.py
r6476 r6672 27 27 # One or all can be changed each time the run_model script is executed 28 28 tide = 0 # difference between MSL and HAT 29 zone = 55 # specify zone of model 29 central_meridian = 147.0 # Central meridian for projection (optional) 30 30 31 31 alpha = 0.1 # smoothing parameter for mesh … … 35 35 36 36 # index is only used when wave = Tb 37 wave = 'Tb' # Bf (sts wave) Tb (index wave) 37 wave = 'Bf' # Bf (sts wave) Tb (index wave) 38 38 39 # for melbourne (west) side: 39 event_number = 58349 # the event number or the mux file name40 ##event_number = 58349 # the event number or the mux file name 40 41 # contributing sources: 41 42 # 51328 (newhebrides), 58349 (puysegur), 64535 (sandwich) 42 43 # other possible: 31966 (kermandec), 65350 (shetland) 43 index = 1983 # index from the PTHA - west side, close to melbourne; 1983 has 1 in 2000 wave height 0.534m44 ##index = 1983 # index from the PTHA - west side, close to melbourne; 1983 has 1 in 2000 wave height 0.534m 44 45 45 46 # for east, Hobart, side - point opposite hobart 46 ##event_number = 51202 # the event number or the mux file name 47 # contributing events: 64434 (sandwich), 58117 (puysegur), 51202 (new hebrides) 48 # Others: 31899 (kermandec), 46657 (midamerica), 57586 (peru), 49 ##index = 2473 # index from PTHA - east, opposite Hobart; 1.112m for 1 in 2000 47 event_number = 64461 # the event number or the mux file name 48 # contributing events: 51450 (newHebrides 7.3e-6), 58078 (puysegur 5.2e-6), 64461 (sandwich 8.2e-6) 49 # Others: 31816 (kermadec 4.7e-7), 46651 (midamerica 3.5e-7), 50 # 51 index = 2479 # index from PTHA - Y2000 1.0977m 52 53 # index = 2460 # for Bichen 50 54 51 55 setup = 'final' # This can be one of three values … … 53 57 # basic - coarse mesh 54 58 # final - fine mesh, slowest 59 60 internal_polygon = False 55 61 56 62 #------------------------------------------------------------------------------- … … 65 71 #------------------------------------------------------------------------------- 66 72 67 output_comment = [setup, tide, event_number, index, wave] 73 if internal_polygon: 74 internal_poly_comment = 'internal' 75 else: 76 internal_poly_comment = None 77 78 output_comment = [setup, tide, event_number, index, wave, internal_poly_comment] 68 79 69 80 #------------------------------------------------------------------------------- … … 77 88 78 89 # Format for point is x,y,elevation (with header) 79 point_filenames = ['mel_grid250m.txt', 80 'bight_250m.txt'] # 250m grid 2005 90 point_filenames = ['grid_250m.txt'] # 250m grid 2005 81 91 82 92 ### Add csv header list to all files in point_filenames … … 94 104 # Used in run_model.py 95 105 # Format for points easting,northing (no header) 96 interior_regions_data = [['tasmania.csv', 1000000]] 106 interior_regions_data = [] 107 108 # add an internal polygon to force different mesh generation 109 # used to test for discretisation error when building elevation 110 # make sure file is in same folder as interior regions and bouding polygon; 111 # format is same (2 column .csv; easting, northing; no header) 112 if internal_polygon: 113 interior_regions_data.append(['internal_polygon.csv', 114 bounding_polygon_maxarea]) 97 115 98 116 # LAND - used to set the initial stage/water to be offcoast only -
anuga_work/production/australia_ph2/hobart/run_model.py
r6511 r6672 156 156 verbose=True) 157 157 domain.set_boundary({'back': Br, 158 'side': B d,158 'side': Bt, 159 159 'ocean': Bf}) 160 160 … … 163 163 164 164 domain.set_boundary({'back': Br, 165 'side': B d,165 'side': Bt, 166 166 'ocean': Tb}) 167 167 else: -
anuga_work/production/australia_ph2/strahan/export_results_max.py
r6619 r6672 22 22 directory = project.output_folder 23 23 24 time_dir1 = '200903 02_104336_run_final_0_27255_mhingee'25 ##time_dir2 = ''26 ##time_dir3 = ''27 ##time_dir4 = ''28 ##time_dir5 = ''29 ##time_dir6 = ''30 ## 31 time_dirs = [time_dir1 ] #time_dir2]#, time_dir4, time_dir5, time_dir6]24 time_dir1 = '20090327_143655_run_final_0_64214_2044_Tb_None_mhingee' 25 time_dir2 = '20090327_143836_run_final_0_58337_2044_Tb_None_mhingee' 26 time_dir3 = '20090330_095007_run_final_0_68779_2044_Tb_None_mhingee' 27 time_dir4 = '20090330_095719_run_final_0_64214_2044_Bf_None_mhingee' 28 time_dir5 = '20090330_095845_run_final_0_58337_2044_Bf_None_mhingee' 29 time_dir6 = '20090330_143221_run_final_0_68779_2044_Bf_None_mhingee' 30 31 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6] 32 32 33 33 … … 36 36 37 37 #timestep = None # None means no timestep! 38 timestep = 038 timestep = None 39 39 40 40 ###### … … 56 56 57 57 # one or more key strings from var_equations above 58 var = [' elevation','stage']58 var = ['stage','speed'] 59 59 60 60 ###### … … 91 91 names = [name1] #, name2]#, name3] 92 92 93 #asc_name = []93 asc_name = [] 94 94 95 95 for name in names: … … 112 112 format = 'asc') 113 113 114 ##asc_name.append(outname + '.asc')115 ##116 ##maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'117 ## 118 ##print 'max asc outname ', maxasc_outname119 ##print 'asc_name ', str(asc_name)120 ##121 ##MaxAsc(maxasc_outname, asc_name)114 asc_name.append(outname + '.asc') 115 116 maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc' 117 118 print 'max asc outname ', maxasc_outname 119 print 'asc_name ', str(asc_name) 120 121 MaxAsc(maxasc_outname, asc_name) -
anuga_work/production/australia_ph2/strahan/project.py
r6619 r6672 28 28 tide = 0 # difference between MSL and HAT 29 29 central_meridian = 147.0 # Central meridian for projection (optional) 30 event_number = # the event number or the mux file name31 # main events: 64214 (sandwich 1.4e-5), 58337 (puysegur 3e-6), 30 event_number = 68779 # the event number or the mux file name 31 # main events: 64214 (sandwich 1.4e-5), 58337 (puysegur 3e-6), 68779 (sumatra 7.53-7) 32 32 # others: 27295 (java 3e-7), 31817 (kermadec 5e-7), 46559 (midAmerica 4.2e-7) 33 33 # 51460 (newHebrides 7e-7), 57477 (peru 5.4e-7), 65322 (shetland 3.6e-7), 34 # 68779 (sumatra 7.53-7) 34 35 35 alpha = 0.1 # smoothing parameter for mesh 36 36 friction=0.01 # manning's friction coefficient … … 40 40 # index is only used when wave = Tb 41 41 index = 2044 # index from the PTHA - Y2000 is 0.60374m 42 wave = ' Tb' # Bf (sts wave) Tb (index wave)42 wave = 'Bf' # Bf (sts wave) Tb (index wave) 43 43 44 44 setup = 'final' # This can be one of three values … … 77 77 78 78 # Format for point is x,y,elevation (with header) 79 point_filenames = ['grid250m_pts.txt', 80 'coast_250.txt'] # 250m grid 2005 79 point_filenames = ['grid_250m.txt'] # 250m grid 2005 81 80 82 81 ### Add csv header list to all files in point_filenames … … 125 124 # Thinned ordering file from Hazard Map (geographic) 126 125 # Format is index,latitude,longitude (with header) 127 urs_order_filename = 'urs_order _austtg.csv'126 urs_order_filename = 'urs_order.csv' 128 127 129 128 # Landward bounding points
Note: See TracChangeset
for help on using the changeset viewer.