Changeset 6672 for anuga_work/production/australia_ph2/strahan

Timestamp:

Mar 31, 2009, 5:39:21 PM (16 years ago)

Author:

myall

Message:

running and exporting scripts

Location:

anuga_work/production/australia_ph2/strahan

Files:

• export_results_max.py (modified) (5 diffs)
• project.py (modified) (4 diffs)

anuga_work/production/australia_ph2/strahan/export_results_max.py

@@ -22 +22 @@
 directory = project.output_folder
 
-time_dir1 = '20090302_104336_run_final_0_27255_mhingee'
-##time_dir2 = ''
-##time_dir3 = ''
-##time_dir4 = ''
-##time_dir5 = ''
-##time_dir6 = ''
-##
-time_dirs = [time_dir1] #time_dir2]# , time_dir4, time_dir5, time_dir6]
+time_dir1 = '20090327_143655_run_final_0_64214_2044_Tb_None_mhingee'
+time_dir2 = '20090327_143836_run_final_0_58337_2044_Tb_None_mhingee'
+time_dir3 = '20090330_095007_run_final_0_68779_2044_Tb_None_mhingee'
+time_dir4 = '20090330_095719_run_final_0_64214_2044_Bf_None_mhingee'
+time_dir5 = '20090330_095845_run_final_0_58337_2044_Bf_None_mhingee'
+time_dir6 = '20090330_143221_run_final_0_68779_2044_Bf_None_mhingee'
+
+time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
   
 
@@ -36 +36 @@
 
 #timestep = None    # None means no timestep!
-timestep = 0
+timestep = None
 
 ######
@@ -56 +56 @@
 
 # one or more key strings from var_equations above
-var = ['elevation','stage']
+var = ['stage','speed']
 
 ######
@@ -91 +91 @@
             names = [name1] #, name2]#, name3]
       
-           # asc_name = []   
+            asc_name = []   
                
             for name in names:
@@ -112 +112 @@
                             format = 'asc')
 
-##                asc_name.append(outname + '.asc')
-##    
-##            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
-##
-##            print 'max asc outname ', maxasc_outname
-##            print 'asc_name ', str(asc_name)
-##            
-##            MaxAsc(maxasc_outname, asc_name)
+                asc_name.append(outname + '.asc')
+    
+            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
+
+            print 'max asc outname ', maxasc_outname
+            print 'asc_name ', str(asc_name)
+            
+            MaxAsc(maxasc_outname, asc_name)

anuga_work/production/australia_ph2/strahan/project.py

@@ -28 +28 @@
 tide = 0                # difference between MSL and HAT
 central_meridian = 147.0 # Central meridian for projection (optional)
-event_number =     # the event number or the mux file name
-# main events: 64214 (sandwich 1.4e-5), 58337 (puysegur 3e-6), 
+event_number = 68779    # the event number or the mux file name
+# main events: 64214 (sandwich 1.4e-5), 58337 (puysegur 3e-6), 68779 (sumatra 7.53-7)
 # others: 27295 (java 3e-7), 31817 (kermadec 5e-7), 46559 (midAmerica 4.2e-7)
 # 51460 (newHebrides 7e-7), 57477 (peru 5.4e-7), 65322 (shetland 3.6e-7),
-# 68779 (sumatra 7.53-7)
+
 alpha = 0.1             # smoothing parameter for mesh
 friction=0.01           # manning's friction coefficient
@@ -40 +40 @@
 # index is only used when wave = Tb
 index = 2044            # index from the PTHA - Y2000 is 0.60374m
-wave = 'Tb'             # Bf (sts wave) Tb (index wave)
+wave = 'Bf'             # Bf (sts wave) Tb (index wave)
 
 setup = 'final'         # This can be one of three values
@@ -77 +77 @@
 
 # Format for point is x,y,elevation (with header)
-point_filenames = ['grid250m_pts.txt',
-                   'coast_250.txt'] # 250m grid 2005
+point_filenames = ['grid_250m.txt'] # 250m grid 2005
 
 ### Add csv header list to all files in point_filenames
@@ -125 +124 @@
 # Thinned ordering file from Hazard Map (geographic)
 # Format is index,latitude,longitude (with header)
-urs_order_filename = 'urs_order_austtg.csv'
+urs_order_filename = 'urs_order.csv'
 
 # Landward bounding points