Changeset 2188 for production/sydney_2006/run_sydney.py

Timestamp:

Jan 6, 2006, 1:54:51 PM (19 years ago)

Author:

sexton

Message:

update Sydney scenario to allow switching between clipping regions. BC = 10min square wave at 6m height on right hand boundary.

File:

• production/sydney_2006/run_sydney.py (modified) (9 diffs)

production/sydney_2006/run_sydney.py

@@ -8 +8 @@
 the elevation data and boundary data obtained from a tsunami simulation done with MOST.
 
-Ole Nielsen, GA - 2005
+Ole Nielsen, GA - 2005 and Jane Sexton, GA - 2006
 """
 
@@ -40 +40 @@
       verbose = True)
 
-#this was going to switch between different clipping polygons but doesn't
-#seem to work, so keep this at zero
-mytest = 0
+#this allows the user to switch between different clipping polygons
+print 'Which total zone are you interested in?'
+mytest = int(raw_input('0 = all, 1 = harbour and 2 = botany bay     '))
 
 #Create Triangular Mesh
@@ -63 +63 @@
                'resolution': 100000,
                'filename': meshname,           
-               'interior_regions': interior_regions},     
+               'interior_regions': interior_regions},
+              evaluate=True,    
               verbose = True)
+    #import sys; sys.exit()
 
 if mytest == 1:
@@ -73 +75 @@
     east = project.heast
 
-    interior_regions = [[project.harbour_polygon, 25000]] # maximal area of per triangle
+    #interior_regions = [[project.harbour_polygon, 25000]] # maximal area of per triangle
                        
     m = cache(create_mesh_from_regions,
@@ -79 +81 @@
               {'boundary_tags': {'bottom': [0], 'top': [2],
                                  'right': [1], 'left': [3]},
-               'resolution': 100000,
-               'filename': meshname,           
-               'interior_regions': interior_regions},     
+               'resolution': 50000,
+               'filename': meshname},           
+     #          'interior_regions': interior_regions},
+              evaluate=True, 
               verbose = True)
 
@@ -91 +94 @@
     east = project.beast
 
-    interior_regions = [[project.botanybay_polygon, 25000]] # maximal area of per triangle
+    #interior_regions = [[project.botanybay_polygon, 25000]] # maximal area of per triangle
                        
     m = cache(create_mesh_from_regions,
@@ -97 +100 @@
               {'boundary_tags': {'bottom': [0], 'top': [2],
                                  'right': [1], 'left': [3]},
-               'resolution': 100000,
-               'filename': meshname,           
-               'interior_regions': interior_regions},     
+               'resolution': 50000,
+               'filename': meshname},
+     #          'interior_regions': interior_regions},
+              evaluate=True, 
               verbose = True)
     
@@ -137 +141 @@
 Bt = Transmissive_boundary(domain)
 Bd = Dirichlet_boundary([0,0,0])
+# 10 min square wave starting at 1 min, 6m high
 Bw = Time_boundary(domain=domain,
-                   f=lambda t: [(6<t<30)*5000, 0, 0])
+                   f=lambda t: [(60<t<660)*6, 0, 0])
 
 
@@ -148 +153 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 1, finaltime = 50):
+for t in domain.evolve(yieldstep = 1, finaltime = 3600):
     domain.write_time()
     domain.write_boundary_statistics(tags = 'right') #quantities = 'stage')