Changeset 6647

Timestamp:

Mar 27, 2009, 3:32:36 PM (15 years ago)

Author:

ted

Message:

Started function for volumetric test

File:

• anuga_core/source/anuga/shallow_water/shallow_water_domain.py (modified) (2 diffs)

anuga_core/source/anuga/shallow_water/shallow_water_domain.py

@@ -781 +781 @@
         
 
+                def compute_volumetric_balance(self):
+                        """Compute volumetric balance at current timestep.
+                        
+                        Quantities computed are:
+                                Total inflow across boundary
+                                Total outflow across boundary
+                                Total inflow through forcing terms
+                                Total outflow through forcing terms
+                                Current total volume in domain
+                                                                
+                        
+                        """
+                
+                        # Run through boundary array and compute for each segment
+                        # the normal momentum ((uh, vh) dot normal) times segment length.
+                        # Based on sign accumulate this into total_boundary_inflow and total_boundary_outflow.
+                        # The go through explicit forcing update and record the rate of change for stage and 
+                        # record into total_forcing_inflow and total_forcing_outflow. Finally compute integral 
+                        # of depth to obtain total volume of domain.
+                        
+                        
+                        
+                        pass
+                
+                
 #=============== End of class Shallow Water Domain ===============================
 
@@ -1559 +1584 @@
     
     return wet_friction
-
-
-
-
-
-