Changeset 704

Timestamp:

Dec 14, 2004, 10:21:40 AM (19 years ago)

Author:

ole

Message:

Added wind field to Merimbula script

Location:

inundation/ga/storm_surge/validation

Files:

• README.txt (added)
• run_merimbula.py (modified) (3 diffs)
• tide_example.txt (added)
• windstress_example.txt (added)

inundation/ga/storm_surge/validation/run_merimbula.py

@@ -13 +13 @@
      Dirichlet_boundary, Wind_stress
 from pmesh2domain import pmesh_to_domain_instance
+from util import File_function
 
 ######################
@@ -32 +33 @@
 
 
-#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))
+#Add time dependent wind field from file
+#Format is time [DD/MM/YY hh:mm:ss], speed [m/s] direction (degrees)
+#See also README.txt
+F = File_function('windstress_example.txt', domain)
+domain.forcing_terms.append(Wind_stress(F))
 
 
@@ -52 +43 @@
 # Boundary conditions
 
-filename = 'Eden_Australia_31082004.txt'
-Bf = File_boundary(filename, domain)
-
-inflow_stage = 20.0
-Bd = Dirichlet_boundary([inflow_stage, 0.0, 0.0])
+Bf = File_boundary('tide_example.txt', domain)
 Br = Reflective_boundary(domain)
 domain.set_boundary({'external': Br, 'open': Bf})