Changeset 694 for inundation/ga/storm_surge/validation

Timestamp:

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

Author:

ole

Message:

 

File:

• inundation/ga/storm_surge/validation/run_merimbula.py (modified) (3 diffs)

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