Changeset 5176

Timestamp:

Mar 19, 2008, 9:02:27 PM (16 years ago)

Author:

steve

Message:

Edge limiter needed protection from small values of h

Files:

• anuga_core/source/anuga/shallow_water/shallow_water_domain.py (modified) (1 diff)
• anuga_work/development/attenuation_near_shore/run_beach_structured.py (modified) (3 diffs)

anuga_core/source/anuga/shallow_water/shallow_water_domain.py

@@ -396 +396 @@
         # Call correct module function
         # (either from this module or C-extension)
-        if self.use_edge_limiter:            
+        if self.use_edge_limiter:
+            protect_against_infinitesimal_and_negative_heights(self)
             for name in self.conserved_quantities:
                 Q = self.quantities[name]

anuga_work/development/attenuation_near_shore/run_beach_structured.py

@@ -79 +79 @@
         thinner=True
     else:
-        yieldstep = 10
+        yieldstep = 1.0
         finaltime = 1600
-        finaltime = 400
+        finaltime = 1000
         maximum_triangle_area = 100
         thinner=True
@@ -143 +143 @@
     domain.set_timestepping_method('rk2')
     domain.set_default_order(2)
-    domain.use_old_limiter = False
+    domain.use_edge_limiter = True
     domain.tight_slope_limiters = True
     
@@ -173 +173 @@
     Bd = Dirichlet_boundary([10.,0,0])
     Bwp = Transmissive_Momentum_Set_Stage_boundary(domain,
-                   lambda t: [(1.2*sin(2*t*pi/50)), 0.0 ,0.0])
+                   lambda t: [(1.2*sin(2*t*pi/500)), 0.0 ,0.0])
     Bwp_velocity = Time_boundary(domain,
                    lambda t: [(1.2*sin(2*t*pi/10)),
                               0.24*sin(2*t*pi/10),0.0])
     Bws = Transmissive_Momentum_Set_Stage_boundary(domain,
-                   lambda t: [1.2*sin(2*t*pi/25),
+                   lambda t: [1.2*sin(2*t*pi/100),
                               0.0, 0.0])
     #domain.set_boundary( {'wall': Br, 'wave': Bwp_velocity} )