Changeset 4722

Timestamp:

Sep 11, 2007, 9:39:35 AM (17 years ago)

Author:

steve

Message:

 

File:

• anuga_validation/convergence_study/convergence_structured.py (modified) (5 diffs)

anuga_validation/convergence_study/convergence_structured.py

@@ -35 +35 @@
 #points, vertices, boundary = rectangular_cross(1000, 100, 20, 3) # Basic mesh
 domain = Domain(points, vertices, boundary) 
-
+domain.set_timestepping_method('rk2')
+domain.set_default_order(2)
 ## # unstructured mesh
 ## poly_domain = [[0,-W],[0,W],[L,W],[L,-W]]
@@ -49 +50 @@
 ## # Create domain
 ## domain = Domain(meshname, use_cache=True, verbose = True)
-domain.set_name('myexample3')                
+domain.set_name('myexample6')                
 domain.set_datadir('.')                     # Use current directory for output
 
@@ -68 +69 @@
 Bd = Dirichlet_boundary([1,0.,0.]) # Constant boundary values
 amplitude = 1
+#Bw = Transmissive_Momentum_Set_Stage_boundary(domain=domain,
 Bw = Time_boundary(domain=domain,     # Time dependent boundary  
 ## Sine wave
@@ -81 +83 @@
 domain.set_boundary({'left': Bw, 'right': Bt, 'top': Br, 'bottom': Br})
 
-vis = RealtimeVisualiser(domain)
-vis.render_quantity_height("elevation", dynamic=False)
-vis.render_quantity_height("stage", dynamic=True)
-vis.colour_height_quantity('stage', (0.0, 0.0, 0.8))
-vis.start()
+#from anuga.visualiser import RealtimeVisualiser
+#vis = RealtimeVisualiser(domain)
+#vis.render_quantity_height("elevation", dynamic=False)
+#vis.render_quantity_height("stage", dynamic=True)
+#vis.colour_height_quantity('stage', (0.0, 0.0, 0.8))
+#vis.start()
 
 #------------------------------------------------------------------------------
@@ -93 +96 @@
 for t in domain.evolve(yieldstep = 20.0, finaltime = 40*60.):
     domain.write_time()
-    vis.update()
+    #vis.update()
     
-vis.evolveFinished()
-vis.join()
+#vis.evolveFinished()
+#vis.join()