Changeset 2615 for production/onslow_2006/project_new.py

Timestamp:

Mar 28, 2006, 11:00:20 AM (19 years ago)

Author:

nick

Message:

update onslow

File:

• production/onslow_2006/project_new.py (moved) (moved from production/onslow_2006/project.py) (4 diffs)

production/onslow_2006/project_new.py

@@ -9 +9 @@
 
 from pmesh.create_mesh import convert_points_from_latlon_to_utm
+
+from coordinate_transforms.redfearn import degminsec2decimal_degrees
+
+from time import localtime, strftime
+
+from os import getcwd
                 
 #Making assumptions about the location of scenario data
 scenario_dir_name = 'onslow_tsunami_scenario_2006'
 
-
 # 250m data to be provided
 coarsename = 'onsl_bathydem250' # get from Neil/Ingo (DEM or topo data)
-"""
+
 # 30m data to be provided
-finename = 'onshore_30' # get from Neil/Ingo (DEM or topo data) 
-"""
+onshore_name = 'onslow_onshore_30m_dted' # get from Neil/Ingo (DEM or topo data) 
 
-# clipping region for fine elevation data
-west = 240000
-east = 340000
-south = 7570000
-north = 7645000
+offshore_name = 'onslow_offshore_points'
+
+boundary_basename = 'SU-AU'
 
 #swollen/ all data output
 basename = 'source'
+
+codename = 'project.py'
 
 if sys.platform == 'win32':
@@ -36 +40 @@
 
 #Derive subdirectories and filenames
+timedir = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
+outputdir = home+sep+scenario_dir_name+sep+'output'+sep+timedir+sep
 meshdir = home+sep+scenario_dir_name+sep+'meshes'+sep
 datadir = home+sep+scenario_dir_name+sep+'topographies'+sep
-outputdir = home+sep+scenario_dir_name+sep+'output'+sep
+
 polygondir = home+sep+scenario_dir_name+sep+'polygons'+sep
 boundarydir = home+sep+scenario_dir_name+sep+'boundaries'+sep
+
+codedir = getcwd()+sep
+                                
+codedirname = codedir + 'project.py'
 
 meshname = meshdir + basename
 
 coarsedemname = datadir + coarsename
-"""
-finedemname = datadir + finename
-combineddemname = datadir + 'onslow_combined_elevation'
-"""
+
+onshore_dem_name = datadir + onshore_name
+
+offshore_dem_name = datadir + offshore_name
+
+combined_dem_name = datadir + 'onslow_combined_elevation'
 
 outputname = outputdir + basename  #Used by post processing
@@ -54 +66 @@
 #!gauge_filename = outputdir + 'onslow_gauges.xya'
 #!gauge_outname = outputdir + 'gauges_max_output.xya'
+
+# clipping region for fine elevation data
+eastingmin = 250000
+eastingmax = 330000
+northingmin = 7580000
+northingmax = 7635000
+
+south = degminsec2decimal_degrees(-22,00,0)
+north = degminsec2decimal_degrees(-21,10,0)
+west = degminsec2decimal_degrees(114,30,0)
+east = degminsec2decimal_degrees(115,30,0)
+
+# region for visualisation
+eminviz = 260000
+emaxviz = 320000
+nminviz = 7590000
+nmaxviz = 7630000
 
 #Georeferencing
@@ -97 +126 @@
 poly_direction = [k0, k1, k2, k3]
 
-#!slump_origin = [385000.0, 6255000.0] #Absolute UTM