Changeset 7011 for anuga_work/production

Timestamp:

May 11, 2009, 3:43:32 PM (16 years ago)

Author:

kristy

Message:

 

Location:

anuga_work/production/australia_ph2

Files:

• carnarvon/Arc_asc2raster_GDA94z50.py (modified) (3 diffs)
• gold_coast/export_results_max.py (modified) (2 diffs)
• townsville/Arc_asc2raster_GDA94z50.py (modified) (1 diff)
• wyndham/Arc_asc2raster_GDA94z50.py (modified) (2 diffs)

anuga_work/production/australia_ph2/carnarvon/Arc_asc2raster_GDA94z50.py

@@ -25 +25 @@
 output_dir="anuga\\outputs\\"
 
-time_dir1 = '20090429_170742_run_trial_0_27285_84_Tb__kvanputt'
-time_dir2 = ''
+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_dirs = [time_dir1] #, time_dir2] #4, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6]
+time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]
  
 for time_dir in time_dirs:
@@ -35 +37 @@
     folder = scenario_dir + output_dir + time_dir +'\\'
     raster_gbd = folder + 'raster.gdb'
-    land = scenario_dir + "map_work\\Busselton.gdb\\Internal_polygons\\initial_conditions_extend_Di"
-    ocean = scenario_dir + "map_work\\Busselton.gdb\\input_boundaries\\Ocean"
-
-    
-    print 'Process: Create File GDB'
-    gp.CreateFileGDB_management(folder, "raster")
+        
+##    print 'Process: Create File GDB'
+##    gp.CreateFileGDB_management(folder, "raster")
 
     gp.Workspace = raster_gbd
@@ -73 +72 @@
                                                    ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',113.5],PARAMETER['Scale_Factor',0.9996]"
                                                    ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")
-##        output_extract = output_DEM + '_E'
-##        print 'Output Extract ',output_extract
-##        print 'Process: Extract by Mask'
-##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
 
-

anuga_work/production/australia_ph2/gold_coast/export_results_max.py

@@ -22 +22 @@
 directory = project.output_folder
 
-##time_dir1 = '20090311_152848_run_final_0_70385_1203_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 = '20090317_144100_run_final_0_51314_3727_Tb_jgriffin'
-time_dir1 = '20090421_150448_run_final_0_50883_3727_Tb_internal_mhingee'
-time_dirs = [time_dir1]#,time_dir2,time_dir3,time_dir4]  
+time_dir1 = '20090505_145250_run_final_0_31854_3727_Tb__lfountai'
+time_dir2 = '20090506_045043_run_final_0_64479_3727_Tb__lfountai'
+time_dir3 = '20090506_190156_run_final_0_50883_3727_Tb__lfountai'
+time_dirs = [time_dir1,time_dir2,time_dir3] #,time_dir4]  
 
 cellsize = 250
@@ -83 +81 @@
 
             name1 = join(directory, time_dir, project.scenario_name)
-            name2 = join(directory, time_dir, project.scenario_name)+'_time_36420_0'
-##            name2 = join(directory, time_dir, project.scenario_name)+'_time_21060_0'
-##            name3 = join(directory, time_dir, project.scenario_name)+'_time_42120_0'
-##            name4 = join(directory, time_dir, project.scenario_name)+'_time_63180_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]
+            name2 = join(directory, time_dir, project.scenario_name)+'_time_36480_0'
+
+            names = [name1, name2]
       
             asc_name = []   

anuga_work/production/australia_ph2/townsville/Arc_asc2raster_GDA94z50.py

@@ -27 +27 @@
 time_dir1 = '20090414_104809_run_final_0_31965_2740_Tb__kvanputt'
 time_dir2 = '20090420_064036_run_final_0_51469_2740_Tb__kvanputt'
-##time_dir3 = '20090429_145415_run_final_0_62890_2740_Tb__mhingee'
+time_dir3 = '20090421_141635_run_final_0_31965_2740_Tb_internal_lfountai'
+time_dir4 = '20090429_145415_run_final_0_62890_2740_Tb__mhingee'
 
-time_dir3 = '20090421_141635_run_final_0_31965_2740_Tb_internal_lfountai'
-
-##time_dirs = [time_dir1, time_dir2, time_dir3]
-time_dirs = [time_dir1, time_dir3]
+time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]
 
 for time_dir in time_dirs:

anuga_work/production/australia_ph2/wyndham/Arc_asc2raster_GDA94z50.py

@@ -25 +25 @@
 output_dir="anuga\\outputs\\"
 
-time_dir1 = '20090415_094009_run_final_0_17918_1823_Tb__kvanputt'
-time_dir2 = '20090418_182759_run_final_0_64962_1823_Tb__kvanputt'
-time_dir3 = '20090421_145541_run_final_0_70920_1823_Tb__kvanputt'
+##time_dir1 = '20090415_094009_run_final_0_17918_1823_Tb__kvanputt'
+##time_dir2 = '20090418_182759_run_final_0_64962_1823_Tb__kvanputt'
+##time_dir3 = '20090421_145541_run_final_0_70920_1823_Tb__kvanputt'
 
-##time_dir1 = '20090428_145948_run_final_0_64962_1823_Tb_internal_mhingee'
+time_dir1 = '20090428_145948_run_final_0_64962_1823_Tb_internal_mhingee'
 
-time_dirs = [time_dir2]#, time_dir2, time_dir3]
+time_dirs = [time_dir1]#, time_dir2, time_dir3]
 
 for time_dir in time_dirs:
@@ -54 +54 @@
 
     generate_filename = []
-    input_ascii = glob.glob(folder + '*elevation.asc')
+    input_ascii = glob.glob(folder + '*max.asc')
 
     for infile in input_ascii: