Changeset 6769 for anuga_work/production/australia_ph2/lockhart_river

Timestamp:

Apr 8, 2009, 1:18:00 PM (16 years ago)

Author:

myall

Message:

running GBR models without LADS data

File:

• anuga_work/production/australia_ph2/lockhart_river/project.py (modified) (10 diffs)

anuga_work/production/australia_ph2/lockhart_river/project.py

@@ -8 +8 @@
 from time import localtime, strftime, gmtime
 from os.path import join, exists
-
+from anuga.lib.add_csv_header.add_csv_header import add_csv_header
 
 #-------------------------------------------------------------------------------
@@ -18 +18 @@
 state = 'australia_ph2'
 scenario_name = 'lockhart_river'
-scenario_folder = scenario_name
 
 #-------------------------------------------------------------------------------
@@ -26 +25 @@
 # Model specific parameters.
 # One or all can be changed each time the run_model script is executed
-tide = 0                # difference between MSL and HAT
-##zone = 55               # specify zone of model
-central_meridian = 143.5 # nonstandard projection
-
-
+
+central_meridian = 143.5 # Central meridian for projection (optional)
+zone = None
 import sys
 if len(sys.argv) > 1:
@@ -36 +33 @@
 else:    
     event_number = 51316    # the event number or the mux file name
-    
+
 event_number_list = [51316, 62786] # To piggy back multiple events
-                                          # using run_multiple_events.py
-
+
+tide = 0                # difference between MSL and HAT
 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 = 1889            # index from the PTHA
-wave = 'Tb'             # Bf (sts wave) Tb (index wave)
-
+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 = 1994            # index from the PTHA - Y2000 0.257m
+wave = 'Tb'             # Bf (sts wave) Tb (index wave)
+
 
 internal_polygon = False
@@ -70 +67 @@
     internal_poly_comment = 'internal'
 else:
-    internal_poly_comment = None
+    internal_poly_comment = ''
 
 output_comment = [setup, tide, event_number, index, wave, internal_poly_comment]
@@ -83 +80 @@
 ascii_grid_filenames = ['grid250m_pro'] # 250m grid 2005
 
-# Format for point is x,y,elevation (with header)
-##point_filenames = []
-point_filenames = ['ahs100_0.txt',
-                   'ahs100_1.txt',
-                   'ahs100_2.txt',
-                   'ahs100_3.txt',
-                   'ahs100_4.txt',
-                   'ahs100_5.txt',
-                   'ahs100_6.txt',
-                   'ahs100_7.txt',
-                   'ahs100_8.txt',
-                   'ahs100_9.txt',
-                   'ahs100_10.txt',
-                   'ahs100_11.txt',
-                   'ahs100_12.txt',
-                   'ahs100_13.txt'] # LADS data
+# Format for points is x,y,elevation (with header)
+point_filenames = []
+##point_filenames = ['ahs100_0.txt',
+##                   'ahs100_1.txt',
+##                   'ahs100_2.txt',
+##                   'ahs100_3.txt',
+##                   'ahs100_4.txt',
+##                   'ahs100_5.txt',
+##                   'ahs100_6.txt',
+##                   'ahs100_7.txt',
+##                   'ahs100_8.txt',
+##                   'ahs100_9.txt',
+##                   'ahs100_10.txt',
+##                   'ahs100_11.txt',
+##                   'ahs100_12.txt',
+##                   'ahs100_13.txt'] # LADS data
+
+### Add csv header list to all files in point_filenames
+##headerlist = ['x', 'y', 'elevation'] 
+##for f in point_filenames:
+##    add_csv_header(join(topographies_folder, f), headerlist)
 
 # BOUNDING POLYGON - for data clipping and estimate of triangles in mesh
 # Used in build_elevation.py
-# Format for points easting,northing (no header)
+# Format for points: easting,northing (no header)
 bounding_polygon_filename = 'bounding_polygon.csv'
 bounding_polygon_maxarea = 125000
@@ -131 +133 @@
 # Used in run_building_inundation.py
 # Format latitude,longitude etc (geographic)
-building_exposure_filename = 'busselton_res_clip.csv' # from NEXIS
+building_exposure_filename = '.csv' # from NEXIS
 
 # BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
@@ -149 +151 @@
 # 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'
 
 
@@ -196 +202 @@
 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')
@@ -235 +248 @@
 event_sts = join(event_folder, scenario_name)
 
-
-
 # The absolute pathname for the output folder names
 # Used for build_elevation.py