Changeset 6971 for anuga_work/production/australia_ph2/melbourne_east

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/melbourne_east/Arc_asc2raster_GDA94z50.py (modified) (3 diffs)

anuga_work/production/australia_ph2/melbourne_east/Arc_asc2raster_GDA94z50.py

@@ -34 +34 @@
 time_dir7 = '20090421_152359_run_final_0_58115_2649_Tb_east_internal_mhingee'
 
-time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir7]
-##time_dirs = [time_dir4]
+##time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir7]
+time_dirs = [time_dir1, time_dir7]
     
 for time_dir in time_dirs:
@@ -44 +44 @@
 ##    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
@@ -58 +58 @@
 
     generate_filename = []
-    input_ascii = glob.glob(folder + '*max.asc')
+    input_ascii = glob.glob(folder + '*elevation.asc')
 
     for infile in input_ascii: