Changeset 2902 for production/pt_hedland_2006/project.py

Timestamp:

May 18, 2006, 10:01:58 AM (18 years ago)

Author:

nick

Message:

update onslow and pt hedland

File:

• production/pt_hedland_2006/project.py (modified) (8 diffs)

production/pt_hedland_2006/project.py

@@ -3 +3 @@
 """
 
-from os import sep, environ
+from os import sep, environ, getenv, getcwd
 from os.path import expanduser
 from utilities.polygon import read_polygon
@@ -10 +10 @@
 from coordinate_transforms.redfearn import degminsec2decimal_degrees
 from time import localtime, strftime
-from os import getcwd
 
 #Making assumptions about the location of scenario data
@@ -32 +31 @@
 
 if sys.platform == 'win32':
-    home = environ['INUNDATIONHOME']     #Sandpit's parent dir
+    home = getenv('INUNDATIONHOME')
+#    python_home = getenv('PWD')     
+#    home = environ['INUNDATIONHOME']     #Sandpit's parent dir
 else:    
-    home = expanduser('~')
+#    home = environ['INUNDATIONHOME']     #Sandpit's parent dir
+#    home = getenv('INUNDATIONHOME')     
+#    python_home = getenv('PWD')     
+    home = getenv('INUNDATIONHOME', sep+'d'+sep+'cit'+sep+'1'+sep+'cit'+sep+'risk_assessment_methods_project'+sep+'inundation')     
+#    home = sep+'d'+sep+'cit'+sep+'1'+sep+'cit'+sep+'risk_assessment_methods_project'+sep+'inundation'
+#    home = expanduser('~')
 
 #Derive subdirectories and filenames
 time = strftime('%Y%m%d_%H%M%S',localtime()) #gets time for new dir
+#print 'home', home
 outputtimedir = home+sep+scenario_dir_name+sep+'output'+sep+time+sep
+#print 'outputtimedir', outputtimedir
 meshdir = home+sep+scenario_dir_name+sep+'meshes'+sep
 datadir = home+sep+scenario_dir_name+sep+'topographies'+sep
@@ -48 +56 @@
 tidedir = home+sep+scenario_dir_name+sep+'tide_data'+sep
 
-print'bound', boundarydir
+#print'bound', boundarydir
 
 #gauge_filename = gaugedir + 'onslow_gauges.xya'
@@ -59 +67 @@
 #MOST_dir = 'f:'+sep+'3'+sep+'ehn'+sep+'users'+sep+'davidb'+sep+'tsunami'+sep+'WA_project'+sep+'SU-AU_90'+sep+'most_2'+sep+'detailed'+sep
 
-codedir = getcwd()+sep
-                                
-codedirname = codedir + 'project.py'
+#print 'name', __name__
+#print 'path', __file__
+#codedir = getcwd()+sep
+
+#project_code_name = __name__ 
+                           
+#project_code_dir_name = __file__
 
 meshname = meshdir + basename
@@ -81 +93 @@
 
 # for ferret2sww
-#south = 
-#north = 
-#west = 
-#east = 
+south = degminsec2decimal_degrees(-20,30,0)
+north = degminsec2decimal_degrees(-17,10,0)
+west = degminsec2decimal_degrees(117,00,0)
+east = degminsec2decimal_degrees(120,00,0)
 
 # region to export (used from export_results.py)
@@ -96 +108 @@
 
 # bounding polygon provided by Hamish
-d0 = [818732.55, 8062768.27]
+#d0 = [818732.55, 8062768.27]
+d0 = [755000.0, 8025000.0]
 d1 = [708940.32, 7750510.33]
 d2 = [656561.15, 7732615.11]
@@ -102 +115 @@
 d4 = [517682.34, 7899310.22]
 
-polyAll = [d0, d1, d2, d3, d4]
+polyAll = [d0, d4, d3, d2, d1]
 
 #Interior region - Pt Hedland town
-i0 = [690000, 7750510.33] # these are just currently close to d1, d2 and d3
-i1 = [656561.15, 7750000]
-i2 = [620000, 7733013.56]
+i0 = [668000, 7757000]
+i1 = [659000, 7755000]
+i2 = [660000, 7749000]
+i3 = [667000, 7746000]
+i4 = [678000, 7751000]
 
-poly_pt_hedland = [i0, i1, i2]
+poly_pt_hedland = [i0, i1, i2, i3, i4]
 
 #Are there other significant features?
-j0 = []
-j1 = []
-j2 = []
-j3 = []
+j0 = [670000, 7780000]
+j1 = [620000, 7745000]
+j2 = [665000, 7735000]
+j3 = [700000, 7755000]
 
-poly_ = [j0, j1, j2, j3]
+poly_region = [j0, j1, j2, j3]