Changeset 7089

Timestamp:

May 26, 2009, 3:52:17 PM (16 years ago)

Author:

Leharne

Message:

updates to gold coast inundation scenario

Location:

anuga_work/production/gold_coast_2009

Files:

• Arc_asc2raster_GDA94z56.py (modified) (4 diffs)
• export_results_max.py (modified) (4 diffs)
• project.py (modified) (3 diffs)

anuga_work/production/gold_coast_2009/Arc_asc2raster_GDA94z56.py

@@ -25 +25 @@
 output_dir = "anuga\\outputs\\"
 
-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_dir1 = '20090508_150215_run_final_0_51469_lfountai'
+time_dir2 = '20090511_165539_run_final_0_50863_lfountai'
+time_dir3 = '20090511_161526_run_final_1.1_51469_kvanputt'
+time_dir4 = '20090518_154710_run_final_0_50994_lfountai'
+##time_dir5 = '20090519_160510_run_final_1.1_50994_lfountai'
+##time_dir6 = '20090521_220101_run_final_1.1_50863_kvanputt'
+##time_dir7 = '20090522_164526_run_final_0_51392_lfountai'
+##time_dir8 = '20090522_164640_run_final_1.1_51392_lfountai'
+##time_dir9 = '20090522_164948_run_final_0_51423_lfountai'
+##time_dir10 = '20090522_165600_run_final_1.1_51423_lfountai'
 
-time_dirs = [time_dir1] #[time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
+time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4] #, time_dir5, time_dir6]
 
 for time_dir in time_dirs:
@@ -40 +44 @@
     raster_gbd = folder + 'raster.gdb'
     #contour = raster_gbd + '\\contour_dep'
-    land = scenario_dir + "ArcGIS\\gold_coast.gdb\\Coast_250"
+##    land = scenario_dir + "ArcGIS\\gold_coast.gdb\\Coast_250"
 ##    ocean = scenario_dir + "map_work\\Perth.gdb\\Outlines\\initial_condition_ocean"
     
@@ -51 +55 @@
     
     #replication dictionary
-    replicate = (('gold_coast', ''), ('_time_21960_0', 'b'), ('_time_43920_0', 'c'),
-                 ('time_65880_0', 'd'), ('_', ''), ('max','M_'), ('depth','_dep_'),
+    replicate = (('gold_coast', ''), ('_time_29220_0', 'b'), ('_time_58440_0', 'c'),
+                 ('_', ''), ('max','M_'), ('depth','_dep_'),
                  ('speed', '_spe_'), ('elevation', '_ele_'), ('stage','_stage'))
 
     generate_filename = []
-    input_ascii = glob.glob(folder + '*depth_max.asc')
+    input_ascii = glob.glob(folder + '*_max.asc')
     print time_dir
 
@@ -80 +84 @@
 
     
-        output_extract = output_DEM + 'E'
-        print 'Output Extract ',output_extract
-    
-        print 'Process: Extract by Mask'
-        gp.ExtractByMask_sa(output_DEM, land, output_extract)
-
+##        output_extract = output_DEM + 'E'
+##        print 'Output Extract ',output_extract
+##    
+##        print 'Process: Extract by Mask'
+##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
+##
        

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_dir1] #time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]   
+time_dir1 = '20090521_220101_run_final_1.1_50863_kvanputt'
+time_dir2 = '20090522_164526_run_final_0_51392_lfountai'
+time_dir3 = '20090522_164640_run_final_1.1_51392_lfountai'
+time_dir4 = '20090522_164948_run_final_0_51423_lfountai'
+time_dir5 = '20090522_165600_run_final_1.1_51423_lfountai
+
+time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]   
 
 ##cellsize = 250
@@ -41 +41 @@
 
 #area = ['Bunbury', 'Busselton']  # strings must match keys in var_equations below
-area = ['Dune']      # 'All' means no special areas - the whole thing
+area = ['All']      # 'All' means no special areas - the whole thing
 
 ######
@@ -54 +54 @@
 
 # one or more key strings from var_equations above
-var = ['depth']
+var = ['speed', 'depth']
 
 ######
@@ -84 +84 @@
 
             name1 = join(directory, time_dir, project.scenario_name)
-            name2 = join(directory, time_dir, project.scenario_name)+'_time_21960_0'
-            name3 = join(directory, time_dir, project.scenario_name)+'_time_43920_0'
-            #name4 = join(directory, time_dir, project.scenario_name)+'_time_65880_0'
-
+            name2 = join(directory, time_dir, project.scenario_name)+'_time_29220_0'
+            name3 = join(directory, time_dir, project.scenario_name)+'_time_58440_0'
+            
             names = [name1, name2, name3]
       

anuga_work/production/gold_coast_2009/project.py

@@ -34 +34 @@
     event_number = int(sys.argv[1])
 else:    
-    event_number = 51469    # the event number or the mux file name
+    event_number = 51423    # the event number or the mux file name
 
 event_number_list = [51469, 50863] # To piggy back multiple events
 
-tide = 0 
+tide = 1.1
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
@@ -106 +106 @@
 building_exposure_filename = 'gold_coast_res_clip.csv' # from NEXIS
 
+# AREA OF IMAGES - Extent of each image to find out highest runup
+# Headerr: easting, northing, id, value
+# Used in get_runup.py
+images_filename = 'images.csv'
+
 # BOUNDING POLYGON - used in build_boundary.py and run_model.py respectively
 # NOTE: when files are put together the points must be in sequence
 # For ease go clockwise!
 # Check the run_model.py for boundary_tags
+
 
 # Thinned ordering file from Hazard Map (geographic)
@@ -244 +250 @@
     building_exposure = join(gauges_folder, building_exposure_filename)
 
+# The absolute pathname for the image file
+# Used for get_runup.py
+if images_filename:
+    images = join(polygons_folder, images_filename)
+
 # full path to where MUX files (or meta-files) live
 mux_input = join(event_folder, mux_input_filename)