Changeset 7385 for anuga_work/production/dampier_2009/Arc_asc2raster_GDA94z50.py

Timestamp:

Aug 20, 2009, 9:20:31 AM (15 years ago)

Author:

kristy

Message:

 

File:

• anuga_work/production/dampier_2009/Arc_asc2raster_GDA94z50.py (modified) (3 diffs)

anuga_work/production/dampier_2009/Arc_asc2raster_GDA94z50.py

@@ -22 +22 @@
 gp.overwriteoutput = 1
 
-scenario_dir = "\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\onslow_tsunami_scenario_2009\\"
+scenario_dir = "\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\dampier_tsunami_scenario_2009\\"
 output_dir="anuga\\outputs\\"
 
-time_dir1 = '20090616_154908_run_final_0.0_27283_mhingee'
-##time_dir2 = '20090416_072538_run_final_0_27338_1307_Tb__kvanputt'
-##time_dir3 = '20090417_021634_run_final_0_70337_1307_Tb__kvanputt'
+time_dir1 = '20090817_150959_run_final_0.0_70408_kvanputt'
 
-##time_dir2 = '20090422_102457_run_final_0_70844_1307_Tb_internal_mhingee'
-
-time_dirs = [time_dir1]#, time_dir2]#, time_dir3]
+time_dirs = [time_dir1]
 
 for time_dir in time_dirs:
@@ -48 +44 @@
     
     #replication dictionary
-    replicate = (('onslow', ''),('_', ''),('max','_M'),
-                 ('CBD', 'CDB'),('All',''),
+    replicate = (('dampier', ''),('_', ''),('max','_M'),
+                 ('CBD', 'CDB'),('Dampier', 'Damp'),('Karratha', 'Karr')
                  ('depth','_depth'),('speed', '_speed'),
-                 ('elevation', '_ele_'), ('stage','_stage'))
+                 ('elevation', '_ele'), ('stage','_stage'))
 
     generate_filename = []
-    input_ascii = glob.glob(folder + '*elevation.asc')
+    input_ascii = glob.glob(folder + '*.asc')
 
     for infile in input_ascii:
@@ -72 +68 @@
 
         print 'Process: Define Projection'
-        #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]"