Changeset 5112

Timestamp:

Mar 4, 2008, 3:58:41 PM (15 years ago)

Author:

ole

Message:

Moved work on Will Power's attenuation problem from
runup demo in the user_manual to the debug area so that demo
matches description in manual and also runs.

Used:

svn merge -r 4883:4882 runup.py

to revert modification.

Files:

• anuga_core/documentation/user_manual/demos/runup.py (modified) (5 diffs)
• anuga_work/debug/will_power_attenuation_problem (added)
• anuga_work/debug/will_power_attenuation_problem/runup.py (added)

anuga_core/documentation/user_manual/demos/runup.py

@@ -10 +10 @@
 #------------------------------------------------------------------------------
 
-from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular_cross
+from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular
 from anuga.shallow_water import Domain
 from anuga.shallow_water import Reflective_boundary
@@ -22 +22 @@
 #------------------------------------------------------------------------------
 
-points, vertices, boundary = rectangular_cross(10, 10,len1=1.0, len2= 1.0 ) # Basic mesh
+points, vertices, boundary = rectangular(10, 10) # Basic mesh
 
 domain = Domain(points, vertices, boundary) # Create domain
@@ -37 +37 @@
 
 def topography(x,y):
-    return -x/2                              # linear bed slope
+    return -x/2                             # linear bed slope
     #return x*(-(2.0-x)*.5)                  # curved bed slope
 
@@ -54 +54 @@
 Bd = Dirichlet_boundary([-0.2,0.,0.]) # Constant boundary values
 Bw = Time_boundary(domain=domain,     # Time dependent boundary  
-                   f=lambda t: [(.1*sin(t*2*pi)-0.3) , 0.0, 0.0])
+                   f=lambda t: [(.1*sin(t*2*pi)-0.3) * exp(-t), 0.0, 0.0])
 
 # Associate boundary tags with boundary objects
 domain.set_boundary({'left': Br, 'right': Bw, 'top': Br, 'bottom': Br})
 
-#===============================================================================
-from anuga.visualiser import RealtimeVisualiser
-vis = RealtimeVisualiser(domain)
-vis.render_quantity_height("elevation", zScale=1.0, dynamic=False)
-vis.render_quantity_height("stage", dynamic=True)
-vis.colour_height_quantity('stage', (lambda q:q['stage'], -0.5, 0.5))
-vis.start()
-#===============================================================================
 
-import time
-t0 = time.time()
 #------------------------------------------------------------------------------
 # Evolve system through time
@@ -76 +66 @@
 for t in domain.evolve(yieldstep = 0.1, finaltime = 10.0):
     domain.write_time()
-    vis.update()
+    
 
-vis.evolveFinished()
-print 'That took %.2f seconds' %(time.time()-t0)