Changeset 6529 for anuga_work/production/australia_ph2/sydney/export_results_max.py

Timestamp:

Mar 17, 2009, 2:29:49 PM (14 years ago)

Author:

jgriffin

Message:

updated export_results_max.py

File:

• anuga_work/production/australia_ph2/sydney/export_results_max.py (modified) (6 diffs)

anuga_work/production/australia_ph2/sydney/export_results_max.py

@@ -1 +1 @@
 """
 Generates ascii grids of nominated areas -
-Input: sww file from run_pt_hedland.py
+Input: sww file from run_perth.py
        boundaries for grids from project.py
 Outputs: ascii grids of specified variables
@@ -13 +13 @@
 import project, os
 import sys
+from os.path import join
 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
-from os.path import join
+
 
 directory = project.output_folder
-print directory
 
-time_dir1 = '20090311_145740_run_final_0_58152_2956_Tb_jgriffin'
-time_dir2 = '20090311_145828_run_final_0_58152_2956_Bf_jgriffin'
-#time_dir3 = '20090123_154945_run_final_0_7875_250m_small_jgriffin'
-time_dirs = [time_dir1,time_dir2]#,time_dir2,time_dir3]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
+time_dir1 = '20090316_142137_run_final_0_58326_2956_Bf_jgriffin'
+##time_dir2 = '20090311_152929_run_final_0_70385_1203_Tb_jgriffin'
+##time_dir3 = '20090311_153009_run_final_0_70454_1203_Bf_jgriffin'
+##time_dir4 = '20090311_153054_run_final_0_70339_1203_Bf_jgriffin'
+time_dirs = [time_dir1]#,time_dir2,time_dir3,time_dir4]  
 
-#cellsize = 20
 cellsize = 250
+##cellsize = 5
 
-timestep = None    # None means no timestep!
-#timestep = 0
+#timestep = None    # None means no timestep!
+timestep = 0
 
 ######
@@ -37 +38 @@
 ######
 
-#area = ['']         # 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
 
 ######
@@ -51 +52 @@
 
 # one or more key strings from var_equations above
-#var = ['stage', 'speed','depth']
-var = ['stage','elevation']
+var = ['elevation','stage']
 
 ######
@@ -82 +82 @@
 
             name1 = join(directory, time_dir, project.scenario_name)
-            name2 = join(directory, time_dir, project.scenario_name)+'_time_66360_0'
-            #name3 = join(directory, time_dir, project.scenario_name)+'_time_33360_0'
-            #name4 = join(directory, time_dir, project.scenario_name)+'_time_50040_0'
-            #name5 = join(directory, time_dir, project.scenario_name)+'_time_66720_0'
-
-            names = [name1, name2]# , name3, name4, name5]
+##            name2 = join(directory, time_dir, project.scenario_name)+'_time_66360_0'
+##            name3 = join(directory, time_dir, project.scenario_name)+'_time_26640_0'
+##            name4 = join(directory, time_dir, project.scenario_name)+'_time_39960_0'
+##            name5 = join(directory, time_dir, project.scenario_name)+'_time_53280_0'
+##            name6 = join(directory, time_dir, project.scenario_name)+'_time_66600_0'
+##            name7 = join(directory, time_dir, project.scenario_name)+'_time_79920_0'
+            
+            names = [name1]#, name2]#, name3,name4,name5,name6,name7]
       
             asc_name = []   
@@ -112 +114 @@
                 asc_name.append(outname + '.asc')
     
-            maxasc_outname = directory+time_dir+sep+project.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
+            maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'
 
             print 'max asc outname ', maxasc_outname