Changeset 7409

Timestamp:

Aug 27, 2009, 1:25:29 PM (16 years ago)

Author:

kristy

Message:

updated to new export script

File:

• anuga_work/production/patong/new_version/export_results_max.py (modified) (4 diffs)

anuga_work/production/patong/new_version/export_results_max.py

@@ -32 +32 @@
 
 
-time_dirs = [time_dir7] #, time_dir2,time_dir3, time_dir4,time_dir5, time_dir6, time_dir7]
+time_dirs = [time_dir4] #, time_dir2,time_dir3, time_dir4,time_dir5, time_dir6, time_dir7]
   
 
-#cellsize = 250
-cellsize = 20
+# sww filename extensions ie. hobart_time_17640_0.sww, input into list 17640
+# make sure numbers are in sequential order
+times = []
+            
+#Modify the cellsize value to set the size of the raster you require
+#Take into account mesh size when aplying this paramater
+cellsize = 10 #250
 
-timestep = None    # None means no timestep!
-#timestep = 0
+#Now set the timestep at which you want the raster generated.
+#Either set the actual timestep required or use 'None' to indicate that
+#you want the maximum values in the raster over all timesteps
+timestep = None   #0 
 
-######
-# Set the special areas of interest.  If none, do: area='All'
-######
+# Set the special areas of interest.  If none, do: area=['All']
+# Areas identified below are specified by cooridinates in project.py  
 
-area = ['CBD']  # strings must match keys in var_equations below
-#area = ['All']      # 'All' means no special areas - the whole thing
+area = ['CBD']  
+#area = ['All']      
+
+# one or more key strings from var_equations below
+var = ['speed']
 
 ######
@@ -58 +67 @@
                  'elevation': 'elevation' }
 
-# one or more key strings from var_equations above
-var = ['depth', 'speed']
 
 ######
-# Start running the various conversions we require.
+# Start script, running through variables, area, folder, sww file (determine by times)
+# Then creates a maximum asci if there is more than one sww file
 ######
 
@@ -87 +95 @@
 
         for time_dir in time_dirs:
-
-            name1 = join(directory, time_dir, project.scenario_name)
-            name2 = join(directory, time_dir, project.scenario_name)+'_time_10270_0'
-            name3 = join(directory, time_dir, project.scenario_name)+'_time_14555_0'
-
-            names = [name1, name2, name3]
-      
+            names = [join(directory, time_dir, project.scenario_name)]
+            for time in times:
+                names.append(join(directory, time_dir, project.scenario_name)+ '_time_' + str(time) + '_0')
+            
             asc_name = []   
-               
+                           
             for name in names:
-                         
-                outname = name + '_' + which_area + '_' + which_var
+                outname = name + '_' + which_area + '_' + which_var + '_' + str(cellsize)
                 quantityname = var_equations[which_var]
-
+                
                 print 'start sww2dem: time_dir=%s' % time_dir
-                
+                print 'name=%s' % name
                 sww2dem(name, basename_out = outname,
                             quantity = quantityname,
@@ -116 +120 @@
 
                 asc_name.append(outname + '.asc')
-    
-            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
+            
+            if len(names) >1:
+                maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'
+                +which_area+'_'+which_var+'_'+str(cellsize)+'_max.asc'
 
-            print 'max asc outname ', maxasc_outname
-            print 'asc_name ', str(asc_name)
-            
-            MaxAsc(maxasc_outname, asc_name)
+                print 'max asc outname ', maxasc_outname
+                print 'asc_name ', str(asc_name)
+                
+                MaxAsc(maxasc_outname, asc_name)
+            else:
+                print 'Only one sww file, no Max file required'