Changeset 3254


Ignore:
Timestamp:
Jun 29, 2006, 10:43:49 AM (18 years ago)
Author:
sexton
Message:

updates for student

File:
1 edited

Legend:

Unmodified
Added
Removed
  • documentation/user_manual/examples/runuptest.py

    r3247 r3254  
    2222#------------------------------------------------------------------------------
    2323
    24 waveheight = 1000.0
     24waveheight = 2000.0
    2525depth_east_edge = -4000.
    2626west = 0.
     
    3636                                        'bottom': [2],
    3737                                        'right': [3]},
    38                          maximum_triangle_area=100000,
     38                         maximum_triangle_area=200000,
    3939                         filename=meshname)
    4040                         #interior_regions=interior_regions)
     
    9191def gauge_line(west,east,north,south):
    9292    from Numeric import arange
    93     gaugex = arange(east,west,-1000.)
     93    gaugex = arange(west,east,1000.)
    9494    gauges = []
    9595    gaugey = []
     
    121121    elevation = []
    122122    depths = []
     123    bed = topography(g[0],g[1])
    123124    for i, t in enumerate(f.get_time()):
    124125        w = f(t, point_id = k)[0]
     
    128129        elevation.append(elev)
    129130        depths.append(elev)
    130         if g[0] < (west+east)/2. and depth <= 0 and count == 0:
    131             count += 1
    132             posx = g[0]
     131        if w-bed <= 0:
     132            count+=1
     133            posx=g[0]
    133134            loc = k
     135        #if g[0] < (west+east)/2. and depth <= 0 and count == 0:
     136        #    count += 1
     137        #    posx = g[0]
     138        #    loc = k
     139       
    134140    maxw.append(max(stage))
    135141    minw.append(min(stage))
     
    152158
    153159figure(1)
    154 plot(gaugex[:loc],maxw[:loc],'g-')
     160plot(gaugex,maxw,'g+',gaugex,topography(gaugex,(north+south)/2.),'r-')
    155161xlabel('x')
    156162ylabel('stage')
Note: See TracChangeset for help on using the changeset viewer.