Changeset 694

Timestamp:

Dec 8, 2004, 2:29:01 PM (20 years ago)

Author:

ole

Message:

 

Location:

inundation/ga/storm_surge

Files:

• examples/beach.py (modified) (3 diffs)
• validation/run_merimbula.py (modified) (3 diffs)

inundation/ga/storm_surge/examples/beach.py ¶

@@ -40 +40 @@
 
 domain.store = True
-domain.set_name(name + '5')
+domain.set_name(name + '6')
 domain.default_order = 2
 print "Output being written to " + data_dir + sep + \
@@ -149 +149 @@
     N = len(x)
 
-    tt = sin(2*pi*t/80)
+    tt = sin(2*pi*t/100)
 
     if tt > 0.95:
@@ -157 +157 @@
     
 domain.forcing_terms.append(Wind_stress(gust2, 90))
+
+#Add lateral wind gusts bearing 255 degrees
+def gust3(t,x,y): 
+    from math import sin, pi
+    from Numeric import zeros, ones, Float
+
+    N = len(x)
+
+    tt = sin(2*pi*(t-30)/55)
+
+    if tt > 0.96:
+        return 24000*tt*ones(N, Float)
+    else:
+        return zeros(N, Float)
+    
+domain.forcing_terms.append(Wind_stress(gust3, 255))
 
 

inundation/ga/storm_surge/validation/run_merimbula.py ¶

@@ -11 +11 @@
 
 from shallow_water import Domain, Reflective_boundary, File_boundary,\
-     Dirichlet_boundary
+     Dirichlet_boundary, Wind_stress
 from pmesh2domain import pmesh_to_domain_instance
 
@@ -24 +24 @@
 print "Number of triangles = ", len(domain)
 
-domain.default_order = 1
-domain.filename = filename
+domain.default_order = 2
+domain.store = True
+domain.set_name('merimbula')
 
 domain.set_quantity('friction', 0.07)
+domain.set_quantity('level', 0.5)
+
+
+#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/30)
+
+    if tt > 0.98:
+        return 24000*tt*ones(N, Float)
+    else:
+        return zeros(N, Float)
+    
+domain.forcing_terms.append(Wind_stress(gust, 135))
+
 
 ######################
@@ -33 +53 @@
 
 filename = 'Eden_Australia_31082004.txt'
-Bf = File_boundary(domain, filename)
+Bf = File_boundary(filename, domain)
 
 inflow_stage = 20.0