Changeset 6842 for anuga_work/production/australia_ph2/sydney

Timestamp:

Apr 20, 2009, 4:45:56 PM (16 years ago)

Author:

myall

Message:

using ARC script to turn ascii stage files into rasters, and dividing by max wave height; done for GC clockwise to perth. not checked thoroughly.

File:

• anuga_work/production/australia_ph2/sydney/Arc_asc2raster_GDA94z56.py (modified) (4 diffs)

anuga_work/production/australia_ph2/sydney/Arc_asc2raster_GDA94z56.py

@@ -25 +25 @@
 output_dir = "anuga\\outputs\\"
 
-time_dir1 = '20090331_150428_run_final_0_58152_2956_Tb_lfountai'
-##time_dir2 = '20081210_100528_run_final_0_68693_250m_none_kvanputt'
-##time_dirs = [time_dir1, time_dir2]
+time_dir1 = '20090408_171539_run_final_0_58152_2938_Tb__kvanputt'
+time_dir2 = '20090409_011331_run_final_0_58222_2938_Tb__kvanputt'
+time_dir3 = '20090409_085101_run_final_0_58326_2938_Tb__kvanputt'
+time_dirs = [time_dir1, time_dir2, time_dir3]
 
 ##time_dir1 = '20090311_145740_run_final_0_58152_2956_Tb_jgriffin'
@@ -37 +38 @@
 
 
-time_dirs = [time_dir1]#, time_dir4 , time_dir3, time_dir4, time_dir5]#, time_dir6]
+##time_dirs = [time_dir1]#, time_dir4 , time_dir3, time_dir4, time_dir5]#, time_dir6]
 
 for time_dir in time_dirs:
@@ -59 +60 @@
                  ('_', ''),('Geordie', 'Geo'),('Sorrento', 'Sor'), ('max','M_'),
                  ('Fremantle', 'Fre'),('Rockingham', 'Roc'),('depth','_dep_'),
-                 ('speed', '_spe_'), ('elevation', '_ele_'), ('stage','_stage'))
+                 ('speed', '_spe_'), ('elevation', '_ele_'), ('stage','_st'))
 
     generate_filename = []
@@ -80 +81 @@
 
         print 'Process: Define Projection' 
-        gp.DefineProjection_management(output_DEM, "PROJCS['GDA_1994_MGA_Zone_56',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',153.0],PARAMETER['Scale_Factor',0.9996]"
+        gp.DefineProjection_management(output_DEM, "PROJCS['CM_151',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',151.0],PARAMETER['Scale_Factor',0.9996]"
                                                    ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")