Changeset 8212 for trunk/anuga_validation/validation_tests/dam_break.py

Timestamp:

Sep 13, 2011, 1:44:23 PM (14 years ago)

Author:

steve

Message:

Committing Sudi's changes

File:

• trunk/anuga_validation/validation_tests/dam_break.py (modified) (7 diffs)

trunk/anuga_validation/validation_tests/dam_break.py

@@ -8 +8 @@
 # Import necessary modules
 #------------------------------------------------------------------------------
-
 import sys
 import anuga
-
 from math import cos
 from numpy import zeros, float
@@ -27 +25 @@
 output_file = 'dam_break'
 
-anuga.copy_code_files(output_dir,__file__)
+#anuga.copy_code_files(output_dir,__file__)
 #start_screen_catcher(output_dir+'_')
 
@@ -72 +70 @@
 domain.set_quantity('friction', 0.0)
 
-h0 = 10.0
+h0 = 10000.0
 h1 = 1.0
 
@@ -83 +81 @@
             z[i] = h1
     return z
+domain.set_quantity('stage', height)
 
-
-
-
-
-
-domain.set_quantity('stage', height)
 #-----------------------------------------------------------------------------
 # Setup boundary conditions
@@ -98 +91 @@
 Bd = anuga.Dirichlet_boundary([1,0.,0.]) # Constant boundary values
 
-
 # Associate boundary tags with boundary objects
 domain.set_boundary({'left': Bt, 'right': Bt, 'top': Br, 'bottom': Br})
@@ -104 +96 @@
 
 #===============================================================================
-from anuga.visualiser import RealtimeVisualiser
-vis = RealtimeVisualiser(domain)
-vis.render_quantity_height("stage", zScale =h0*500, dynamic=True)
-vis.colour_height_quantity('stage', (0.0, 0.5, 1.0))
-vis.start()
+##from anuga.visualiser import RealtimeVisualiser
+##vis = RealtimeVisualiser(domain)
+##vis.render_quantity_height("stage", zScale =h0*500, dynamic=True)
+##vis.colour_height_quantity('stage', (0.0, 0.5, 1.0))
+##vis.start()
 #===============================================================================
 
@@ -115 +107 @@
 # Evolve system through time
 #------------------------------------------------------------------------------
-
-for t in domain.evolve(yieldstep = 100.0, finaltime = 60*60.):
+for t in domain.evolve(yieldstep = 1.0, finaltime = 10*10.):
     #print domain.timestepping_statistics(track_speeds=True)
     print domain.timestepping_statistics()
-    vis.update()
+    #vis.update()
 
 
 #test against know data
     
-vis.evolveFinished()
+#vis.evolveFinished()