Changeset 8446 for trunk/anuga_work/development/gareth/tests/runup/runup.py

Timestamp:

Jun 14, 2012, 12:31:32 AM (13 years ago)

Author:

davies

Message:

bal_dev: Updates

File:

• trunk/anuga_work/development/gareth/tests/runup/runup.py (modified) (1 diff)

trunk/anuga_work/development/gareth/tests/runup/runup.py

@@ -67 +67 @@
 #------------------------------
 
-#xwrite=open("xvel.out","wb")
-#ywrite=open("yvel.out","wb")
-## Set print options to be compatible with file writing via the 'print' statement 
-#numpy.set_printoptions(threshold=numpy.nan, linewidth=numpy.nan)
-
-for t in domain.evolve(yieldstep=0.2,finaltime=30.0):
+for t in domain.evolve(yieldstep=0.2,finaltime=60.0):
     print domain.timestepping_statistics()
-
+    uh=domain.quantities['xmomentum'].centroid_values
+    vh=domain.quantities['ymomentum'].centroid_values
+    depth=domain.quantities['stage'].centroid_values - domain.quantities['elevation'].centroid_values
+    depth=depth*(depth>1.0e-06) + 1.0e-06
+    vel=((uh/depth)**2 + (vh/depth)**2)**0.5
+    print 'peak speed is', vel.max()
 
 print 'Finished'