Changeset 2875

Timestamp:

May 16, 2006, 1:35:21 PM (19 years ago)

Author:

sexton

Message:

sydney - export results

Location:

production/sydney_2006

Files:

• export_results.py (modified) (1 diff)
• project.py (modified) (15 diffs)

production/sydney_2006/export_results.py

@@ -4 +4 @@
 from pyvolution.ermapper_grids import read_ermapper_grid
 
-#name = project.newoutputname
-name = project.outputname
+name = project.outputname2
 
 #print 'Which variable do you want to export?'
 #which_var = int(raw_input('Stage = 0, Absolute Momentum = 1, Depth = 2, Speed = 3  '  ))
-which_var = 3
+which_var = 2
 
 if which_var == 0:  # Stage

production/sydney_2006/project.py

@@ -30 +30 @@
 nmaxviz = 6283000
 
-basename = 'slump_friction04'
+basename = 'slump_315res'
+basename2 = 'slump_1000res'
 basename4 = 'slump_poly_ingo_test'
 
@@ -44 +45 @@
 polygondir = home+sep+scenario_dir_name+sep+'polygons'+sep
 gaugedir = home+sep+scenario_dir_name+sep+'gauges'+sep
+#reportdir = home+sep+scenario_dir_name+sep+'reports'+sep
 
 meshname = meshdir + basename
+meshelevname = meshdir + 'test_elev.tsh'
 meshname4 = meshdir + basename4
 coarsedemname = datadir + coarsename
@@ -51 +54 @@
 combineddemname = datadir + 'sydneytopo'
 outputname = outputdir + basename  #Used by post processing
-outputname4 = outputdir + basename4  #Used by post processing
+outputname2 = outputdir+sep+'Coast_Polygons'+sep+'res1000'+sep+basename2  #Used by post processing
 
 #csv file of coastline 50m epsilon belt
@@ -59 +62 @@
 
 gauge_filename = gaugedir + 'sydney_gauges.xya'
+#gauge_filename = gaugedir + 'sydney_gauges_test.xya'
 gauge_outname = gaugedir + 'gauges_max_output.xya'
 #gauge_filename = gaugedir + 'nest_gauges_Manly.xya'
@@ -64 +68 @@
 #gauge_filename = gaugedir + 'GA_gauge.csv' # from Benfield
 #gauge_outname = gaugedir + 'gauges_max_output_next.xya'
-gaugetimeseries = gaugedir + 'gauges_time_series_Benfield'
+gaugetimeseries = gaugedir + 'gauges_time_series'
 polygonptsfile = polygondir + 'poly'
 integraltimeseries = outputdir + 'integral_time_series'
+integraltimeseries2 = outputdir + 'integral_time_series_move_origin'
+
 
 #Georeferencing
@@ -86 +92 @@
 print 'Got refzone', refzone
 
+# original clipping
 dsouth = degminsec2decimal_degrees(-34,05,0)
+dsouth2 = degminsec2decimal_degrees(-34,01,0)
 dnorth = degminsec2decimal_degrees(-33,33,0)
 dnorth1 = degminsec2decimal_degrees(-33,40,0)
 dnorth2 = degminsec2decimal_degrees(-33,58,30)
 dnorth3 = degminsec2decimal_degrees(-33,46,0)
-dwest = degminsec2decimal_degrees(151,2,20)
+dwest = degminsec2decimal_degrees(151,1,0)
 deast1 = degminsec2decimal_degrees(151,20,0)
 deast2 = degminsec2decimal_degrees(151,48,0)
@@ -106 +114 @@
 dp7 = [dnorth3, dwest]
 dp8 = [dnorth, dwest]
+
+dp12 = [dsouth2, deast2]
+dp34 = [dnorth, deast2]
+diffpolygonall, zone = convert_points_from_latlon_to_utm([dp0, dp1, dp2, dp3, dp4, dp5, dp6, dp7])
+# used for new tests 4 April 2006 (ensure slump contained in domain)
+diffpolygonall2, zone = convert_points_from_latlon_to_utm([dp0, dp1, dp12, dp34, dp4, dp5, dp6, dp7])
+
+# clipping used for look at increasingly finer resolution
+j0 = [338000, 6243000]
+j1 = [365000, 6243000]
+j2 = [365000, 6273000]
+j3 = [338000, 6273000]
     
-diffpolygonall, zone = convert_points_from_latlon_to_utm([dp0, dp1, dp2, dp3, dp4, dp5, dp6, dp7])
+diffpolygonall_test = [j0, j1, j2, j3]
+# clipping used for further investigation of black screen of death
+j0 = [328000, 6255000]
+j1 = [355000, 6255000]
+j2 = [355000, 6270000]
+j3 = [328000, 6270000]
+
+diffpolygonall_test2 = [j0, j1, j2, j3]
+
+# clipping used for demo purposes to fit around poly3 and poly4 (Ingo's files)
+j0 = [318000, 6249000]
+j1 = [387000, 6249000]
+j2 = [387000, 6270000]
+j3 = [318000, 6270000]
+
+demopoly = [j0, j1, j2, j3]
+
 # to put chunk back in
 #diffpolygonall = [dp0, dp1, dp2, dp3, dp4, dp8]
