Changeset 519

Timestamp:

Nov 10, 2004, 5:24:16 PM (20 years ago)

Author:

ole

Message:

Committed wind field with examples.
This may need a C-implementation for speed.

Location:

inundation/ga/storm_surge/pyvolution

Files:

• shallow_water.py (modified) (1 diff)
• wind_example.py (added)
• wind_example_rot.py (added)

inundation/ga/storm_surge/pyvolution/shallow_water.py

@@ -963 +963 @@
             ymom_update += S*v        
         
-        
-def wind_stress(domain):
-    """Apply wind stress to water momentum
-    """
-
-    #FIXME: Under construction
-
-    #Hardwired velocity field - should be initialised as a
-    #callable object
-
-    u = -10000.0
-    v =  10000.0
-
-    Cw = 3.0e-3 #Wind stress coefficient
-    rho_a = 1.2e-3 #Atmospheric density
-    rho = 1023 #Density of water
-    
-    from math import sqrt
-
-    xmom_update = domain.quantities['xmomentum'].explicit_update
-    ymom_update = domain.quantities['ymomentum'].explicit_update    
-    
-    N = domain.number_of_elements
-
-    c = Cw*rho_a/rho
-    S = c * sqrt(u**2 + v**2)
-
-    #print 'B4:', xmom_update[10]
-    
-    xmom_update += S*u
-    ymom_update += S*v
-
-    #print 'After:', xmom_update[10]    
-