Changeset 3224

Timestamp:

Jun 24, 2006, 5:11:13 PM (19 years ago)

Author:

sexton

Message:

potential validation test - ref Liu paper

Location:

documentation/user_manual/examples

Files:

• harbour.py (added)
• runuptest.py (modified) (2 diffs)

documentation/user_manual/examples/runuptest.py

@@ -99 +99 @@
 
 f = file_function(gauges)
-minw = 10000.0
-theminw = minw
-minws = []
+maxw = 0.0
+themaxw = maxw
+maxws = []
 for k, g in enumerate(gauges):
     for i, t in enumerate(f.get_time()):
-         if w < minw: minw = w
-    minws.append(minw) 
-    if theminw < minw: 
-       theminw = minw
-       gaugeloc = k
+        if elevation > 0:
+            if w > maxw: maxw = w
+    maxws.append(maxw) 
+    if elevation > 0:
+       if themaxw < maxw: 
+           themaxw = maxw
+           gaugeloc = k
 
 filename = 'maxrunup'+waveheight+'.csv'
@@ -119 +121 @@
 
 from pylab import plot, xlabel, ylabel, title
-plot(gaugex,minws,'g-',gaugestar[0],gaugestar[1],'r*')
+plot(gaugex,maxws,'g-',gaugestar[0],gaugestar[1],'r*')
 xlabel('x')
 ylabel('stage')
-title('Minimum stage for gauge line')
+title('Maximum stage for gauge line')