Changeset 668 for inundation/ga/storm_surge/examples/beach.py

Timestamp:

Dec 2, 2004, 11:40:39 PM (20 years ago)

Author:

ole

Message:

 

File:

• inundation/ga/storm_surge/examples/beach.py (modified) (7 diffs)

inundation/ga/storm_surge/examples/beach.py

@@ -17 +17 @@
 
 from shallow_water import Domain, Reflective_boundary, Dirichlet_boundary,\
-     Transmissive_boundary, Time_boundary
+     Transmissive_boundary, Time_boundary, Wind_stress
 
 from pmesh2domain import pmesh_to_domain_instance
@@ -38 +38 @@
 
 domain = Domain(domain.coordinates, domain.triangles, domain.boundary)
-print "Number of triangles = ", len(domain)
 
 domain.store = True
@@ -59 +58 @@
 
 def riverbed(x, y):
-    return (y-100)/10 + 4*x/100
+    #return (y-100)/10 + 5*x/100
+    return (y-100)/50 + 4.4
     
 
@@ -79 +79 @@
 
 river = [[20,100], [18,90], [20,80], [20,60], [15,40], [11,20], [2,0], [10,0],
-         [14,10], [20,30], [24,45], [27,80], [27,85], [32,90], [39,100]]
+         [14,10], [20,30], [24,45], [27,80], [27,85], [35,90], [39,100]]
               
 print 'Set elevation'    
 domain.set_quantity('elevation',
     Polygon_function( [(shoreline, bathymetry), (land, topography), 
-                      (building1, 4), (building2, 11), 
-                      (building3, 7), (building4, 8),
-                      (building5, 8), (building6, 7),
+                      (building1, 3), (building2, 9), 
+                      (building3, 8), (building4, 7),
+                      (building5, 8), (building6, 6),
                       (river, riverbed)]))
                       
 print 'Set level'                     
 domain.set_quantity('level', 
-                    Polygon_function( [(shoreline, -1.2), 
+                    Polygon_function( [(shoreline, -1.5), 
                     (land, -10)] )) 
                     
 print 'Set friction'                                        
-domain.set_quantity('friction', 0.07)
+domain.set_quantity('friction', 0.03)
 
 print domain.get_extent()
+
+
+
+#Add lateral wind gusts bearing 135 degrees
+def gust(t,x,y): 
+    from math import sin, pi
+    from Numeric import zeros, ones, Float
+
+    N = len(x)
+
+    tt = sin(2*pi*t/50)
+
+    if tt > 0.7:
+        return 10000*tt*ones(N, Float)
+    else:
+        return zeros(N, Float)
+    
+domain.forcing_terms.append(Wind_stress(gust, 135))
+
 
 ######################
@@ -104 +123 @@
 print 'Boundaries'
 Br = Reflective_boundary(domain)
-Bt = Transmissive_boundary(domain)
+Bo = Transmissive_boundary(domain)
 
 #Constant inflow
 Bd = Dirichlet_boundary([0.0, -1.0, 0.0])
-Bt = Time_boundary(domain, lambda t: [ 2.7*(1+sin(2*pi*t/50)), 0.0, 0.0])
+Bt = Time_boundary(domain, lambda t: [ 3.0*(1+sin(2*pi*t/100)), 0.0, 0.0])
 
 print 'Available boundary tags are', domain.get_boundary_tags()
@@ -116 +135 @@
 tags['ocean'] = Bt
 tags['wall'] = Br 
+tags['outflow'] = Bo 
 tags['exterior'] = Br
 tags['external'] = Br 
@@ -125 +145 @@
 ######################
 #Evolution
-for t in domain.evolve(yieldstep = 0.5, finaltime = 100):
+for t in domain.evolve(yieldstep = 0.2, finaltime = 500):
     domain.write_time()