Changeset 7004 for anuga_work/production/gold_coast_2009

Timestamp:

May 8, 2009, 4:25:40 PM (16 years ago)

Author:

Leharne

Message:

Updates to Gold Coast tsunami inundation scenario. Two out of three final events selected and area of interest redefined.

Location:

anuga_work/production/gold_coast_2009

Files:

• Arc_asc2raster_GDA94z56.py (added)
• build_urs_boundary.py (modified) (2 diffs)
• export_results_max.py (modified) (1 diff)
• get_timeseries.py (modified) (1 diff)
• project.py (modified) (5 diffs)
• run_multiple_events.py (added)

anuga_work/production/gold_coast_2009/build_urs_boundary.py

@@ -140 +140 @@
                 weights=mux_weights,
                 zone=project.zone,
-                mean_stage=project.tide,
                 verbose=True)
     else:                           # a single mux stem file, assume 1.0 weight
@@ -158 +157 @@
                 ordering_filename=order_filename,
                 weights=mux_weights,
-                mean_stage=project.tide,
                 verbose=True)
 

anuga_work/production/gold_coast_2009/export_results_max.py

@@ -22 +22 @@
 directory = project.output_folder
 
-#time_dir1 = '20090327_135943_run_final_0_51469_lfountai'
-time_dir2 = '20090327_153206_run_final_0_51348_lfountai'
-time_dir3 = '20090327_153231_run_final_0_51253_lfountai'
-time_dir4 = '20090327_153304_run_final_0_50863_lfountai'
-time_dir5 = '20090327_153403_run_final_0_51423_lfountai'
-time_dir6 = '20090327_153602_run_final_0_62833_lfountai'
-time_dirs = [time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]   
+time_dir1 = '20090327_135943_run_final_0_51469_lfountai'
+##time_dir2 = '20090327_153206_run_final_0_51348_lfountai'
+##time_dir3 = '20090327_153231_run_final_0_51253_lfountai'
+##time_dir4 = '20090327_153304_run_final_0_50863_lfountai'
+##time_dir5 = '20090327_153403_run_final_0_51423_lfountai'
+##time_dir6 = '20090327_153602_run_final_0_62833_lfountai' 
+time_dirs = [time_dir1] #time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]   
 
 ##cellsize = 250

anuga_work/production/gold_coast_2009/get_timeseries.py

@@ -19 +19 @@
 directory = project.output_folder
 
-#time_dir1 = '20090327_135943_run_final_0_51469_lfountai'
-time_dir2 = '20090327_153206_run_final_0_51348_lfountai'
-time_dir3 = '20090327_153231_run_final_0_51253_lfountai'
-time_dir4 = '20090327_153304_run_final_0_50863_lfountai'
+time_dir1 = '20090327_135943_run_final_0_51469_lfountai'
+##time_dir2 = '20090327_153206_run_final_0_51348_lfountai'
+##time_dir3 = '20090327_153231_run_final_0_51253_lfountai'
+##time_dir4 = '20090327_153304_run_final_0_50863_lfountai'
 time_dir5 = '20090327_153403_run_final_0_51423_lfountai'
 time_dir6 = '20090327_153602_run_final_0_62833_lfountai'
-time_dirs = [time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
+time_dirs = [time_dir1, time_dir5, time_dir6] #time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
 
 for time_dir in time_dirs:

anuga_work/production/gold_coast_2009/project.py

@@ -27 +27 @@
 # One or all can be changed each time the run_model script is executed
 
-model_description = 'small_v2'
-
-tide = 0                # difference between MSL and HAT in metres
-zone = 56               # specify zone of model
-##event_number = 51469    # the event number or the mux file name with:
-                          # 2.32m waveheight
-##event_number = 51423    # 1.5m waveheight
-##event_number= 51348    # 1.0m waveheight 
-##event_number = 51253    # 0.5m wavheight
-##event_number = 50863    # 0.3m wavheight from New Hebridies
-event_number = 62833    # 0.3m waveheight from SE Solomons
+               # difference between MSL and HAT in metres
+central_meridian = None # Central meridian for projection (optional)
+zone = 56
+import sys
+if len(sys.argv) > 1:
+    event_number = int(sys.argv[1])
+else:    
+    event_number = 51469    # the event number or the mux file name
+
+event_number_list = [51469, 50863] # To piggy back multiple events
+
+tide = 0 
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
@@ -59 +60 @@
 #-------------------------------------------------------------------------------
 
-output_comment = [setup, tide, model_description, event_number]
+output_comment = [setup, tide, event_number]
 
 #-------------------------------------------------------------------------------
@@ -82 +83 @@
 # Used in build_elevation.py
 # Format for points easting,northing (no header)
-bounding_polygon_filename = 'bounding_polygon_sml_v2.csv'
+bounding_polygon_filename = 'bounding_polygon_sml.csv'
 bounding_polygon_maxarea = 125000
 
@@ -88 +89 @@
 # Used in run_model.py
 # Format for points easting,northing (no header)                    
-interior_regions_data = [['buffer_3km_simplify_sml_v2.csv', 500]]
+interior_regions_data = [['area_of_interest.csv', 500],
+                         ['intermediate.csv', 25000]]
 
 # LAND - used to set the initial stage/water to be offcoast only
 # Used in run_model.py.  Format for points easting,northing (no header)
-land_initial_conditions_filename = [['initial_conditions_sml_v2.csv', 0]]
+land_initial_conditions_filename = [['initial_conditions_sml.csv', 0]]
 
 # GAUGES - for creating timeseries at a specific point
@@ -111 +113 @@
 # Thinned ordering file from Hazard Map (geographic)
 # Format is index,latitude,longitude (with header)
-urs_order_filename = 'urs_order_sml_v2.csv'
+urs_order_filename = 'urs_order_sml.csv'
 
 # Landward bounding points
 # Format easting,northing (no header)
-landward_boundary_filename = 'landward_boundary_sml_v2.csv'
+landward_boundary_filename = 'landward_boundary_sml.csv'
 
 # MUX input filename.