Changeset 2317

Timestamp:

Feb 1, 2006, 4:28:22 PM (19 years ago)

Author:

sexton

Message:

Updates to Sydney example.

Location:

production/sydney_2006

Files:

• export_results.py (modified) (2 diffs)
• project.py (modified) (4 diffs)
• run_sydney_smf.py (modified) (3 diffs)

production/sydney_2006/export_results.py

@@ -8 +8 @@
 
 #print 'Which variable do you want to export?'
-#which_var = int(raw_input('Stage = 0, Absolute Momentum = 1, Depth = 2'  ))
-which_var = 2
+#which_var = int(raw_input('Stage = 0, Absolute Momentum = 1, Depth = 2, Speed = 3  '  ))
+which_var = 3
 
 if which_var == 0:  # Stage
@@ -22 +22 @@
     outname = name + '_depth'
     quantityname = 'stage-elevation'  #Depth
+
+if which_var == 3:  # Speed
+    outname = name + '_speed'
+    #quantityname = '((xmomentum/(stage-elevation))**2 + (ymomentum/(stage-elevation))**2)**0.5'  #Speed
+    #quantityname = 'xmomentum/(stage-elevation)'  #Speed
+    quantityname = '(xmomentum**2 + ymomentum**2 + 1.e-30)**0.5/(stage-elevation)'  #Speed
+
     
 sww2dem(name, basename_out = outname,

production/sydney_2006/project.py

@@ -37 +37 @@
 nminviz = 6231000
 nmaxviz = 6283000
-   
+# for harbour region
+#nminviz = 6248000
+#nmaxvix = 6270000
+# for botany bay region
+#nmaxviz = 6248000
+#nminvix = 6283000
+
 basename = 'slump'
 
@@ -60 +66 @@
 #newoutputname = outputdir+'Combined Data v1'+sep+'d400'+sep+basename
 gauge_filename = outputdir + 'sydney_gauges.xya'
+gauge_outname = outputdir + 'gauges_max_output.xya'
 
 #Georeferencing
@@ -131 +138 @@
 harbour_7x = degminsec2decimal_degrees(-33,34,30)
 harbour_7y = degminsec2decimal_degrees(151,20,20)
+#harbour_8x = degminsec2decimal_degrees(-33,33,30)
+#harbour_8y = degminsec2decimal_degrees(151,15,0)
 harbour_8x = degminsec2decimal_degrees(-33,33,30)
-harbour_8y = degminsec2decimal_degrees(151,15,0)
+harbour_8y = degminsec2decimal_degrees(151,17,0)
+#harbour_9x = degminsec2decimal_degrees(-33,45,30)
+#harbour_9y = degminsec2decimal_degrees(151,15,0)
 harbour_9x = degminsec2decimal_degrees(-33,45,30)
-harbour_9y = degminsec2decimal_degrees(151,15,0)
+harbour_9y = degminsec2decimal_degrees(151,17,0)
 harbour_10x = degminsec2decimal_degrees(-33,45,10)
 harbour_10y = degminsec2decimal_degrees(151,11,40)
@@ -197 +208 @@
 j92 = [bb_10x, bb_10y]
 
-botanybay_polygon_2 = [j92, j12, j22, j62, j82, j72, j42, j02]
+botanybay_polygon_2 = [j92, j12, j22, j62, j82, j72, j42] # worked
 #botanybay_polygon_2 = [j02, j12, j22, j32, j42, j52] # worked
 

production/sydney_2006/run_sydney_smf.py

@@ -80 +80 @@
 if mytest == 0:
     # for whole region
-    interior_res = 2000 # run at 2000 for final one
+    interior_res = 5000 
     interior_regions = [[project.harbour_polygon_2, interior_res],
                         [project.botanybay_polygon_2, interior_res]] # maximal area of per triangle
@@ -134 +134 @@
     # as advised by Adrian, this can be as much as 15000-30000 metres
     # len = 4500.0
-    len = 15000.0
+    len = 30000.0
     #thk = 670 #thk = 760.0
-    thk = 75.0
+    dep = 400.0
+    thk = 0.44*dep #thk = 120.0
     wid = 15000.0 # for slump scenario, wid=len
     #dep = 400 #1000.0 #d = 400 doesn't make sense with thk = 670!
-    dep = 400.0
     rad = 3330
     dp = 0.23
@@ -190 +190 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 120, finaltime = 10800): #120 10800 = 3hrs
+for t in domain.evolve(yieldstep = 120, finaltime = 18000): #120 10800 = 3hrs, 18000 = 5hrs
     domain.write_time()
     domain.write_boundary_statistics(tags = 'bottom') #quantities = 'stage')