Changeset 5952

Timestamp:

Nov 12, 2008, 3:35:44 PM (16 years ago)

Author:

kristy

Message:

Running smoothly, for all var and area

File:

• anuga_work/production/perth/export_results_max.py (modified) (5 diffs)

anuga_work/production/perth/export_results_max.py

@@ -20 +20 @@
 directory = project.output_dir
 
-time_dir1 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'
+#time_dir1 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'
 #time_dir1 = '20080924_123626_run_final_0_27283_250m_all_kvanputt' # This uses the 250m bathymetry with all interior polygons
 #time_dir2 = '20080912_154716_run_final_0_27283_alpha0.1_kvanputt' # This uses original bathyemetry data
 #time_dir3 = '20080924_123601_run_final_0_27283_250m_none_kvanputt' # This uses the 250m bathymetry without any interior polygons
-time_dirs = [time_dir1] #, time_dir2]#, time_dir3]
+#time_dirs = [time_dir1] #, time_dir2]#, time_dir3]
 
-#cellsize = 20
-cellsize = 250
+time_dir1 = '20081031_133353_run_final_0.6_68693_alpha0.1_kvanputt'
+time_dir2 = '20081031_133511_run_final_0_68693_alpha0.1_kvanputt'
+time_dir3 = '20081031_133624_run_final_0_27255_alpha0.1_kvanputt'
+time_dir4 = '20081031_133735_run_final_0.6_27255_alpha0.1_kvanputt'
+time_dir5 = '20081031_133841_run_final_0_27283_alpha0.1_kvanputt'
+time_dir6 = '20081031_133925_run_final_0.6_27283_alpha0.1_kvanputt'
+
+time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
+
+  
+
+cellsize = 20
+#cellsize = 250
 
 timestep = None    # None means no timestep!
@@ -36 +47 @@
 ######
 
-#area = ['Geordie', 'Sorrento', 'Fremantle', 'Rockingham']  # strings must match keys in var_equations below
-area = ['All']      # 'All' means no special areas - the whole thing
+area = ['Rockingham', 'Sorrento', 'Fremantle','Geordie' ]  # strings must match keys in var_equations below
+#area = ['All']      # 'All' means no special areas - the whole thing
 
 ######
@@ -50 +61 @@
 
 # one or more key strings from var_equations above
-var = ['stage', 'elevation']
+var = ['depth', 'speed']
 
 ######
 # Start running the various conversions we require.
 ######
+
 for which_var in var:
     if which_var not in var_equations:
@@ -60 +72 @@
         break
 
-    asc_name = []
+    for which_area in area:
+        if which_area == 'All':
+            easting_min = None
+            easting_max = None
+            northing_min = None
+            northing_max = None
+        else:
+            try:
+                easting_min = eval('project.xmin%s' % which_area)
+                easting_max = eval('project.xmax%s' % which_area)
+                northing_min = eval('project.ymin%s' % which_area)
+                northing_max = eval('project.ymax%s' % which_area)
+            except AttributeError:
+                print 'Unrecognized area name: %s' % which_area
+                break
 
-    for time_dir in time_dirs:
-        name1 = directory+time_dir+sep+project.scenario_name
-        name2 = directory+time_dir+sep+project.scenario_name+'_time_39900_0'
-        name3 = directory+time_dir+sep+project.scenario_name+'_time_79800_0'
-        names = [name1, name2, name3]
-       
-        for name in names:
-            for which_area in area:
-                if which_area == 'All':
-                    easting_min = None
-                    easting_max = None
-                    northing_min = None
-                    northing_max = None
-                else:
-                    try:
-                        easting_min = eval('project.xmin%s' % which_area)
-                        easting_max = eval('project.xmax%s' % which_area)
-                        northing_min = eval('project.ymin%s' % which_area)
-                        northing_max = eval('project.ymax%s' % which_area)
-                    except AttributeError:
-                        print 'Unrecognized area name: %s' % which_area
-                        break
-                  
+        for time_dir in time_dirs:
+
+            name1 = directory+time_dir+sep+project.scenario_name
+            name2 = directory+time_dir+sep+project.scenario_name+'_time_39900_0'
+            name3 = directory+time_dir+sep+project.scenario_name+'_time_79800_0'
+
+            names = [name1, name2, name3]
+      
+            asc_name = []   
+               
+            for name in names:
+                         
                 outname = name + '_' + which_area + '_' + which_var
                 quantityname = var_equations[which_var]
@@ -104 +119 @@
                 asc_name.append(outname + '.asc')
     
-        maxasc_outname = directory+time_dir+sep+project.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
+            maxasc_outname = directory+time_dir+sep+project.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
 
-        maxasc.MaxAsc(maxasc_outname, asc_name)
+            print 'max asc outname ', maxasc_outname
+            print 'asc_name ', str(asc_name)
+            
+            maxasc.MaxAsc(maxasc_outname, asc_name)