Changeset 4258

Timestamp:

Feb 14, 2007, 3:03:11 PM (18 years ago)

Author:

ole

Message:

Added set_minimum_allowed_height and updated user manual.
Played with Brad's dam break example.

Files:

• anuga_core/documentation/user_manual/anuga_user_manual.tex (modified) (1 diff)
• anuga_core/source/anuga/shallow_water/shallow_water_domain.py (modified) (2 diffs)
• anuga_work/development/dam_test_from_brad_2007/dam_sample.py (modified) (2 diffs)

anuga_core/documentation/user_manual/anuga_user_manual.tex

@@ -1726 +1726 @@
     \code{config.py}.
 \end{methoddesc}
+
+
+\begin{methoddesc} {set\_minimum_allowed_height}{}
+    Module: \module{shallow\_water.shallow\_water\_domain}
+
+    Set the minimum depth (in meters) that will be recognised in 
+    the numerical scheme (including limiters and flux computations)
+    
+    Default value is $10^{-3}$ m, but by setting this to a greater value, 
+    e.g.\ for large scale simulations, the computation time can be 
+    significantly reduced. 
+\end{methoddesc}
+
 
 \begin{methoddesc} {set\_minimum_storable_height}{}

anuga_core/source/anuga/shallow_water/shallow_water_domain.py

@@ -146 +146 @@
                                 number_of_full_triangles=number_of_full_triangles) 
 
-        self.minimum_allowed_height = minimum_allowed_height
-        self.H0 = minimum_allowed_height  # Minimal height for flux computation        
+        #self.minimum_allowed_height = minimum_allowed_height
+        #self.H0 = minimum_allowed_height
+        self.set_minimum_allowed_height(minimum_allowed_height)
+        
         self.maximum_allowed_speed = maximum_allowed_speed
         self.g = g
@@ -214 +216 @@
         """
         self.minimum_storable_height = minimum_storable_height
+
+
+    def set_minimum_allowed_height(self, minimum_allowed_height):
+        """
+        Set the minimum depth that will be recognised in the numerical
+        scheme
+
+        The minimum allowed depth is in meters.
+
+        The parameter H0 (Minimal height for flux computation)
+        is also set by this function
+        """
+        
+        self.minimum_allowed_height = minimum_allowed_height
+        self.H0 = minimum_allowed_height   
         
 

anuga_work/development/dam_test_from_brad_2007/dam_sample.py

@@ -28 +28 @@
 domain = Domain(points, vertices, boundary)   # Create domain
 domain.set_name('dam_sample')                 # Output name
-domain.set_store_vertices_uniquely(True)      # Look at the 'real' triangles
-domain.beta_h = 0
+#domain.set_store_vertices_uniquely(True)      # Look at the 'real' triangles
+domain.beta_h = 0.0
 domain.limit2007 = 1
+domain.set_minimum_allowed_height(0.01)
 
 #------------------------------------------------------------------------------
@@ -92 +93 @@
 #------------------------------------------------------------------------------
 for t in domain.evolve(yieldstep = 0.2, finaltime = 20):
-    domain.write_time(track_speeds=True)
+    domain.write_time(track_speeds=False)