Changeset 7179
- Timestamp:
- Jun 10, 2009, 6:32:03 PM (16 years ago)
- Location:
- anuga_work/production/Broome_2009
- Files:
-
- 2 edited
- 1 copied
-
. (copied) (copied from anuga_work/production/Broome_2008)
-
project.py (modified) (9 diffs)
-
run_model.py (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/Broome_2009/project.py
r7161 r7179 10 10 11 11 #------------------------------------------------------------------------------- 12 # Fitting Scenario13 #-------------------------------------------------------------------------------14 15 extent = None #'small' #None #'sw' # 'aos1' ##'aos1'16 extent_maxarea = 3000017 18 #-------------------------------------------------------------------------------19 12 # Directory setup 20 13 #------------------------------------------------------------------------------- … … 22 15 # this section needs to be updated to reflect the modelled community. 23 16 # Note, the user needs to set up the directory system accordingly 24 state = ' tasmania'25 scenario_name = ' hobart'26 scenario_folder = ' hobart_tsunami_scenario_2009'17 state = 'western_australia' 18 scenario_name = 'broome' 19 scenario_folder = 'broome_tsunami_scenario_2009' 27 20 28 21 #------------------------------------------------------------------------------- … … 32 25 # Model specific parameters. 33 26 # One or all can be changed each time the run_model script is executed 34 tide = 0. 8# difference between MSL and HAT in metres35 zone = 5 5# specify UTM zone of model36 event_number = 58292# 58280, 64477 the event number or the mux file name27 tide = 0.0 # difference between MSL and HAT in metres 28 zone = 51 # specify UTM zone of model 29 event_number = 27279 # 58280, 64477 the event number or the mux file name 37 30 alpha = 0.1 # smoothing parameter for mesh 38 31 friction=0.01 # manning's friction coefficient 39 32 starttime=0 # start time for simulation 40 finaltime= 65000 # final time for simulation41 42 setup = ' final' # This can be one of three values33 finaltime=1000 # final time for simulation 34 35 setup = 'trial' # This can be one of three values 43 36 # trial - coarsest mesh, fast 44 37 # basic - coarse mesh … … 56 49 #------------------------------------------------------------------------------- 57 50 58 output_comment = [setup, tide, event_number , extent]51 output_comment = [setup, tide, event_number] 59 52 60 53 #------------------------------------------------------------------------------- … … 64 57 # ELEVATION DATA 65 58 # Used in build_elevation.py 66 if extent == None: 67 ascii_grid_filenames = ['topo_grid_notsw', 'fitting_problem_sw_elevation'] 68 point_filenames = [] 69 elif extent == 'sw': 70 # Format for ascii grids, as produced in ArcGIS + a projection file 71 ascii_grid_filenames = ['hob3_pro_extract', 'grid_250m_project', 'ldr_bbay', 72 'ldr_bellerive', 'ldr_brunyislandneck', 73 'ldr_calvertsbeach', 'ldr_carltonbeach', 74 'ldr_connelly', 'ldr_cremorne', 'ldr_dodge1', 75 'ldr_dodge2', 'ldr_kingston', 'ldr_lauderdale', 76 'ldr_linersfarne', 'ldr_oppossumbay', 'ldr_primrose', 77 'ldr_rosny', 'ldr_sandybay', 'ldr_selfspoint', 78 'ldr_sevenmile', 'ldr_southarm1', 'ldr_southarm2', 79 'ldr_sullivanscove', 'ldr_tranmere', 'tasdem_mask'] 80 81 # Format for point is x,y,elevation (with header) 82 point_filenames = ['Topo_1.txt','Topo_2.txt','Topo_3.txt', 83 'tasmania_data.txt', # The data from Tasmania 84 'hydro_data.txt'] # Data from Hydro 59 ascii_grid_filenames = ['town_topo_10m', 60 'cable_250m', 61 'inferred_north', 62 'inferred_south', 63 'north_250m', 64 'other_topo_250m', 65 'south_250m'] 66 point_filenames = ['Broome_coastline.txt','Broome_Bathymetry.txt'] 85 67 86 68 ### Add csv header list to all files in point_filenames … … 98 80 # Used in run_model.py 99 81 # Format for points easting,northing (no header) 100 if extent == 'sw': 101 interior_regions_data = [['aos1.csv', 1500], 102 ['aos2.csv', 1500]] 103 PriorityArea_filename = 'PriorityAreas.csv' 104 elif extent == 'aos1': 105 interior_regions_data = [] 106 PriorityArea_filename = 'PriorityAreas_aos1.csv' 107 elif extent == 'aos2': 108 interior_regions_data = [] 109 PriorityArea_filename = 'PriorityAreas_aos2.csv' 110 elif extent == 'small': 111 interior_regions_data = [['aoi_bruny.csv', 500], 112 ['aoi_S_arms.csv', 500], 113 ['aos1_small.csv', 1500], 114 ['aos2.csv', 1500], 115 ['sw.csv', 30000]] 116 PriorityArea_filename = None 117 elif extent == None: 118 interior_regions_data = [['aos1.csv', 1500], 119 ['aos2.csv', 1500], 120 ['sw.csv', 30000]] 121 PriorityArea_filename = 'PriorityAreas.csv' 82 interior_regions_data = [['Shallow_Water_MH.csv', 1500], 83 ['AoS_MH', 800], 84 ['AoI_MH', 500], 85 ['AoS_north_MH', 800], 86 ['AoI_north_MH', 500], 87 ['bay_coast_MH', 1500], 88 ['bay_coast_small_MH', 800]] 89 PriorityArea_filename = 'PriorityAreas.csv' 122 90 123 91 # LAND - used to set the initial stage/water to be offcoast only 124 92 # Used in run_model.py. Format for points easting,northing (no header) 125 land_initial_conditions_filename = 'initial_conditions .txt'93 land_initial_conditions_filename = 'initial_conditions_mainland.csv' 126 94 127 95 # GAUGES - for creating timeseries at a specific point 128 96 # Used in get_timeseries.py. 129 97 # Format easting,northing,name,elevation (with header) 130 gauges_filename = 'time_of_arrival_hobart.csv' #'TideGaugesPoints.csv' #'tsunamipointsMGA.csv' 131 98 gauges_filename = '' 132 99 # BUILDINGS EXPOSURE - for identifying inundated houses 133 100 # Used in run_building_inundation.py … … 138 105 # Header - easting,northing,id,value 139 106 # Used in get_runup.py 140 images_filename = ' images.csv'107 images_filename = '' 141 108 142 109 # BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively … … 147 114 # Thinned ordering file from Hazard Map (geographic) 148 115 # Format is index,latitude,longitude (with header) 149 urs_order_filename = ' urs_order.csv'116 urs_order_filename = 'thinned_boundary_ordering.csv' 150 117 151 118 # Landward bounding points 152 119 # Format easting,northing (no header) 153 landward_boundary_filename = 'landward_boundary .csv'120 landward_boundary_filename = 'landward_boundary_MH.csv' 154 121 155 122 # MUX input filename. … … 168 135 169 136 # ASCII export grid for Bruny 170 xminBruny = 523900 171 xmaxBruny = 533200 172 yminBruny = 5204300 173 ymaxBruny = 5213100 174 175 # ASCII export grid for South Arms 176 xminSArms = 532497 177 xmaxSArms = 534326 178 yminSArms = 5237028 179 ymaxSArms = 5238465 180 181 # ASCII export grid for Hobart Large 182 xminHobart = 520000 183 xmaxHobart = 545000 184 yminHobart = 5230000 185 ymaxHobart = 5260000 186 187 # ASCII export grid for North West Communities 188 xminNW = 548000 189 xmaxNW = 561000 190 yminNW = 5250000 191 ymaxNW = 5258000 192 193 # ASCII export grid for South Communities 194 xminSouth = 523000 195 xmaxSouth = 533300 196 yminSouth = 5197000 197 ymaxSouth = 5214000 137 ##xminBruny = 523900 138 ##xmaxBruny = 533200 139 ##yminBruny = 5204300 140 ##ymaxBruny = 5213100 198 141 199 142 -
anuga_work/production/Broome_2009/run_model.py
r7159 r7179 89 89 # 4 points equals 5 segments start at N 90 90 boundary_tags={'back': range(num_ocean_segments+1, 91 num_ocean_segments+num_land_points ),91 num_ocean_segments+num_land_points-1), #last segment will be side 92 92 'side': [num_ocean_segments, 93 num_ocean_segments+num_land_points-1, # three side segments total 93 94 num_ocean_segments+num_land_points], 94 95 'ocean': range(num_ocean_segments)} … … 141 142 domain, mean_stage=project.tide, 142 143 time_thinning=1, 143 default_boundary= Bd,144 default_boundary=Dirichlet_boundary([0, 0, 0]), 144 145 boundary_polygon=bounding_polygon_sts, 145 146 use_cache=True,
Note: See TracChangeset
for help on using the changeset viewer.
