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

Timestamp:

Mar 12, 2009, 10:57:08 AM (14 years ago)

Author:

jgriffin

Message:

changed event number; added Arc_asc2raster_GDA94z56.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

@@ -2 +2 @@
 Generates ascii grids of nominated areas -
 Input: sww file from run_pt_hedland.py
-       boundaries for grids from project_250m.py
+       boundaries for grids from project.py
 Outputs: ascii grids of specified variables
 Stored in the 'outputs_dir' folder for respective .sww file
@@ -11 +11 @@
 """
 
-import project_250m, os
+import project, os
 import sys
 from anuga.lib.maxasc.maxasc import MaxAsc
@@ -17 +17 @@
 from anuga.shallow_water.data_manager import sww2dem
 from os import sep
+from os.path import join
 
-directory = project_250m.output_dir
+directory = project.output_folder
+print directory
 
-time_dir1 = '20090206_150250_run_final_0_250m_large_jgriffin'
-#time_dir2 = '20090123_155045_run_final_0_7875_250m_large_jgriffin'
+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_dir3]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
+time_dirs = [time_dir1,time_dir2]#,time_dir2,time_dir3]#, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
 
 #cellsize = 20
@@ -69 +71 @@
         else:
             try:
-                easting_min = eval('project_250m.xmin%s' % which_area)
-                easting_max = eval('project_250m.xmax%s' % which_area)
-                northing_min = eval('project_250m.ymin%s' % which_area)
-                northing_max = eval('project_250m.ymax%s' % which_area)
+                easting_min = eval('project.xmin%s' % which_area)
+                easting_max = eval('project.xmax%s' % which_area)
+                northing_min = eval('project.ymin%s' % which_area)
+                northing_max = eval('project.ymax%s' % which_area)
             except AttributeError:
                 print 'Unrecognized area name: %s' % which_area
@@ -79 +81 @@
         for time_dir in time_dirs:
 
-            name1 = directory+time_dir+sep+project_250m.scenario_name
-            #name2 = directory+time_dir+sep+project_250m.scenario_name+'_time_52620_0'
-           # name3 = directory+time_dir+sep+project_250m.scenario_name+'_time_78360_0'
+            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]
+            names = [name1, name2]# , name3, name4, name5]
       
             asc_name = []   
@@ -108 +112 @@
                 asc_name.append(outname + '.asc')
     
-            maxasc_outname = directory+time_dir+sep+project_250m.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
+            maxasc_outname = directory+time_dir+sep+project.scenario_name+'_'+which_area+'_'+which_var+'_max.asc'
 
             print 'max asc outname ', maxasc_outname