Changeset 4243 for anuga_core/source/anuga

Timestamp:

Feb 8, 2007, 4:43:35 AM (18 years ago)

Author:

ole

Message:

Tried new balanced limiter with examples

Location:

anuga_core/source/anuga/examples

Files:

• island.py (modified) (4 diffs)
• limiter_test.py (modified) (1 diff)

anuga_core/source/anuga/examples/island.py

@@ -33 +33 @@
                                            'top': [2],
                                            'left': [3]},
-                          maximum_triangle_area = 1.0,
+                          maximum_triangle_area = 10.0,
                           filename = 'island.msh' ,
-                          interior_regions=[ ([[50,25], [70,25], [70,75], [50,75]], 1.0)]
+                          interior_regions=[ ([[50,25], [70,25], [70,75], [50,75]], 10.0)]
                           #interior_holes=[[[50,25], [70,25], [70,75], [50,75]]],
                           )
@@ -55 +55 @@
 # beta_h == 0.0 means that constant (1st order) gradients are introduced
 # on h. This is equivalent to the constant depth used previously.
-domain.beta_h     = 0.5
-domain.beta_w_dry = 0.0
-domain.alpha_balance = 10.0
-domain.minimum_allowed_height = 1.0e-4 
-domain.maximum_allowed_speed = 100.0 
-domain.minimum_storable_height = 1.0e-4 
+#domain.beta_h     = 0.5
+#domain.beta_w_dry = 0.0
+#domain.alpha_balance = 10.0
+#domain.minimum_allowed_height = 1.0e-4 
+#domain.maximum_allowed_speed = 100.0 
+#domain.minimum_storable_height = 1.0e-4 
+
+domain.beta_h     = 0.0
+domain.limit2007 = 1
 
 #------------------------------------------------------------------------------
@@ -90 +93 @@
 
 
+
 #------------------------------------------------------------------------------
 # Setup boundary conditions (all reflective)
@@ -105 +109 @@
 
 import time
-for t in domain.evolve(yieldstep = 1, finaltime = 100):
+for t in domain.evolve(yieldstep = 1, finaltime = 10):
     domain.write_time()
     #if allclose(t, 100):

anuga_core/source/anuga/examples/limiter_test.py

@@ -27 +27 @@
 domain.set_store_vertices_uniquely(True)
 domain.set_default_order(2)
+domain.limit2007 = 1
 
 #------------------------------------------------------------------------------