Changeset 5987 for anuga_work/production/geraldton/export_results_max.py

Timestamp:

Nov 20, 2008, 1:46:04 PM (16 years ago)

Author:

kristy

Message:

Ready to run export_results_max.py, copied from perth.
Set for event 27283 at 0.6 for one sww

File:

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

anuga_work/production/geraldton/export_results_max.py

@@ -20 +20 @@
 directory = project.output_dir
 
-time_dir1 = '20081102_104418_run_final_0_27283_alpha0.1_kvanputt'
-time_dir2 = '20081002_110950_run_final_0_27283_250m_all_kvanputt'
+time_dir1 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt'
+##time_dir2 = '20081114_155508_run_final_0.6_27283_alpha0.1_kvanputt' Perth
+##time_dir3 = '20081114_155610_run_final_0.6_27255_alpha0.1_kvanputt'
+##time_dir4 = '20081114_155642_run_final_0_27255_alpha0.1_kvanputt'
+##time_dir5 = '20081114_155725_run_final_0_68693_alpha0.1_kvanputt'
+##time_dir6 = '20081114_155802_run_final_0.6_68693_alpha0.1_kvanputt'
 
-time_dirs = [time_dir1, time_dir2]#, time_dir3]
+time_dirs = [time_dir1] #, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
 
-#cellsize = 20
-cellsize = 250
+  
+
+cellsize = 20
+#cellsize = 250
 
 timestep = None    # None means no timestep!
@@ -35 +41 @@
 ######
 
-#area = ['Geordie', 'Sorrento', 'Fremantle', 'Rockingham']  # strings must match keys in var_equations below
-area = ['All']      # 'All' means no special areas - the whole thing
+area = ['CBD']  # strings must match keys in var_equations below
+#area = ['All']      # 'All' means no special areas - the whole thing
 
 ######
@@ -49 +55 @@
 
 # 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:
@@ -59 +66 @@
         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_41700_0'
-        names = [name1, name2]        
-       
-        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]
@@ -99 +110 @@
                             verbose = True,
                             format = 'asc')
-
-                asc_name.append(outname + '.asc')
-    
-        maxasc_outname = directory+time_dir+sep+project.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
-
-        maxasc.MaxAsc(maxasc_outname, asc_name)
+##
+##                asc_name.append(outname + '.asc')
+##    
+##            maxasc_outname = directory+time_dir+sep+project.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
+##
+##            print 'max asc outname ', maxasc_outname
+##            print 'asc_name ', str(asc_name)
+##            
+##            maxasc.MaxAsc(maxasc_outname, asc_name)