Changeset 6041 for anuga_work/production/geraldton/Arc_asc2raster_GDA94z50.py

Timestamp:

Dec 1, 2008, 11:48:14 AM (16 years ago)

Author:

kristy

Message:

change of name

File:

• anuga_work/production/geraldton/Arc_asc2raster_GDA94z50.py (moved) (moved from anuga_work/production/geraldton/asc2raster_GDA94z50.py) (5 diffs)

anuga_work/production/geraldton/Arc_asc2raster_GDA94z50.py

@@ -26 +26 @@
 
 
-time_dir1 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt'
-##time_dir2 = ''
-##time_dir3 = ''
-##time_dir4 = ''
-##time_dir5 = ''
-##time_dir6 = ''
+time_dir1 = '20081117_141353_run_final_0_27283_alpha0.1_kvanputt'
+time_dir2 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt' 
+time_dir3 = '20081117_141520_run_final_0_27255_alpha0.1_kvanputt'
+time_dir4 = '20081117_141558_run_final_0.6_27255_alpha0.1_kvanputt'
+time_dir5 = '20081117_141634_run_final_0_68693_alpha0.1_kvanputt'
+time_dir6 = '20081117_141717_run_final_0.6_68693_alpha0.1_kvanputt'
 
-time_dirs = [time_dir1] #, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
- 
+time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
+
 for time_dir in time_dirs:
 
@@ -43 +43 @@
 
     # Process: Create File GDB...
-    gp.CreateFileGDB_management(folder, "raster")
+    #gp.CreateFileGDB_management(folder, "raster")
 
     gp.Workspace = raster_gbd
@@ -56 +56 @@
 
     generate_filename = []
-    input_ascii = glob.glob(folder + '*.asc')
+    input_ascii = glob.glob(folder + '*max.asc')
 
     for infile in input_ascii:
@@ -70 +70 @@
         print 'Output DEM ',output_DEM
         output_extract = output_DEM + 'E'
-##        print 'Output Extract ',output_extract
+        print 'Output Extract ',output_extract
         
         print 'Process: ASCII to Raster'
@@ -80 +80 @@
                                                    ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")
 
-##        print 'Process: Extract by Mask'
-##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
+        print 'Process: Extract by Mask'
+        gp.ExtractByMask_sa(output_DEM, land, output_extract)