Changeset 3224


Ignore:
Timestamp:
Jun 24, 2006, 5:11:13 PM (18 years ago)
Author:
sexton
Message:

potential validation test - ref Liu paper

Location:
documentation/user_manual/examples
Files:
1 added
1 edited

Legend:

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

    r3223 r3224  
    9999
    100100f = file_function(gauges)
    101 minw = 10000.0
    102 theminw = minw
    103 minws = []
     101maxw = 0.0
     102themaxw = maxw
     103maxws = []
    104104for k, g in enumerate(gauges):
    105105    for i, t in enumerate(f.get_time()):
    106          if w < minw: minw = w
    107     minws.append(minw)
    108     if theminw < minw:
    109        theminw = minw
    110        gaugeloc = k
     106        if elevation > 0:
     107            if w > maxw: maxw = w
     108    maxws.append(maxw)
     109    if elevation > 0:
     110       if themaxw < maxw:
     111           themaxw = maxw
     112           gaugeloc = k
    111113
    112114filename = 'maxrunup'+waveheight+'.csv'
     
    119121
    120122from pylab import plot, xlabel, ylabel, title
    121 plot(gaugex,minws,'g-',gaugestar[0],gaugestar[1],'r*')
     123plot(gaugex,maxws,'g-',gaugestar[0],gaugestar[1],'r*')
    122124xlabel('x')
    123125ylabel('stage')
    124 title('Minimum stage for gauge line')
     126title('Maximum stage for gauge line')
Note: See TracChangeset for help on using the changeset viewer.