Changeset 2186 for production/sydney_2006/run_sydney.py

Timestamp:

Jan 4, 2006, 4:33:02 PM (19 years ago)

Author:

sexton

Message:

Update script for Sydney scenario

File:

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

production/sydney_2006/run_sydney.py

@@ -40 +40 @@
       verbose = True)
 
-south = project.south 
-north = project.north
-west = project.west
-east = project.east
+#this was going to switch between different clipping polygons but doesn't
+#seem to work, so keep this at zero
+mytest = 0
 
 #Create Triangular Mesh
 from pmesh.create_mesh import create_mesh_from_regions
 
-interior_regions = [[project.harbour_polygon, 25000]] # maximal area of per triangle
-print project.harbour_polygon
-                   
-m = cache(create_mesh_from_regions,
-          project.polygon,
-          {'boundary_tags': {'bottom': [0], 'top': [2],
-                             'right': [1], 'left': [3]},
-           'resolution': 100000,
-           'filename': meshname,           
-           'interior_regions': interior_regions},     
-          verbose = True)
+if mytest == 0:
+    # for whole region
+    south = project.south 
+    north = project.north
+    west = project.west
+    east = project.east
 
+    interior_regions = [[project.harbour_polygon, 25000],
+                        [project.botanybay_polygon, 25000]] # maximal area of per triangle
+                       
+    m = cache(create_mesh_from_regions,
+              project.polygonall,
+              {'boundary_tags': {'bottom': [0], 'top': [2],
+                                 'right': [1], 'left': [3]},
+               'resolution': 100000,
+               'filename': meshname,           
+               'interior_regions': interior_regions},     
+              verbose = True)
 
+if mytest == 1:
+    # for harbour region
+    south = project.hsouth 
+    north = project.hnorth
+    west = project.hwest
+    east = project.heast
+
+    interior_regions = [[project.harbour_polygon, 25000]] # maximal area of per triangle
+                       
+    m = cache(create_mesh_from_regions,
+              project.polygon_h,
+              {'boundary_tags': {'bottom': [0], 'top': [2],
+                                 'right': [1], 'left': [3]},
+               'resolution': 100000,
+               'filename': meshname,           
+               'interior_regions': interior_regions},     
+              verbose = True)
+
+if mytest == 2:
+    # for botany bay region
+    south = project.bsouth 
+    north = project.bnorth
+    west = project.bwest
+    east = project.beast
+
+    interior_regions = [[project.botanybay_polygon, 25000]] # maximal area of per triangle
+                       
+    m = cache(create_mesh_from_regions,
+              project.polygon_bb,
+              {'boundary_tags': {'bottom': [0], 'top': [2],
+                                 'right': [1], 'left': [3]},
+               'resolution': 100000,
+               'filename': meshname,           
+               'interior_regions': interior_regions},     
+              verbose = True)
+    
 #Setup domain
 
@@ -97 +138 @@
 Bd = Dirichlet_boundary([0,0,0])
 Bw = Time_boundary(domain=domain,
-                   f=lambda t: [4*sin(2*pi*t/100)+4, 0, 0])
+                   f=lambda t: [(6<t<30)*5000, 0, 0])
 
 
-domain.set_boundary( {'top': Bd, 'bottom': Bd, 'right': Bd, 'left': Bd} )
+domain.set_boundary( {'top': Br, 'bottom': Br, 'right': Bw, 'left': Br} )
 
 
@@ -107 +148 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 1, finaltime = 4):
+for t in domain.evolve(yieldstep = 1, finaltime = 50):
     domain.write_time()
     domain.write_boundary_statistics(tags = 'right') #quantities = 'stage')