Changeset 1839 for inundation/analytical solutions

Timestamp:

Sep 16, 2005, 12:11:41 PM (19 years ago)

Author:

steve

Message:

Investigating Circular Hydraulic jump.

File:

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

inundation/analytical solutions/Analytical_solution_circular_hydraulic_jump.py

@@ -113 +113 @@
 # Order of accuracy
 domain.default_order = 2
+domain.CFL = 1.0
+#domain.beta_w = 0.5
+#domain.beta_h = 0.2
 domain.smooth = True
+
+
 
 
@@ -119 +124 @@
     #domain.visualiser.coloring['stage'] = True
 domain.visualiser.scale_z['stage'] = 2.0
-#domain.visualiser.scale_z['elevation'] = 0.05
+domain.visualiser.scale_z['elevation'] = 0.05
+
+
+from realtime_visualisation_new import Visualiser
+vxmom = Visualiser(domain,title='xmomentum',scale_z=10.0)
+vymom = Visualiser(domain,title='ymomentum',scale_z=10.0)
+
+
 #----------
 # Evolution
 import time
+
+print 'PROBLEM WITH INSTABILITY AFTER t'
 t0 = time.time()
 for t in domain.evolve(yieldstep = .01, finaltime = 10):
     domain.write_time()
+    vxmom.update_quantity('xmomentum')
+    vymom.update_quantity('ymomentum')
 
 print 'That took %.2f seconds' %(time.time()-t0)