Changeset 7995

Timestamp:

Sep 2, 2010, 4:01:03 PM (14 years ago)

Author:

steve

Message:

Added stats

Location:

trunk/anuga_core/source/anuga/structures

Files:

• boyd_box_operator.py (modified) (2 diffs)
• structure_operator.py (modified) (3 diffs)
• testing_wide_bridge.py (modified) (1 diff)

trunk/anuga_core/source/anuga/structures/boyd_box_operator.py

@@ -48 +48 @@
         self.culvert_height = self.get_culvert_height()
 
+
         self.sum_loss = 0.0
         self.max_velocity = 10.0
         self.log_filename = None
-        
+
         self.inlets = self.get_inlets()
+
+
+        # Stats
+        
+        self.discharge = 0.0
+        self.velocity = 0.0
         
         
@@ -111 +118 @@
         #print Q, Q*timestep, barrel_speed, outlet_depth, Qstar, factor, timestep_star
         #print '  ', loss, gain
+
+        # Stats
+        self.discharge  = outflow_extra_height*self.outflow.get_area()/timestep
+        self.velocity = self.discharge/outlet_depth/self.width
 
         new_outflow_height = self.outflow.get_average_height() + outflow_extra_height

trunk/anuga_core/source/anuga/structures/structure_operator.py

@@ -6 +6 @@
 
 class Structure_operator:
-    """Culvert flow - transfer water from one rectangular box to another.
+    """Structure Operator - transfer water from one rectangular box to another.
     Sets up the geometry of problem
     
@@ -43 +43 @@
         self.enquiry_gap = enquiry_gap
         self.verbose = verbose
+
+        self.discharge = 0.0
+        self.velocity = 0.0
         
         self.__create_exchange_polygons()
@@ -149 +152 @@
         print '====================================='
 
-                        
+
+    def structure_statistics(self):
+
+        message = '---------------------------\n'
+        message += 'Structure report:\n'
+        message += '--------------------------\n'
+        message += 'Discharge [m^3/s]: %.2f\n' % self.discharge
+        message += 'Velocity  [m/s]: %.2f\n' % self.velocity
+#        message += 'Total boundary outflow [m^3/s]: %.2f\n' % total_boundary_outflow
+#        message += 'Net boundary flow by tags [m^3/s]\n'
+#        for tag in boundary_flows:
+#            message += '    %s [m^3/s]: %.2f\n' % (tag, boundary_flows[tag])
+#
+#        message += 'Total net boundary flow [m^3/s]: %.2f\n' % \
+#                    (total_boundary_inflow + total_boundary_outflow)
+#        message += 'Total volume in domain [m^3]: %.2f\n' % \
+#                    self.compute_total_volume()
+#
+#        # The go through explicit forcing update and record the rate of change
+#        # for stage and
+#        # record into forcing_inflow and forcing_outflow. Finally compute
+#        # integral of depth to obtain total volume of domain.
+#
+        # FIXME(Ole): This part is not yet done.
+
+        return message
+
     def get_inlets(self):
         

trunk/anuga_core/source/anuga/structures/testing_wide_bridge.py

@@ -256 +256 @@
     print domain.timestepping_statistics()
     print domain.volumetric_balance_statistics()
+    for culvert in culverts:
+        print culvert.structure_statistics()