Changeset 6796 for anuga_work/production/australia_ph2/ceduna
- Timestamp:
- Apr 15, 2009, 8:51:00 AM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/australia_ph2/ceduna/project.py
r6726 r6796 8 8 from time import localtime, strftime, gmtime 9 9 from os.path import join, exists 10 10 from anuga.lib.add_csv_header.add_csv_header import add_csv_header 11 11 12 12 #------------------------------------------------------------------------------- … … 18 18 state = 'australia_ph2' 19 19 scenario_name = 'ceduna' 20 scenario_folder = scenario_name21 20 22 21 #------------------------------------------------------------------------------- … … 26 25 # Model specific parameters. 27 26 # One or all can be changed each time the run_model script is executed 28 tide = 0 # difference between MSL and HAT 29 ##zone = 53 # specify zone of model 30 central_meridian = 133.0 # nonstandard projection 31 27 28 central_meridian = 133 # Central meridian for projection (optional) 29 zone = None 32 30 import sys 33 31 if len(sys.argv) > 1: 34 32 event_number = int(sys.argv[1]) 35 33 else: 36 34 event_number = 27347 # the event number or the mux file name 37 35 38 36 event_number_list = [27347, 64448, 58331] # To piggy back multiple events … … 43 41 # 57483 (peru), 65304 (shetland), 68792 (sumatra) 44 42 43 tide = 0 # difference between MSL and HAT 45 44 alpha = 0.1 # smoothing parameter for mesh 46 45 friction=0.01 # manning's friction coefficient 47 46 starttime=0 # start time for simulation 48 finaltime=1000 # final time for simulation 49 50 # index is only used when wave = Tb 51 index = 1889 # index from the PTHA 52 wave = 'Bf' # Bf (sts wave) Tb (index wave) 53 54 setup = 'trial' # This can be one of three values 47 finaltime=60000 #60000 # final time for simulation 48 setup = 'final' # This can be one of three values 55 49 # trial - coarsest mesh, fast 56 50 # basic - coarse mesh 57 51 # final - fine mesh, slowest 52 53 # index is only used when wave = Tb 54 index = 1884 # index from the PTHA - Y2000 0.257m 55 wave = 'Tb' # Bf (sts wave) Tb (index wave) 56 58 57 59 58 internal_polygon = False … … 73 72 internal_poly_comment = 'internal' 74 73 else: 75 internal_poly_comment = None74 internal_poly_comment = '' 76 75 77 76 output_comment = [setup, tide, event_number, index, wave, internal_poly_comment] … … 84 83 # Used in build_elevation.py 85 84 # Format for ascii grids, as produced in ArcGIS + a projection file 86 ascii_grid_filenames = [] # 250m grid 2005 87 88 # Format for point is x,y,elevation (with header) 89 point_filenames = ['grid_250m.txt'] # 250m grid 2005 85 ascii_grid_filenames = ['grid_250m_simple'] 86 87 # Format for points is x,y,elevation (with header) 88 point_filenames = [] 89 ##point_filenames = ['brisbane_250m.txt', 90 ## 'GBR_250m.txt', 91 ## 'Sydney_250m.txt'] # 250m grid 2005 90 92 91 93 ### Add csv header list to all files in point_filenames … … 94 96 ## add_csv_header(join(topographies_folder, f), headerlist) 95 97 96 98 # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh 97 99 # Used in build_elevation.py 98 # Format for points easting,northing (no header)99 bounding_polygon_filename = 'bounding_polygon .csv'100 bounding_polygon_maxarea = 1 00000100 # Format for points: easting,northing (no header) 101 bounding_polygon_filename = 'bounding_polygon_simple.csv' 102 bounding_polygon_maxarea = 125000 101 103 102 104 # INTERIOR REGIONS - for designing the mesh … … 125 127 # Used in run_building_inundation.py 126 128 # Format latitude,longitude etc (geographic) 127 building_exposure_filename = ' busselton_res_clip.csv' # from NEXIS129 building_exposure_filename = '.csv' # from NEXIS 128 130 129 131 # BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively … … 134 136 # Thinned ordering file from Hazard Map (geographic) 135 137 # Format is index,latitude,longitude (with header) 136 urs_order_filename = 'urs_order .csv'138 urs_order_filename = 'urs_order_simple.csv' 137 139 138 140 # Landward bounding points 139 141 # Format easting,northing (no header) 140 landward_boundary_filename = 'landward_boundary .csv'142 landward_boundary_filename = 'landward_boundary_simple.csv' 141 143 142 144 # MUX input filename. 143 145 # If a meta-file from EventSelection is used, set 'multi-mux' to True. 144 146 # If a single MUX stem filename (*.grd) is used, set 'multi-mux' to False. 145 ##mux_input_filename = event_number # to be found in event_folder 146 # (ie boundaries/event_number/) 147 ##mux_input_filename = 'Java-0016-z.grd' 147 148 ##multi_mux = False 148 149 mux_input_filename = 'event.list' 149 150 multi_mux = True 151 152 # Specify if share cache is to be used 153 # Whatever is specified here will be relative to INUNDATION_HOME/.cache 154 # If nothing is specified, local cache will be used. 155 cachedir = '.python_cache_phII' 150 156 151 157 … … 175 181 # Output filename for elevation 176 182 # this is a combination of all the data generated in build_elevation.py 177 combined_elevation_basename = scenario_name + ' _combined_elevation'183 combined_elevation_basename = scenario_name + 'simple_combined_elevation' 178 184 179 185 #------------------------------------------------------------------------------- … … 190 196 home = join(os.getenv(ENV_INUNDATIONHOME), 'data') # Absolute path for data folder 191 197 muxhome = os.getenv(ENV_MUXHOME) 198 199 # Create absolute pathname for cache directory 200 # and change caching to use it 201 if 'cachedir' in dir(): 202 cachedir = join(os.getenv(ENV_INUNDATIONHOME), '.cache', cachedir) 203 from anuga.caching import caching 204 caching.set_option('cachedir', cachedir) 192 205 193 206 # check various directories/files that must exist 194 anuga_folder = join(home, state, scenario_ folder, 'anuga')207 anuga_folder = join(home, state, scenario_name, 'anuga') 195 208 topographies_folder = join(anuga_folder, 'topographies') 196 209 polygons_folder = join(anuga_folder, 'polygons') … … 229 242 event_sts = join(event_folder, scenario_name) 230 243 231 232 233 244 # The absolute pathname for the output folder names 234 245 # Used for build_elevation.py
Note: See TracChangeset
for help on using the changeset viewer.