Changeset 2487 for production/sydney_2006/project.py

Timestamp:

Mar 7, 2006, 10:33:14 AM (19 years ago)

Author:

sexton

Message:

New run_sydney file incorporating (1) reading in polygons (created by GIS guys) to create the interior regions, (2) allocating a test data set (defined in project.py) for elevation to test for the black screen of death and (3) track the volume in time (plot_integral.py is the script to plot the output)

File:

• production/sydney_2006/project.py (modified) (12 diffs)

production/sydney_2006/project.py

@@ -10 +10 @@
 from pmesh.create_mesh import convert_points_from_latlon_to_utm
                 
-
-
 
 #Making assumptions about the location of scenario data
@@ -20 +18 @@
 finename = 'bathy_dem25' # get from Neil/Ingo (DEM or topo data) Wed 25 Jan
 
-
-
 # creating easting and northing max and min for fine data - region of interest
 eastingmin = 332090
@@ -34 +30 @@
 nmaxviz = 6283000
 
-basename = 'slump_TEST'
-basename2 = 'slump_newfriction'
+basename = 'slump_friction04'
+basename4 = 'slump_poly_ingo_test'
 
 if sys.platform == 'win32':
@@ -42 +38 @@
     home = expanduser('~')
     
-
-
 #Derive subdirectories and filenames
 meshdir = home+sep+scenario_dir_name+sep+'meshes'+sep
@@ -49 +43 @@
 outputdir = home+sep+scenario_dir_name+sep+'output'+sep
 polygondir = home+sep+scenario_dir_name+sep+'polygons'+sep
+gaugedir = home+sep+scenario_dir_name+sep+'gauges'+sep
 
 meshname = meshdir + basename
+meshname4 = meshdir + basename4
 coarsedemname = datadir + coarsename
 finedemname = datadir + finename
 combineddemname = datadir + 'sydneytopo'
 outputname = outputdir + basename  #Used by post processing
-outputname2 = outputdir + basename2  #Used by post processing
+outputname4 = outputdir + basename4  #Used by post processing
 
 #csv file of coastline 50m epsilon belt
@@ -62 +58 @@
 #print manly_polygon
 
-gauge_filename = outputdir + 'sydney_gauges.xya'
-gauge_outname = outputdir + 'gauges_max_output.xya'
-polygonptsfile = outputdir + 'poly'
+gauge_filename = gaugedir + 'sydney_gauges.xya'
+gauge_outname = gaugedir + 'gauges_max_output.xya'
+#gauge_filename = gaugedir + 'nest_gauges_Manly.xya'
+#gauge_filename = gaugedir + 'west_of_quay_yprofile.xya'
+#gauge_filename = gaugedir + 'GA_gauge.csv' # from Benfield
+#gauge_outname = gaugedir + 'gauges_max_output_next.xya'
+gaugetimeseries = gaugedir + 'gauges_time_series_Benfield'
+polygonptsfile = polygondir + 'poly'
+integraltimeseries = outputdir + 'integral_time_series'
 
 #Georeferencing
@@ -181 +183 @@
 np1 = [318200, 6253000]
 np2 = [327000, 6250500]
+np2east = [327200, 6250500]
+np2west = [326800, 6250500]
 np3 = [333000, 6249000]
 np4 = [342000, 6249000]
@@ -186 +190 @@
 np6 = [343000, 6256000]
 np7 = [344000, 6258000]
-np8 = [343000, 6260000]
+np8 = [343000, 6260300]
+#np8 = [343000, 6260000]
 np9 = [343000, 6264000]
 np10 = [345000, 6265000]
@@ -216 +221 @@
 np34 = [337000, 6256000]
 np35 = [338000, 6255000]
-np36 = [319000, 6254000]
+np36 = [320000, 6254000]
 np37 = [322000, 6252000]
 np38 = [324000, 6253000]
@@ -234 +239 @@
 #             np31, np32, \
 #             np23]
-newpoly1 = [np2, np3, \
-            np28, np29, np30, \
+newpoly1 = [np28, np29, np30, \
             np5, np6, np7, np8, np9, np10, \
             np11, np12, np13, np14, np15, \
             np27, np26, np25, np24, np17, \
-            np18, np19, np20, np21, np31]
+            np33, \
+            np19, np20, np21, np22, np31, np32, \
+            np36, np37, np38, np39, np2, np3]
+# last two lines for second run
+
+north1 = [np2, np3, np30, np5, np24, np17, np33, np19, np20, np21, \
+          np31, np32, np37, np38, np39]
+north2 = [np8, np9, np19, np11, np12, np13, np14, np15, np27, np27, np25]
+npinsert = [332000, 6255000]
+parrariver = [np3, np18, npinsert, np19, np32, np36, np38, np39, np2] #first run
+#parrariver = [np32, np36, np37, np38, np39, np2west, np31] #second run
 
 # south
@@ -266 +280 @@
 sp23 = [339000, 6236000]
 
-newpoly2 = [sp1, sp2, sp3, \
+south1 = [sp1, sp2, sp3, \
              sp16, sp17, sp18, sp19, sp20, sp21, sp22, sp23, \
              sp4, sp5, sp6, sp7, sp8, sp9, sp10, \
@@ -446 +460 @@
 #slump_origin = [x0+314036.58727982, y0+6224951.2960092] #Absolute UTM
 slump_origin = [385000.0, 6255000.0] #Absolute UTM
+
+a = [340000, 6255000]
+b = [340000, 6270000]
+c = [318000, 6255000]
+d = [318000, 6270000]
+e = [316000, 6255000]
+f = [395000, 6255000]
+g = [355000, 6280000]
+h = [355000, 6224000]
+test_pts = [a, b, c, d, e, f, g, h]
+test_elev = [1.0, 4.0, 3.0, 0.1, 5, -10.0, -20, -15]