Changeset 6016

Timestamp:

Nov 28, 2008, 11:31:45 AM (16 years ago)

Author:

kristy

Message:

Changed from Maxasc in folder to in Lib

File:

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

anuga_work/production/busselton/export_results_max.py

@@ -13 +13 @@
 import project, os
 import sys
+from anuga.lib.maxasc.maxasc import MaxAsc
 from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts
 from anuga.shallow_water.data_manager import sww2dem
 from os import sep
-import maxasc
+
 
 directory = project.output_dir
 
-#time_dir1 = '20081002_111432_run_final_0_27283_250m_none_kvanputt'
-#time_dir1 = '20081002_111258_run_final_0_27283_250m_all_kvanputt'
-time_dir1 = '20081031_120642_run_final_0_27283_alpha0.1_kvanputt'
-#time_dir3 = '20080815_103818_run_final_0_polyline_newExtent_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]
+time_dir1 = '20081009_184721_run_final_0.6_68693_alpha0.1_kvanputt'
+time_dir2 = '20081031_120515_run_final_0_27255_alpha0.1_kvanputt'
+time_dir3 = '20081031_120753_run_final_0.6_27283_alpha0.1_kvanputt'
+time_dir4 = '20081031_120642_run_final_0_27283_alpha0.1_kvanputt'
+time_dir5 = '20081031_120316_run_final_0_68693_alpha0.1_kvanputt'
+time_dir6 = '20081111_113832_run_final_0.6_27255_alpha0.1_kvanputt'
 
-#cellsize = 20
-cellsize = 250
+time_dirs = [time_dir6] #1, time_dir2, time_dir3, time_dir4, time_dir5]
+  
+
+cellsize = 20
+#cellsize = 250
 
 timestep = None    # None means no timestep!
@@ -36 +46 @@
 # Set the special areas of interest.  If none, do: area='All'
 ######
-#area = ['Geordie', 'Sorrento', 'Fremantle', 'Rockingham']  # strings must match keys in var_equations below
-area = ['All']      # 'All' means no special areas - the whole thing
+
+area = ['Bunbury', 'Busselton']  # 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_39600_0' 
-        name3 = directory+time_dir+sep+project.scenario_name+'_time_79200_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_39600_0'
+            name3 = directory+time_dir+sep+project.scenario_name+'_time_79200_0'
+
+            names = [name1, name2, name3]
+      
+            asc_name = []   
+               
+            for name in names:
+                         
                 outname = name + '_' + which_area + '_' + which_var
                 quantityname = var_equations[which_var]
@@ -103 +118 @@
 
                 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)
+    
+            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_outname, asc_name)