Changeset 7351

Timestamp:

Aug 11, 2009, 9:03:37 AM (15 years ago)

Author:

jgriffin

Message:

 

Location:

anuga_work/production/new_south_wales

Files:

• batemans_bay/export_results_max.py (modified) (4 diffs)
• gosford/export_results_max.py (modified) (4 diffs)

anuga_work/production/new_south_wales/batemans_bay/export_results_max.py

@@ -24 +24 @@
 #Specify output directories
 time_dir1 = '20090529_143527_run_final_0.0_51424_jgriffin'
-#time_dir2 = '20090529_143458_run_final_0.0_58346_jgriffin'
+time_dir2 = '20090529_143458_run_final_0.0_58346_jgriffin'
 
-time_dirs = [time_dir1]#, time_dir2]   
+time_dirs = [time_dir1, time_dir2]   
 
+# sww filename extensions ie. batemans_bay_time_37860_0.sww, input into list 17640
+# make sure numbers are in sequential order
+times = [37860]
+
+#Modify the cellsize value to set the size of the raster you require
+#Take into account mesh size when aplying this paramater
 cellsize = 20   #dependent on data resolution in area of interest.
 
+#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    # None means no timestep!
 ##timestep = 0
@@ -54 +63 @@
 var = ['elevation']
 ######
-# 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 +89 @@
 
         for time_dir in time_dirs:
-
-            name1 = join(directory, time_dir, project.scenario_name)
-            name2 = join(directory, time_dir, project.scenario_name)+'_time_37860_0'
-
-            
-
-            names = [name1, name2]
+            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 = []   
@@ -110 +116 @@
                 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+'_max.asc'
 
-            print 'max asc outname ', maxasc_outname
-            print 'asc_name ', str(asc_name)
+                print 'max asc outname ', maxasc_outname
+                print 'asc_name ', str(asc_name)
             
-            MaxAsc(maxasc_outname, asc_name)
+                MaxAsc(maxasc_outname, asc_name)
+            else:
+                print 'only one sww input'

anuga_work/production/new_south_wales/gosford/export_results_max.py

@@ -23 +23 @@
 
 #Specify output directories
-time_dir1 = '20090615_161420_run_final_0.0_58025_jgriffin'
-#time_dir2 = '20090615_161318_run_final_0.0_58187_jgriffin'
+time_dir1 = '20090529_143527_run_final_0.0_51424_jgriffin'
+time_dir2 = '20090529_143458_run_final_0.0_58346_jgriffin'
 
-time_dirs = [time_dir1]#, time_dir2]
- 
+time_dirs = [time_dir1, time_dir2]   
 
+# sww filename extensions ie. batemans_bay_time_37860_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 = 20   #dependent on data resolution in area of interest.
 
+#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    # None means no timestep!
 ##timestep = 0
@@ -52 +60 @@
 
 # one or more key strings from var_equations above
+#var = ['depth', 'speed','stage']
 var = ['elevation']
-#var = ['depth', 'speed', 'stage']
 ######
-# 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
 ######
 
@@ -80 +89 @@
 
         for time_dir in time_dirs:
-
-            name1 = join(directory, time_dir, project.scenario_name)
-           
-            names = [name1]
+            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 = []   
@@ -108 +116 @@
                 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+'_max.asc'
 
-            print 'max asc outname ', maxasc_outname
-            print 'asc_name ', str(asc_name)
+                print 'max asc outname ', maxasc_outname
+                print 'asc_name ', str(asc_name)
             
-            MaxAsc(maxasc_outname, asc_name)
+                MaxAsc(maxasc_outname, asc_name)
+            else:
+                print 'only one sww input'