Changeset 1043 for inundation/ga/storm_surge/pyvolution/domain.py

Timestamp:

Mar 8, 2005, 5:54:25 PM (20 years ago)

Author:

ole

Message:

Added diagnostics about boundary forcing

File:

• inundation/ga/storm_surge/pyvolution/domain.py (modified) (2 diffs)

inundation/ga/storm_surge/pyvolution/domain.py

@@ -253 +253 @@
 
         self.boundary_objects = []
+        self.boundary_map = boundary_map  #Store for use with eg. boundary_stats.
 
         #FIXME: Try to remove the sorting and fix test_mesh.py
@@ -316 +317 @@
                     self.number_of_first_order_steps)
 
+    def boundary_stats(self, quantities = None, tag = None):
+        """Output statistics about boundary forcing
+
+        
+        """
+        
+        if quantities is None:
+            quantities = self.conserved_quantities
+
+
+        print 'Boundary values at time %.4f:' %self.time
+        for name in quantities:
+            q = self.quantities[name]
+
+            if tag is None:
+                #Take entire boundary
+                print '    Quantity %s: min = %12.8f, max = %12.8f'\
+                      %(name, min(q.boundary_values), max(q.boundary_values))
+            else:
+                #Take only boundary associated with tag
+                maxval = minval = None
+                for i, ((vol_id, edge_id), B) in enumerate(self.boundary_objects):
+                    if self.boundary[(vol_id, edge_id)] == tag:
+                        v = q.boundary_values[i]
+                        if minval is None or v < minval: minval = v
+                        if maxval is None or v > maxval: maxval = v                        
+
+                if minval is None or maxval is None:
+                    print 'Sorry no information about tag %s' %tag
+                else:    
+                    print '    Quantity %s, tag %s: min = %12.8f, max = %12.8f'\
+                          %(name, tag, minval, maxval)
+                        
+
+
+                
+
+        
+        
+           
+        
 
     def get_name(self):