Changeset 5176 for anuga_work/development/attenuation_near_shore

Timestamp:

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

Author:

steve

Message:

Edge limiter needed protection from small values of h

File:

• anuga_work/development/attenuation_near_shore/run_beach_structured.py (modified) (3 diffs)

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} )