Changeset 5643

Timestamp:

Aug 13, 2008, 8:14:27 AM (15 years ago)

Author:

kristy

Message:

Updated script so that get_timeseries.py can still be acquired even though some of the gauges could be out of the domain.

File:

• anuga_core/source/anuga/abstract_2d_finite_volumes/util.py (modified) (1 diff)

anuga_core/source/anuga/abstract_2d_finite_volumes/util.py

@@ -2262 +2262 @@
                             vel = 0.0
                         else:
-                            momentum = sqrt(point_quantities[2]**2 +\
-                                            point_quantities[3]**2)
-#                            vel = momentum/depth              
-                            vel = momentum/(point_quantities[0] - point_quantities[1])
-#                            vel = momentum/(depth + 1.e-6/depth)              
-                        
+                            if point_quantities[2] < 1.0e6:
+                                momentum = sqrt(point_quantities[2]**2 +\
+                                                point_quantities[3]**2)
+    #                            vel = momentum/depth              
+                                vel = momentum/(point_quantities[0] - point_quantities[1])
+    #                            vel = momentum/(depth + 1.e-6/depth)
+                            else:
+                                momentum = 0
+                                vel = 0
+                            
                         points_list.append(vel)