Changeset 5439

Timestamp:

Jun 26, 2008, 10:58:08 AM (15 years ago)

Author:

ole

Message:

Updated 'two_levels' example to show that tight_slope_limiters take care of upwards momentum or 'creep'. This used to be treated with the h-limiter with beta_h > 0 but this is no longer needed.

File:

• anuga_work/development/demos/twolevels.py (modified) (2 diffs)

anuga_work/development/demos/twolevels.py

@@ -9 +9 @@
 #
 from os import sep, path
-from mesh_factory import rectangular
-from shallow_water import Domain, Reflective_boundary, Dirichlet_boundary,\
-     Constant_height
+from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular
+from shallow_water import Domain, Reflective_boundary, Dirichlet_boundary
 from Numeric import array
-from anuga.pyvolution.util import Polygon_function, read_polygon
+from anuga.utilities.polygon import Polygon_function, read_polygon
 
 
@@ -33 +32 @@
 #PLAY WITH THIS [0;1]: 
 #
-# beta_h == 0.0 reveals the problem
-# beta_h > 0.2 alleviates it
-domain.beta_h = 0.2    
+# To reveal the problem set
+# domain.tight_slope_limiters = False
+domain.tight_slope_limiters = False
+# To fix it set
+#domain.tight_slope_limiters = True
+domain.tight_slope_limiters = True
+
 
 #IC