Changeset 7354 for anuga_work/production/hobart_2009

Timestamp:

Aug 12, 2009, 2:54:13 PM (16 years ago)

Author:

kristy

Message:

Updated all scripts since new runs

File:

• anuga_work/production/hobart_2009/Arc_asc2raster_GDA94z55.py (modified) (6 diffs)

anuga_work/production/hobart_2009/Arc_asc2raster_GDA94z55.py

@@ -26 +26 @@
 output_dir = "anuga\\outputs\\"
 
-##time_dir1 = '20090601_133348_run_final_0.8_58280_None_kvanputt'
-##time_dir2 = '20090601_133503_run_final_0.8_58292_None_kvanputt'
-##time_dir3 = '20090505_150711_run_final_0_58280_None_kvanputt'
-##time_dir4 = '20090505_150517_run_final_0_58292_None_kvanputt'
-##time_dir5 = '20090505_151322_run_final_0.8_64477_None_kvanputt'
-##time_dir6 = '20090505_151447_run_final_0_64477_None_kvanputt'
-##
-##time_dirs = [time_dir4] #1, time_dir2] #, time_dir3, time_dir4, time_dir5, time_dir6]
+time_dir1 = '20090714_162136_run_final_0_58260_None_kvanputt'
+time_dir2 = '20090714_162041_run_final_0.8_58292_None_kvanputt'
+time_dir3 = '20090714_161944_run_final_0.8_58280_None_kvanputt'
+time_dir4 = '20090714_161740_run_final_0.8_58260_None_kvanputt'
 
-time_dir1 = '20090619_165107_run_final_0_58260_None_kvanputt'
-time_dir2 = '20090619_165136_run_final_0.8_58260_None_kvanputt'
+time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4]
 
-time_dirs = [time_dir1, time_dir2]
 for time_dir in time_dirs:
     # Local variables...
@@ -44 +38 @@
     raster_gbd = folder + 'raster.gdb'
     land_folder = scenario_dir + "ArcGIS\\data.gdb\\Land"
-    areas = ['NW', 'South', 'Hobart']  
+    areas = ['NW', 'South', 'Hobart'] 
     print time_dir
 
@@ -55 +49 @@
         #replication dictionary
         replicate = (('hobart', ''), ('fitting_problem', ''),
-                     ('_', ''), ('max','_M'), ('depth','_depth'),
+                     ('_', ''), ('Hobart','H'), ('max','_M'), ('depth','_depth'),
                      ('speed', '_speed'), ('elevation', '_ele'), ('stage','_stage'))
 
@@ -70 +64 @@
             sys.exit(10)
         generate_filename.append(output_DEM)
+        output_DEM = raster_gbd + "\\" + output_DEM 
         print 'Output DEM ',output_DEM
        
@@ -82 +77 @@
                                        ",PARAMETER['Central_Meridian',147.0],PARAMETER['Scale_Factor',0.9996]"
                                        ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")
-
+               
 ##        land = land_folder + '_' + area
 ##        print 'Process: Extract by Mask'
@@ -88 +83 @@
 ##        print 'Output Extract ',output_extract
 ##        gp.ExtractByMask_sa(output_DEM, land, output_extract)
-##
-##
+##        print 'Process: Delete'
+##        gp.Delete_management(output_DEM, "RasterDataset")
+