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

Timestamp:

Apr 14, 2007, 4:05:58 PM (16 years ago)

Author:

ole

Message:

Prepared tests for making new slope limiters default. Still a few unprepared tests remaining and
one example in test_data_manager where small timesteps are generated.

File:

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

anuga_core/source/anuga/abstract_2d_finite_volumes/domain.py

@@ -630 +630 @@
 
             x, y = self.get_centroid_coordinates()[k]
-
-            msg += '  Triangle #%d with centroid (%.4f, %.4f) ' %(k, x, y)
-            msg += 'had the largest computed speed: %.4f m/s\n' %(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 += 'had the largest computed speed: %.6f m/s ' %(max_speed)
+            if max_speed > 0.0:
+                msg += '(timestep=%.6f)\n' %(radius/max_speed)
+            else:
+                msg += '(timestep=%.6f)\n' %(0)     
             
             # Report all quantity values at vertices