Changeset 6985 for anuga_work/production/australia_ph2/eucla_motel/Arc_asc2raster_GDA94z50.py

Timestamp:

May 6, 2009, 4:06:54 PM (15 years ago)

Author:

myall

Message:

making rasters

File:

• anuga_work/production/australia_ph2/eucla_motel/Arc_asc2raster_GDA94z50.py (modified) (4 diffs)

anuga_work/production/australia_ph2/eucla_motel/Arc_asc2raster_GDA94z50.py

@@ -28 +28 @@
 time_dir2 = '20090409_073114_run_final_0_64449_1837_Tb__mhingee'
 time_dir3 = '20090409_222301_run_final_0_68755_1837_Tb__mhingee'
+time_dir4 = '20090421_164640_run_final_0_64449_1837_Tb_internal_mhingee'
 
-events = [[time_dir1,0.265183],[time_dir2,0.27574],[time_dir3,0.267985]]
+##time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]
 
-for event in events:
-##for time_dir in time_dirs:
-    time_dir = event[0]
-    max_wave = event[1]
-    print time_dir
-    print max_wave
+time_dirs = [time_dir2]
+
+for time_dir in time_dirs:
     # Local variables...
     folder = scenario_dir + output_dir +  time_dir +'\\'
@@ -43 +41 @@
 ##    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
@@ -57 +55 @@
 
     generate_filename = []
-    input_ascii = glob.glob(folder + '*max.asc')
+    input_ascii = glob.glob(folder + '*elevation.asc')
 
     for infile in input_ascii:
@@ -83 +81 @@
 ##        print 'Process: Extract by Mask'
 ##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
-        # do this bit only if there are only stage asc files
-        div_file = output_DEM.replace('stage','div_stage')
-        if not div_file == output_DEM:
-            print 'divide', output_DEM,' by', max_wave,' and call',div_file
-            gp.Divide_sa(output_DEM,max_wave,div_file)
-
+