Changeset 5180

Timestamp:

Mar 28, 2008, 2:13:13 PM (17 years ago)

Author:

steve

Message:

Added an example which has soliton type solutions (tilted.py)

Files:

• anuga_core/source/anuga/shallow_water/test_eq.py (modified) (2 diffs)
• anuga_validation/convergence_study/tilted.py (added)
• anuga_validation/convergence_study/wave.py (modified) (4 diffs)

anuga_core/source/anuga/shallow_water/test_eq.py

@@ -23 +23 @@
         test = -2205.904774487
 
-        print 'depth', depth
+        #print 'depth', depth
         z = Okada_func(length=length, width=width, dip=dip, \
                        x0=x0, y0=y0, strike=strike, depth=depth, \
@@ -38 +38 @@
         assert allclose(z.rake, rake)
 
-        print 'in test', z.test
+        #print 'in test', z.test
         assert allclose(z.test, -2205.904774487)
         
-        print 'hello finished okada'
+        #print 'hello finished okada'
 
     def test_earthquake_tsunami(self):

anuga_validation/convergence_study/wave.py

@@ -42 +42 @@
 # Setup domain
 #------------------------------------------------------------------------------
-dx = 200.
+dx = 100.
 dy = dx
 L = 100000.
@@ -76 +76 @@
 domain.beta_h    = 0.0
 
+interactive_visualisation = False
+
 
 #------------------------------------------------------------------------------
@@ -107 +109 @@
 
 #===============================================================================
-from anuga.visualiser import RealtimeVisualiser
-vis = RealtimeVisualiser(domain)
-vis.render_quantity_height("stage", zScale =10000, dynamic=True)
-vis.colour_height_quantity('stage', (1.0, 0.5, 0.5))
-vis.start()
+if interactive_visualisation:
+    from anuga.visualiser import RealtimeVisualiser
+    vis = RealtimeVisualiser(domain)
+    vis.render_quantity_height("stage", zScale =10000, dynamic=True)
+    vis.colour_height_quantity('stage', (1.0, 0.5, 0.5))
+    vis.start()
 #===============================================================================
 
@@ -121 +124 @@
 for t in domain.evolve(yieldstep = 50.0, finaltime = 60*60.):
     domain.write_time()
-    vis.update()
-    
-vis.evolveFinished()
+    if interactive_visualisation:
+        vis.update()
 
+if interactive_visualisation:
+    vis.evolveFinished()
+