Changeset 7355 for anuga_work/production/hobart_2009/export_results_max.py

Timestamp:

Aug 12, 2009, 2:56:52 PM (15 years ago)

Author:

kristy

Message:

 

File:

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

anuga_work/production/hobart_2009/export_results_max.py

@@ -21 +21 @@
 
 directory = project.output_folder
-time_dir1 = '20090619_165107_run_final_0_58260_None_kvanputt'
-time_dir2 = '20090619_165136_run_final_0.8_58260_None_kvanputt'
+time_dir1 = '20090714_162136_run_final_0_58260_None_kvanputt'
+time_dir2 = '20090714_162041_run_final_0.8_58292_None_kvanputt'
+time_dir3 = '20090714_161944_run_final_0.8_58280_None_kvanputt'
+time_dir4 = '20090714_161740_run_final_0.8_58260_None_kvanputt'
 
-time_dirs = [time_dir1, time_dir2] #, time_dir2] #, time_dir3, time_dir4, time_dir5, time_dir6]
-  
+time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]  
 
-#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 = [16980, 33960, 50940]
+            
+#Modify the cellsize value to set the size of the raster you require
+#Take into account mesh size when aplying this paramater
+cellsize = 20 #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']
-######
+# Areas identified below are specified by cooridinates in project.py  
 
-area = ['Hobart', 'NW', 'South']  
-#area = ['All']      # 'All' means no special areas - the whole thing
+area = ['South', 'NW', 'Hobart'] 
+#area = ['All']      
+
+# one or more key strings from var_equations below
+var = ['depth']
 
 ######
@@ -50 +60 @@
                  '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
 ######
 
@@ -79 +88 @@
 
         for time_dir in time_dirs:
-            name1 = join(directory, time_dir, project.scenario_name)
-            name2 = join(directory, time_dir, project.scenario_name)+'_time_17640_0'
-##            name3 = join(directory, time_dir, project.scenario_name)+'_time_33960_0'
-##            name4 = join(directory, time_dir, project.scenario_name)+'_time_50940_0'
+            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')
             
-            names = [name1, name2] #, name3, name4]
-      
             asc_name = []   
-               
+            print 'start sww2dem: time_dir = %s' % time_dir               
             for name in names:
-                         
+
                 outname = name + '_' + which_area + '_' + which_var
                 quantityname = var_equations[which_var]
 
-                print 'start sww2dem: time_dir=%s' % time_dir
+                print 'start sww2dem: name = %s' % os.path.basename(name)
                 
                 sww2dem(name, basename_out = outname,
@@ -108 +114 @@
 
                 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)
+            #print 'list = %s' % asc_name
+            #print 'max asc outname ', maxasc_outname
+            #print 'name=%s' % name
+            print 'which area=%s' % which_area
+            print 'which var=%s' % which_var
+            #print 'asc_name ', str(asc_name)
+            if times == []:     
+                print 'Only one sww file, no Max file required'
+            else:
+                print 'Multiple sww files'
+                maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
+                print 'maxasc_outname = ', maxasc_outname
+                MaxAsc(maxasc_outname, asc_name)