Changeset 907 for inundation/ga/storm_surge/pyvolution/quantity.py

Timestamp:

Feb 17, 2005, 5:36:29 PM (20 years ago)

Author:

ole

Message:

Addressed problem with artificial momentum generated by discontinuous water depths in the presence of steep slopes.
Now very shallow water is limited with a separate h-limiter controlled by beta_h
(see config.py) for details.

File:

• inundation/ga/storm_surge/pyvolution/quantity.py (modified) (2 diffs)

inundation/ga/storm_surge/pyvolution/quantity.py

@@ -813 +813 @@
     N = quantity.domain.number_of_elements
     
-    beta = quantity.domain.beta
+    beta_w = quantity.domain.beta_w
         
     qc = quantity.centroid_values
@@ -852 +852 @@
             if (dq[k,i] < 0): r = dqmin[k]/dq[k,i]
             
-            phi = min( min(r*beta, 1), phi )    
+            phi = min( min(r*beta_w, 1), phi )    
 
         #Then update using phi limiter