Ignore:
Timestamp:
May 26, 2009, 11:33:53 AM (15 years ago)
Author:
rwilson
Message:

More changes - make text agree with changed Cairns code.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/documentation/user_manual/demos/runup.py

    r7064 r7086  
    2929# Setup initial conditions
    3030#------------------------------------------------------------------------------
    31 def topography(x,y):
     31def topography(x, y):
    3232    return -x/2                              # linear bed slope
    3333    #return x*(-(2.0-x)*.5)                  # curved bed slope
     
    3535domain.set_quantity('elevation', topography) # Use function for elevation
    3636domain.set_quantity('friction', 0.1)         # Constant friction
    37 domain.set_quantity('stage', -.4)            # Constant negative initial stage
     37domain.set_quantity('stage', -0.4)           # Constant negative initial stage
    3838
    3939#------------------------------------------------------------------------------
     
    4444Bd = Dirichlet_boundary([-0.2,0.,0.]) # Constant boundary values
    4545Bw = Time_boundary(domain=domain,     # Time dependent boundary 
    46                    f=lambda t: [(.1*sin(t*2*pi)-0.3) * exp(-t), 0.0, 0.0])
     46                   f=lambda t: [(0.1*sin(t*2*pi)-0.3)*exp(-t), 0.0, 0.0])
    4747
    4848# Associate boundary tags with boundary objects
     
    5252# Evolve system through time
    5353#------------------------------------------------------------------------------
    54 for t in domain.evolve(yieldstep = 0.1, finaltime = 10.0):
     54for t in domain.evolve(yieldstep=0.1, finaltime=10.0):
    5555    print domain.timestepping_statistics()
Note: See TracChangeset for help on using the changeset viewer.