Changeset 1929 for inundation/analytical solutions/Analytical_solution_circular_hydraulic_jump.py

Timestamp:

Oct 14, 2005, 5:22:01 PM (19 years ago)

Author:

ole

Message:

Suggested using the new expression functionality for printing the radial momentum.
I have checked the calculated results are the same.
The old code is still there but has been commented out.

File:

• inundation/analytical solutions/Analytical_solution_circular_hydraulic_jump.py (modified) (3 diffs)

inundation/analytical solutions/Analytical_solution_circular_hydraulic_jump.py

@@ -17 +17 @@
 from math import pi, sqrt
 from mesh_factory import strang_mesh
+
 
 #---------
@@ -134 +135 @@
 #vymom = Visualiser(domain,title='ymomentum',scale_z=10.0)
 
+w = domain.quantities['stage']
 
 #----------
@@ -142 +144 @@
 for t in domain.evolve(yieldstep = .01, finaltime = 10):
     domain.write_time()
+
+    exp = '(xmomentum**2 + ymomentum**2)**0.5'
+    radial_momentum = domain.create_quantity_from_expression(exp)
+
+    print 'outer stage      ', w.get_values(location='centroids',
+                                            indices=[typical_outer[0]])
+    print '      radial mom ', \
+          radial_momentum.get_values(location='centroids',
+                                     indices=[typical_outer[0]])
+
+    print 'inner stage      ', w.get_values(location='centroids',
+                                            indices=[typical_inner[0]])
+    print '      radial mom ', \
+          radial_momentum.get_values(location='centroids',
+                                     indices=[typical_inner[0]])    
+    
+
+    
+
+    
+    
     #vxmom.update_quantity('xmomentum')
     #vymom.update_quantity('ymomentum')
-    print 'outer stage      ',domain.quantities['stage'].get_values(location='centroids',indices=[typical_outer[0]])
-    print '      radial mom ',sqrt(domain.quantities['xmomentum'].get_values(location='centroids',indices=[typical_outer[0]])[0]**2 +
-                                   domain.quantities['ymomentum'].get_values(location='centroids',indices=[typical_outer[0]])[0]**2)
-    print 'inner stage      ',domain.quantities['stage'].get_values(location='centroids',indices=[typical_inner[0]])
-    print '      radial mom ',sqrt(domain.quantities['xmomentum'].get_values(location='centroids',indices=[typical_inner[0]])[0]**2 +
-                                   domain.quantities['ymomentum'].get_values(location='centroids',indices=[typical_inner[0]])[0]**2)
+
+    #Old version
+    #print 'outer stage      ',domain.quantities['stage'].get_values(location='centroids',indices=[typical_outer[0]])
+    #print '      radial mom ',sqrt(domain.quantities['xmomentum'].get_values(location='centroids',indices=[typical_outer[0]])[0]**2 +
+    #                               domain.quantities['ymomentum'].get_values(location='centroids',indices=[typical_outer[0]])[0]**2)
+    #print 'inner stage      ',domain.quantities['stage'].get_values(location='centroids',indices=[typical_inner[0]])
+    #print '      radial mom ',sqrt(domain.quantities['xmomentum'].get_values(location='centroids',indices=[typical_inner[0]])[0]**2 +
+    #                               domain.quantities['ymomentum'].get_values(location='centroids',indices=[typical_inner[0]])[0]**2)