Changeset 6985 for anuga_work/production/australia_ph2/carnarvon

Timestamp:

May 6, 2009, 4:06:54 PM (16 years ago)

Author:

myall

Message:

making rasters

Location:

anuga_work/production/australia_ph2/carnarvon

Files:

• export_results.py (added)
• export_results_max.py (modified) (5 diffs)

anuga_work/production/australia_ph2/carnarvon/export_results_max.py

@@ -22 +22 @@
 directory = project.output_folder
 
-time_dir1 = '20090429_170742_run_trial_0_27285_84_Tb__kvanputt'
-#time_dir2 = '20090324_161452_run_final_0_27285_109_Tb_mhingee'
-##time_dir3 = ''
-##time_dir4 = ''
+time_dir1 = '20090504_114635_run_final_0_68754_84_Tb__mhingee'
+time_dir2 = '20090504_114715_run_final_0_70881_84_Tb__mhingee'
+time_dir3 = '20090505_100250_run_final_0_27285_84_Tb__mhingee'
+time_dir4 = '20090505_173023_run_final_0_27285_84_Tb_internal_mhingee'
 ##time_dir5 = ''
 ##time_dir6 = ''
 ##
-time_dirs = [time_dir1] #, time_dir2] # , time_dir4, time_dir5, time_dir6]
+time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]
   
 
@@ -35 +35 @@
 ##cellsize = 5
 
-##timestep = None    # None means no timestep!
-timestep = 0
+timestep = None    # None means no timestep!
+##timestep = 0
 
 ######
@@ -56 +56 @@
 
 # one or more key strings from var_equations above
-var = ['elevation']
+var = ['stage']
 
 ######
@@ -86 +86 @@
 
             name1 = join(directory, time_dir, project.scenario_name)
-##            name2 = join(directory, time_dir, project.scenario_name)+'_time_8760_0'
-##            name3 = join(directory, time_dir, project.scenario_name)+'_time_17520_0'
-##            name4 = join(directory, time_dir, project.scenario_name)+'_time_26280_0'
-##            name5 = join(directory, time_dir, project.scenario_name)+'_time_35040_0'
-##            name6 = join(directory, time_dir, project.scenario_name)+'_time_43800_0'
-##            name7 = join(directory, time_dir, project.scenario_name)+'_time_52560_0'
-##            name8 = join(directory, time_dir, project.scenario_name)+'_time_61320_0'
-##            name9 = join(directory, time_dir, project.scenario_name)+'_time_70080_0'
-##            name10 = join(directory, time_dir, project.scenario_name)+'_time_78840_0'
-##            
-            names = [name1] #, name2, name3, name4, name5, name6, name7, name8, name9, name10]
-      
-##            asc_name = []   
+            name2 = join(directory, time_dir, project.scenario_name)+'_time_31980_0'
+
+            names = [name1, name2]
+            
+            asc_name = []   
                
             for name in names:
@@ -119 +111 @@
                             format = 'asc')
 
-##                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)
+                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)