Changeset 6942 for anuga_work/production/australia_ph2/albany

Timestamp:

May 1, 2009, 1:37:13 PM (16 years ago)

Author:

myall

Message:

 

File:

• anuga_work/production/australia_ph2/albany/Arc_asc2raster_GDA94z50.py (modified) (5 diffs)

anuga_work/production/australia_ph2/albany/Arc_asc2raster_GDA94z50.py

@@ -29 +29 @@
 time_dir3 = '20090409_010309_run_final_0_46697_1330_Tb__mhingee'
 
-events = [[time_dir1,0.395764],[time_dir2,0.390969],[time_dir3,0.392536]]
+time_dir4 = '20090421_165204_run_final_0_64344_1330_Tb_internal_mhingee'
 
-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_dir1, time_dir2]#, time_dir3]
+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 + '*stage.asc')
+    input_ascii = glob.glob(folder + '*elevation.asc')
 
     for infile in input_ascii:
@@ -75 +73 @@
 
         print 'Process: Define Projection'
-        # GDA_1994_MGA_Zone_54
+        #GDA_1994_MGA_Zone_54
         gp.DefineProjection_management(output_DEM, "PROJCS['GDA_1994_MGA_Zone_50',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"
                                                    ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',117.0],PARAMETER['Scale_Factor',0.9996]"
@@ -83 +81 @@
 ##        print 'Process: Extract by Mask'
 ##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
-        
-        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)
 
+