Changeset 6949 for anuga_work/production/australia_ph2/carnarvon

Timestamp:

May 5, 2009, 9:35:20 AM (16 years ago)

Author:

myall

Message:

re running carnarvon, and exporting results

Location:

anuga_work/production/australia_ph2/carnarvon

Files:

• Arc_asc2raster_GDA94z50.py (modified) (3 diffs)
• project.py (modified) (2 diffs)

anuga_work/production/australia_ph2/carnarvon/Arc_asc2raster_GDA94z50.py

@@ -25 +25 @@
 output_dir="anuga\\outputs\\"
 
-time_dir1 = '20090324_150124_run_final_0_27285_109_Bf_mhingee'
-time_dir2 = '20090324_161452_run_final_0_27285_109_Tb_mhingee'
+time_dir1 = '20090429_170742_run_trial_0_27285_84_Tb__kvanputt'
+time_dir2 = ''
 
-time_dirs = [time_dir2] #, time_dir2] #4, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6]
+time_dirs = [time_dir1] #, time_dir2] #4, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6]
  
 for time_dir in time_dirs:
@@ -52 +52 @@
 
     generate_filename = []
-    input_ascii = glob.glob(folder + '*.asc')
+    input_ascii = glob.glob(folder + '*elevation.asc')
 
     for infile in input_ascii:
@@ -70 +70 @@
 
         print 'Process: Define Projection' 
-        gp.DefineProjection_management(output_DEM, "PROJCS['GDA_1994_MGA_Zone_49',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',111.0],PARAMETER['Scale_Factor',0.9996]"
+        gp.DefineProjection_management(output_DEM, "PROJCS['CM_113.5',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',113.5],PARAMETER['Scale_Factor',0.9996]"
                                                    ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]")
 ##        output_extract = output_DEM + '_E'

anuga_work/production/australia_ph2/carnarvon/project.py

@@ -40 +40 @@
 friction=0.01           # manning's friction coefficient
 starttime=0             # start time for simulation
-finaltime=1000 # final time for simulation
-setup = 'trial'         # This can be one of three values
+finaltime=60000 # final time for simulation
+setup = 'final'         # This can be one of three values
                         #    trial - coarsest mesh, fast
                         #    basic - coarse mesh
@@ -128 +128 @@
 # Thinned ordering file from Hazard Map (geographic)
 # Format is index,latitude,longitude (with header)
-urs_order_filename = 'urs_order_simple.csv'
+urs_order_filename = 'urs_order.csv'
 
 # Landward bounding points