Changeset 6797 for anuga_work/production/australia_ph2/dampier/project.py

Timestamp:

Apr 15, 2009, 8:51:36 AM (15 years ago)

Author:

kristy

Message:

 

File:

• anuga_work/production/australia_ph2/dampier/project.py (modified) (8 diffs)

anuga_work/production/australia_ph2/dampier/project.py

@@ -29 +29 @@
 tide = 0                # difference between MSL and HAT
 zone = 50               # specify zone of model
-event_number = 70385    # the event number or the mux file name
+central_meridian = None
+
+import sys
+if len(sys.argv) > 1:
+    event_number = int(sys.argv[1])
+else:    
+    event_number = 70844    # the event number or the mux file name
+
+event_number_list = [70844, 27338, 70337] # To piggy back multiple events
+
 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
+finaltime=1000          # final time for simulation
+
+# for sensetivity study
+model = 'simple' #60km 90km 120km all
 
 # index is only used when wave = Tb
-index = 1203            # index from the PTHA
-wave = 'Bf'             # Bf (sts wave) Tb (index wave)
-
-setup = 'final'         # This can be one of three values
+index = 1307            # index from the PTHA
+wave = 'Tb'             # Bf (sts wave) Tb (index wave)
+
+setup = 'trial'         # This can be one of three values
                         #    trial - coarsest mesh, fast
                         #    basic - coarse mesh
@@ -55 +67 @@
 #-------------------------------------------------------------------------------
 
-output_comment = [setup, tide, event_number, index, wave]
+output_comment = [setup, tide, event_number, wave, model]
 
 #-------------------------------------------------------------------------------
@@ -64 +76 @@
 # Used in build_elevation.py
 # Format for ascii grids, as produced in ArcGIS + a projection file
-ascii_grid_filenames = [] # 250m grid 2005
+ascii_grid_filenames = ['grid250m_simple'] # 250m grid 2005
 
 # Format for point is x,y,elevation (with header)
-point_filenames = ['grid250m.txt'] # 250m grid 2005
+point_filenames = [] # 250m grid 2005
 
 ### Add csv header list to all files in point_filenames
@@ -77 +89 @@
 # Used in build_elevation.py
 # Format for points easting,northing (no header)
-bounding_polygon_filename = 'bounding_polygon.csv'
+bounding_polygon_filename = 'poly_all_' + model + '.csv'
 bounding_polygon_maxarea = 100000
 
@@ -106 +118 @@
 # 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_' + model + '.csv'
 
 # Landward bounding points
 # Format easting,northing (no header)
-landward_boundary_filename = 'landward_boundary.csv'
+landward_boundary_filename = 'landward_boundary_' + model + '.csv'
 
 # MUX input filename.
@@ -121 +133 @@
 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'
 
 ################################################################################
@@ -164 +178 @@
 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')
@@ -201 +222 @@
 
 # The absolute pathname for the .sts file, generated in build_boundary.py
-event_sts = join(event_folder, scenario_name)
+event_sts = join(event_folder, scenario_name)+ '_' +model
 
 # The absolute pathname for the output folder names