Changeset 7273

Timestamp:

Jun 29, 2009, 4:27:55 PM (16 years ago)

Author:

kristy

Message:

 

Location:

anuga_work/production/exmouth_2009

Files:

• Arc_asc2raster_GDA94z50.py (modified) (4 diffs)
• build_elevation.py (modified) (1 diff)
• export_results.py (modified) (4 diffs)

anuga_work/production/exmouth_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\\exmouth_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 = '20090624_145617_run_final_0.0_27283_mhingee'
+##time_dir2 = ''
+##time_dir3 = ''
 
-##time_dir2 = '20090422_102457_run_final_0_70844_1307_Tb_internal_mhingee'
+##time_dir2 = ''
 
 time_dirs = [time_dir1]#, time_dir2]#, time_dir3]
@@ -48 +48 @@
     
     #replication dictionary
-    replicate = (('onslow', ''),('_', ''),('max','_M'),
+    replicate = (('exmouth', ''),('_', ''),('max','_M'),
                  ('CBD', 'CDB'),('All',''),
                  ('depth','_depth'),('speed', '_speed'),
@@ -54 +54 @@
 
     generate_filename = []
-    input_ascii = glob.glob(folder + '*elevation.asc')
+    input_ascii = glob.glob(folder + '*.asc')
 
     for infile in input_ascii:
@@ -73 +73 @@
         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]"
+        gp.DefineProjection_management(output_DEM, "PROJCS['CM_114',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',114.0],PARAMETER['Scale_Factor',0.9996]"
                                                    ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")
 ##        output_extract = output_DEM + 'E'

anuga_work/production/exmouth_2009/build_elevation.py

@@ -45 +45 @@
 #------------------------------------------------------------------------------
 
-print 'project.bounding_polygon', project.bounding_polygon
-print 'project.combined_elevation_basename', project.combined_elevation_basename
+##print 'project.bounding_polygon', project.bounding_polygon
+##print 'project.combined_elevation_basename', project.combined_elevation_basename
 
 # Create Geospatial data from ASCII files

anuga_work/production/exmouth_2009/export_results.py

@@ -22 +22 @@
 directory = project.output_folder
 
-time_dir1 = '20090616_154908_run_final_0.0_27283_mhingee'
-time_dir2 = '20090617_110255_run_final_0.0_27283_mhingee'
+time_dir1 = '20090624_145617_run_final_0.0_27283_mhingee'
+##time_dir2 = ''
 ##time_dir3 = ''
 ##time_dir4 = ''
@@ -29 +29 @@
 ##time_dir6 = ''
 ##
-time_dirs = [time_dir2]#, time_dir2]# , time_dir4, time_dir5, time_dir6]
+time_dirs = [time_dir1]#, time_dir2]# , time_dir4, time_dir5, time_dir6]
   
 
@@ -42 +42 @@
 ######
 
-#area = ['Bunbury', 'Busselton']  # strings must match keys in var_equations below
+#area = ['Exmouth']  # strings must match keys in var_equations below
 area = ['All']      # 'All' means no special areas - the whole thing
 
@@ -56 +56 @@
 
 # one or more key strings from var_equations above
-var = ['elevation']
+var = ['elevation','stage']
 
 ######