Changeset 4218 for anuga_core/source/anuga/shallow_water/shallow_water_domain.py

Timestamp:

Feb 6, 2007, 11:47:41 AM (18 years ago)

Author:

ole

Message:

Implemented new Flux Limiter which I believe is more robust and will pave the way for new balanced slope limiters.
For backwards compatibility using fluxes as they were before today (5/2/7) set
domain.H0 = 0 instead of the new default which is minimum_allowable_height (0.001).

File:

• anuga_core/source/anuga/shallow_water/shallow_water_domain.py (modified) (3 diffs)

anuga_core/source/anuga/shallow_water/shallow_water_domain.py

@@ -147 +147 @@
 
         self.minimum_allowed_height = minimum_allowed_height
+        self.H0 = minimum_allowed_height  # Minimal height for flux computation        
         self.maximum_allowed_speed = maximum_allowed_speed
         self.g = g
@@ -832 +833 @@
          compute_fluxes_ext_central as compute_fluxes_ext
     
-    domain.timestep = compute_fluxes_ext(timestep, domain.epsilon, domain.g,
+
+    domain.timestep = compute_fluxes_ext(timestep, domain.epsilon,
+                                         domain.H0,                                         
+                                         domain.g,
                                          domain.neighbours,
                                          domain.neighbour_edges,
@@ -2082 +2086 @@
                     protect_against_infinitesimal_and_negative_heights_c
 
-
     #distribute_to_vertices_and_edges =\
     #              distribute_to_vertices_and_edges_c #(like MH's)