Changeset 4702 for anuga_core/source/anuga/abstract_2d_finite_volumes/domain.py

Timestamp:

Sep 5, 2007, 10:06:48 AM (17 years ago)

Author:

ole

Message:

Work on ticket:192 (and untabify)

File:

• anuga_core/source/anuga/abstract_2d_finite_volumes/domain.py (modified) (6 diffs)

anuga_core/source/anuga/abstract_2d_finite_volumes/domain.py

@@ -573 +573 @@
         # FIXME (Ole): This is under construction. See ticket:192
 
+        from anuga.abstract_2d_finite_volumes.util import\
+             apply_expression_to_dictionary        
+
         if q is None:
             self.quantities_to_be_monitored = None
@@ -582 +585 @@
             q = [q] # Turn argument into a list
 
-        # Check correcness
+        # Check correcness and initialise
+        self.quantities_to_be_monitored = {}
         for quantity_name in q:
             msg = 'Quantity %s is not a valid conserved quantity'\
                   %quantity_name
             
-            assert quantity_name in self.conserved_quantities, msg
+
+            if not quantity_name in self.quantities:
+                # See if this expression is valid
+                apply_expression_to_dictionary(quantity_name, self.quantities)
+
+            # Initialise extrema
+            self.quantities_to_be_monitored[quantity_name] = [None, None]
+            
 
         if polygon is not None:
@@ -598 +609 @@
 
 
-        self.quantities_to_be_monitored = q
+        
         self.monitor_polygon = polygon
         self.monitor_time_interval = time_interval        
@@ -694 +705 @@
 
             x, y = self.get_centroid_coordinates()[k]
-            radius = self.get_radii()[k]
-            area = self.get_areas()[k]      
-            max_speed = self.max_speed[k]           
+            radius = self.get_radii()[k]
+            area = self.get_areas()[k]      
+            max_speed = self.max_speed[k]           
 
             msg += '  Triangle #%d with centroid (%.4f, %.4f), ' %(k, x, y)
-            msg += 'area = %.4f and radius = %.4f ' %(area, radius)
+            msg += 'area = %.4f and radius = %.4f ' %(area, radius)
             msg += 'had the largest computed speed: %.6f m/s ' %(max_speed)
-            if max_speed > 0.0:
+            if max_speed > 0.0:
                 msg += '(timestep=%.6f)\n' %(radius/max_speed)
             else:
-                msg += '(timestep=%.6f)\n' %(0)     
+                msg += '(timestep=%.6f)\n' %(0)     
             
             # Report all quantity values at vertices
@@ -1015 +1026 @@
         """Go through list of boundary objects and update boundary values
         for all conserved quantities on boundary.
-        It is assumed that the ordering of conserved quantities is 
-        consistent between the domain and the boundary object, i.e. 
-        the jth element of vector q must correspond to the jth conserved
-        quantity in domain.
+        It is assumed that the ordering of conserved quantities is 
+        consistent between the domain and the boundary object, i.e. 
+        the jth element of vector q must correspond to the jth conserved
+        quantity in domain.
         """
 
@@ -1024 +1035 @@
         #FIXME: Boundary objects should not include ghost nodes.
         for i, ((vol_id, edge_id), B) in enumerate(self.boundary_objects):
-            if B is None:
-                print 'WARNING: Ignored boundary segment %d (None)'
-            else:
+            if B is None:
+                print 'WARNING: Ignored boundary segment %d (None)'
+            else:
                 q = B.evaluate(vol_id, edge_id)