Changeset 6773

Timestamp:

Apr 8, 2009, 4:47:05 PM (15 years ago)

Author:

kristy

Message:

 

File:

• anuga_work/production/australia_ph2/sydney/project.py (modified) (12 diffs)

anuga_work/production/australia_ph2/sydney/project.py

@@ -8 +8 @@
 from time import localtime, strftime, gmtime
 from os.path import join, exists
-from os import sep
-
+from anuga.lib.add_csv_header.add_csv_header import add_csv_header
 
 #-------------------------------------------------------------------------------
@@ -19 +18 @@
 state = 'australia_ph2'
 scenario_name = 'sydney'
-scenario_folder = scenario_name
 
 #-------------------------------------------------------------------------------
@@ -27 +25 @@
 # Model specific parameters.
 # One or all can be changed each time the run_model script is executed
+
+central_meridian = 151 # Central meridian for projection (optional)
+zone = None
+import sys
+if len(sys.argv) > 1:
+    event_number = int(sys.argv[1])
+else:    
+    event_number = 58152    # the event number or the mux file name
+
+event_number_list = [58152, 58222, 58326] # To piggy back multiple events
+
+# events: 27346 (java 2.1e-6), 58334 (puysegur 8.8e-7), 64322 (sandwich 8.2e-6), 
+# other events: 51466 (newHebrides 7.2e-7), 57537 (peru 4.7e-7), 46591 (midAmerica 1.3e-6),
+# 68791 (sumatra 7.5e-7), 11100 (andaman 4.9e-7), 65368 (shetland 2.5e-6), 7867 (altiplano 1.0e-6)
 tide = 0                # difference between MSL and HAT
-zone = 56              # specify zone of model
-event_number = 58152    # the event number or the mux file name
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
 starttime=0             # start time for simulation
-finaltime=80000          # final time for simulation
-
-# index is only used when wave = Tb
-index = 2956            # index from the PTHA
-wave = 'Tb'             # Bf (sts wave) Tb (index wave)
-
-setup = 'final'         # This can be one of three values
+finaltime=1000   #60000 # final time for simulation
+setup = 'trial'         # This can be one of three values
                         #    trial - coarsest mesh, fast
                         #    basic - coarse mesh
                         #    final - fine mesh, slowest
+
+# index is only used when wave = Tb
+index = 2938            # index from the PTHA - Y2000 0.257m
+wave = 'Tb'             # Bf (sts wave) Tb (index wave)
+
+
+internal_polygon = False
 
 #-------------------------------------------------------------------------------
@@ -55 +67 @@
 #-------------------------------------------------------------------------------
 
-output_comment = [setup, tide, event_number, index, wave]
+if internal_polygon:
+    internal_poly_comment = 'internal'
+else:
+    internal_poly_comment = ''
+
+output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
 
 #-------------------------------------------------------------------------------
@@ -64 +81 @@
 # Used in build_elevation.py
 # Format for ascii grids, as produced in ArcGIS + a projection file
-ascii_grid_filenames = [] # 250m grid 2005
-
-# Format for point is x,y,elevation (with header)
-point_filenames = ['sydney_elevation.txt',
-                   'melbourne_elevation.txt'] # 250m grid 2005
+ascii_grid_filenames = ['grid_250m_simple']
+
+# Format for points is x,y,elevation (with header)
+point_filenames = []
+##point_filenames = ['sydney_elevation.txt',
+##                   'melbourne_elevation.txt'] # 250m grid 2005
 
 ### Add csv header list to all files in point_filenames
@@ -75 +93 @@
 ##    add_csv_header(join(topographies_folder, f), headerlist)
 
-    # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
+# BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
 # Used in build_elevation.py
-# Format for points easting,northing (no header)
-bounding_polygon_filename = 'bounding_polygon.csv'
-bounding_polygon_maxarea = 100000
+# Format for points: easting,northing (no header)
+bounding_polygon_filename = 'bounding_polygon_simple.csv'
+bounding_polygon_maxarea = 125000
 
 # INTERIOR REGIONS -  for designing the mesh
@@ -86 +104 @@
 interior_regions_data = []
 
+# add an internal polygon to force different mesh generation
+# used to test for discretisation error when building elevation
+# make sure file is in same folder as interior regions and bouding polygon;
+# format is same (2 column .csv; easting, northing; no header)
+if internal_polygon:
+    interior_regions_data.append(['internal_polygon.csv',
+                                  bounding_polygon_maxarea])
+
 # LAND - used to set the initial stage/water to be offcoast only
 # Used in run_model.py.  Format for points easting,northing (no header)
@@ -93 +119 @@
 # Used in get_timeseries.py.  
 # Format easting,northing,name,elevation (with header)
-gauges_filename = 'MH_gauges_10m2.csv' #'gauges.csv'
+gauges_filename = 'gauges.csv'
 
 # BUILDINGS EXPOSURE - for identifying inundated houses
 # Used in run_building_inundation.py
 # Format latitude,longitude etc (geographic)
-building_exposure_filename = 'sydney_res_clip.csv' # from NEXIS
+building_exposure_filename = '.csv' # from NEXIS
 
 # BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
@@ -107 +133 @@
 # Thinned ordering file from Hazard Map (geographic)
 # Format is index,latitude,longitude (with header)
-urs_order_filename = 'urs_order.csv'
+urs_order_filename = 'urs_order_simple.csv'
 
 # Landward bounding points
 # Format easting,northing (no header)
-landward_boundary_filename = 'landward_boundary.csv'
+landward_boundary_filename = 'landward_boundary_simple.csv'
 
 # MUX input filename.
 # If a meta-file from EventSelection is used, set 'multi-mux' to True.
 # If a single MUX stem filename (*.grd) is used, set 'multi-mux' to False.
-##mux_input_filename = event_number # to be found in event_folder
-                                    # (ie boundaries/event_number/)
+##mux_input_filename = 'Java-0016-z.grd'
 ##multi_mux = False
 mux_input_filename = 'event.list'
 multi_mux = True
 
-
+# Specify if share cache is to be used
+# Whatever is specified here will be relative to INUNDATION_HOME/.cache
+# If nothing is specified, local cache will be used.
+cachedir = '.python_cache_phII'
 
 
@@ -150 +178 @@
 # Output filename for elevation
 # this is a combination of all the data generated in build_elevation.py
-combined_elevation_basename = scenario_name + '_combined_elevation'
+combined_elevation_basename = scenario_name + 'simple_combined_elevation'
 
 #-------------------------------------------------------------------------------
@@ -165 +193 @@
 home = join(os.getenv(ENV_INUNDATIONHOME), 'data') # Absolute path for data folder
 muxhome = os.getenv(ENV_MUXHOME)
+
+# Create absolute pathname for cache directory
+# and change caching to use it
+if 'cachedir' in dir():
+    cachedir = join(os.getenv(ENV_INUNDATIONHOME), '.cache', cachedir)
+    from anuga.caching import caching
+    caching.set_option('cachedir', cachedir) 
     
 # check various directories/files that must exist
-anuga_folder = join(home, state, scenario_folder, 'anuga')
+anuga_folder = join(home, state, scenario_name, 'anuga')
 topographies_folder = join(anuga_folder, 'topographies')
 polygons_folder = join(anuga_folder, 'polygons')
@@ -223 +258 @@
 mux_input = join(event_folder, mux_input_filename)
 
-boundary_csv = event_folder + sep + 'sts_gauge_' + str(index) +'.csv'
+# sts gauge with desired index number - used for wave 'Tb'
+boundary_csv = join(event_folder, 'sts_gauge_' + str(index) + '.csv')
+
+