Changeset 6971 for anuga_work/production/australia_ph2/geraldton

Timestamp:

May 5, 2009, 5:45:54 PM (16 years ago)

Author:

myall

Message:

converting ascii to raster and making images

File:

• anuga_work/production/australia_ph2/geraldton/Arc_asc2raster_GDA94z50.py (modified) (3 diffs)

anuga_work/production/australia_ph2/geraldton/Arc_asc2raster_GDA94z50.py

@@ -26 +26 @@
 
 time_dir1 = '20090414_114540_run_final_0_64475_576_Tb__kvanputt'
-time_dir2 = '20090415_014342_run_final_0_70921_576_Tb__kvanputt'
-time_dir3 = '20090415_153319_run_final_0_26985_576_Tb__kvanputt'
+##time_dir2 = '20090415_014342_run_final_0_70921_576_Tb__kvanputt'
+##time_dir3 = '20090415_153319_run_final_0_26985_576_Tb__kvanputt'
 
-time_dirs = [time_dir1, time_dir2, time_dir3]
+time_dir2 = '20090421_172030_run_final_0_64475_576_Tb_internal_mhingee'
+
+time_dirs = [time_dir1, time_dir2]#, time_dir3]
 
 for time_dir in time_dirs:
@@ -38 +40 @@
 ##    ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_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
@@ -52 +54 @@
 
     generate_filename = []
-    input_ascii = glob.glob(folder + '*max.asc')
+    input_ascii = glob.glob(folder + '*elevation.asc')
 
     for infile in input_ascii: