Changeset 3486

Timestamp:

Aug 11, 2006, 8:30:52 PM (18 years ago)

Author:

ole

Message:

Attempt to make runup example more interesting and also using the symmetric mesh generated by rectangular_cross

File:

• documentation/user_manual/examples/runup.py (modified) (4 diffs)

documentation/user_manual/examples/runup.py

@@ -10 +10 @@
 #------------------------------------------------------------------------------
 
-from pyvolution.mesh_factory import rectangular
+from pyvolution.mesh_factory import rectangular_cross
 from pyvolution.shallow_water import Domain
 from pyvolution.shallow_water import Reflective_boundary
@@ -22 +22 @@
 #------------------------------------------------------------------------------
 
-points, vertices, boundary = rectangular(10, 10) # Basic mesh
+points, vertices, boundary = rectangular_cross(10, 10) # Basic mesh
 
 domain = Domain(points, vertices, boundary) # Create domain
-domain.set_name('bedslope')                 # Output to bedslope.sww
+domain.set_name('runup')                    # Output to bedslope.sww
 domain.set_datadir('.')                     # Use current directory for output
-domain.set_quantities_to_be_stored(['stage',# Store all conserved quantities
-                                    'xmomentum',
-                                    'ymomentum'])   
+#domain.set_quantities_to_be_stored(['stage',# Store all conserved quantities
+#                                    'xmomentum',
+#                                    'ymomentum'])   
 
 
@@ -52 +52 @@
 Br = Reflective_boundary(domain)      # Solid reflective wall
 Bt = Transmissive_boundary(domain)    # Continue all values on boundary 
-Bd = Dirichlet_boundary([-0.5,0.,0.]) # Constant boundary values
+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), 0.0, 0.0])
 
 # Associate boundary tags with boundary objects
@@ -64 +64 @@
 #------------------------------------------------------------------------------
 
-for t in domain.evolve(yieldstep = 0.1, finaltime = 4.0):
+for t in domain.evolve(yieldstep = 0.1, finaltime = 10.0):
     domain.write_time()