Changeset 2317 for production/sydney_2006/run_sydney_smf.py

Timestamp:

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

Author:

sexton

Message:

Updates to Sydney example.

File:

• production/sydney_2006/run_sydney_smf.py (modified) (3 diffs)

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')