Changeset 4347 for anuga_work/production/pt_hedland_2006

Timestamp:

Apr 2, 2007, 6:23:29 PM (18 years ago)

Author:

sexton

Message:

update export_results script with updates calculation for speed

File:

• anuga_work/production/pt_hedland_2006/export_results.py (modified) (1 diff)

anuga_work/production/pt_hedland_2006/export_results.py

@@ -44 +44 @@
 if which_var == 3:  # Speed
     outname = name + '_speed'
-    #quantityname = '((xmomentum/(stage-elevation))**2 + (ymomentum/(stage-elevation))**2)**0.5'  #Speed
-    #quantityname = 'xmomentum/(stage-elevation)'  #Speed
-    quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-30)'  #Speed
+    quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
 
 if which_var == 4:  # Elevation