@@ -141 +177 @@
 poly2 = [qp0, qp1, qp2, qp3, qp4, qp5, qp6, qp7, qp8, qp9, qp10, qp11, qp12]
 #poly2 = [qp5, qp6, qp7, qp8, qp9, qp10]
-
-# didn't like this one - poly2 from Ingo
-# warning generated about vertex blah doesn't belong to an element
-# is there a maximum number of vertices per polygon?
-p0 = [343488.9455,      6270644.956]
-p1 = [343681.208,       6270500.759]
-p2 = [343997.0678,      6270253.564]
-p3 = [343482.0789,      6269930.838]
-p4 = [343550.7441,      6269546.313]
-p5 = [343111.287,       6269099.989]
-p6 = [342630.6307,      6267534.423]
-p7 = [342870.9588,      6266380.848]
-p8 = [342479.5673,      6266126.787]
-p9 = [342174.4759,      6266317.303]
-p10 = [342248.668,      6266727.125]
-p11 = [342101.9088,     6267108.699]
-p12 = [340986.2201,     6267454.634]
-p13 = [340663.6936,     6266113.602]
-p14 = [340120.4909,     6265892.926]
-p15 = [339560.3132,     6266317.303]
-p16 = [339237.7867,     6266928.406]
-p17 = [339000.1355,     6267912.961]
-p18 = [339169.8863,     6268642.889]
-p19 = [340001.6654,     6269016.341]
-p20 = [340392.0922,     6269423.743]
-p21 = [340867.3945,     6268863.566]
-p22 = [341614.2982,     6268659.865]
-p23 = [341495.4726,     6269220.042]
-p24 = [341750.0988,     6269287.943]
-p25 = [341902.8746,     6269882.07]
-p26 = [342700.7034,     6269780.22]
-p27 = [343125.0805,     6269746.27]
-p28 = [343125.0805,     6270221.572]
-p29 = [343057.1801,     6270527.124]
-p30 = [343488.9455,     6270644.956]
-
-testpoly = [p0,p1,p2,p3,p4,p5,p6,p7,p8,p9,p10,p11,p12,p13,p14,p15,p16,p17,p18,p19,p20,p21,p22,p23,p24,p25,p26,p27,p28,p29,p30]
 
 # test from Jane drawing
@@ -194 +193 @@
 np9 = [343000, 6264000]
 np10 = [345000, 6265000]
+np10_2 = [344500, 6265000]
+np10_3 = [344200, 6265000]
 np11 = [343000, 6266000]
 np12 = [344000, 6269000]
@@ -225 +226 @@
 np38 = [324000, 6253000]
 np39 = [325000, 6251000]
+nptest = [339000, 6266000]
+nptest2 = [344000, 6265000]
 
 #testpoly1 = [np1, \
@@ -239 +242 @@
 #             np31, np32, \
 #             np23]
+
+# reduce this polygon for black screen of death testing 24/04/06
+#newpoly1_refine = [nptest, np27, np16] # worked
+#newpoly1_refine = [np9, np10, np11] # worked
+#newpoly1_refine = [np9, np10, np11, nptest2] # worked
+#newpoly1_refine = [np9, np10_2, np11, nptest2] # worked
+newpoly1_refine = [np9, np10_3, np11, nptest2] # worked
+
+# used for refined interior regions and fine mesh
 newpoly1 = [np28, np29, np30, \
             np5, np6, np7, np8, np9, np10, \
@@ -246 +258 @@
             np19, np20, np21, np22, np31, np32, \
             np36, np37, np38, np39, np2, np3]
+
 # last two lines for second run
 
@@ -280 +293 @@
 sp23 = [339000, 6236000]
 
+#original
 south1 = [sp1, sp2, sp3, \
              sp16, sp17, sp18, sp19, sp20, sp21, sp22, sp23, \
              sp4, sp5, sp6, sp7, sp8, sp9, sp10, \
-             sp11, sp12, sp13, sp14, sp15]
+             sp11, sp12, sp13, sp14, sp15] 
 
 m1 = [340000, 6256000]
@@ -298 +312 @@
     
 finepolyquay = [q1, q2, q3, q4]
+
+# refined 13/03/06 to look at effect of finer resolution on no IC scenario
+q2south = [340000, 6249000]
+newpoly1 = [np28, np29, q2south, sp7, sp8, sp5, sp6, np30, \
+            np5, np6, np7, np8, np9, np10, \
+            np11, np12, np13, np14, np15, \
+            np27, np26, np25, np24, np17, \
+            np33, \
+            np19, np20, np21, np22, np31, np32, \
+            np36, np37, np38, np39, np2, np3]
 
 #Interior regions - the Harbour
@@ -460 +484 @@
 #slump_origin = [x0+314036.58727982, y0+6224951.2960092] #Absolute UTM
 slump_origin = [385000.0, 6255000.0] #Absolute UTM
+# move 10km west
+slump_origin2 = [375000.0, 6255000.0] #Absolute UTM
 
 a = [340000, 6255000]