Changeset 6561

Timestamp:

Mar 20, 2009, 4:03:56 PM (14 years ago)

Author:

kristy

Message:

 

File:

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

anuga_work/production/busselton/Arc_asc2raster_GDA94z50.py

@@ -25 +25 @@
 output_dir="anuga\\outputs\\"
 
-##time_dir1 = '20081209_155610_run_final_0_27255_250m_none_lfountai'
-##time_dir2 = '20081209_160607_run_final_0_68693_250m_none_lfountai'
-
-##time_dir1 = '20081002_111432_run_final_0_27283_250m_none_kvanputt'
-##time_dir2 = '20081217_115336_run_final_0_27283_250m_none_dp_kvanputt'
-
-
-time_dir1 = '20081211_154006_run_final_0.6_27255_alpha0.1_kvanputt'
+time_dir1 = '20090130_165532_run_final_0.6_27255_extend_dt_kvanputt'
 ####time_dir2 = '20081211_162311_run_final_0_27255_alpha0.1_kvanputt'
 ####time_dir3 = '20081211_162346_run_final_0_68693_alpha0.1_kvanputt'
@@ -46 +39 @@
     folder = scenario_dir + output_dir + time_dir +'\\'
     raster_gbd = folder + 'raster.gdb'
-    land = scenario_dir + "map_work\\Busselton.gdb\\Internal_polygons\\initial_conditions_extend_Di"
+    land = scenario_dir + "map_work\\Busselton.gdb\\Internal_polygons\\initial_conditions_extend"
     ocean = scenario_dir + "map_work\\Busselton.gdb\\input_boundaries\\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 +51 @@
     
     #replication dictionary
-    replicate = (('busselton', ''),('_', ''),('max','_M_5'),
+    replicate = (('busselton', ''),('_', ''),('max','_M'),
                  ('Busselton', 'Bus'),('Bunbury', 'Bun'),
                  ('depth','_depth'),('speed', '_speed'),
@@ -64 +57 @@
 
     generate_filename = []
-    input_ascii = glob.glob(folder + '*Bunbury_depth_max.asc')
+    input_ascii = glob.glob(folder + '*.asc')
 
     for infile in input_ascii:
@@ -85 +78 @@
                                                    ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',117.0],PARAMETER['Scale_Factor',0.9996]"
                                                    ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")
-        output_extract = output_DEM + '_E'
-        print 'Output Extract ',output_extract
-        print 'Process: Extract by Mask'
-        gp.ExtractByMask_sa(output_DEM, land, output_extract)
+##        output_extract = output_DEM + '_E'
+##        print 'Output Extract ',output_extract
+##        print 'Process: Extract by Mask'
+##        gp.ExtractByMask_sa(output_DEM, land, output_extract)