Changeset 3254

Timestamp:

Jun 29, 2006, 10:43:49 AM (19 years ago)

Author:

sexton

Message:

updates for student

File:

• documentation/user_manual/examples/runuptest.py (modified) (6 diffs)

documentation/user_manual/examples/runuptest.py

@@ -22 +22 @@
 #------------------------------------------------------------------------------
 
-waveheight = 1000.0
+waveheight = 2000.0
 depth_east_edge = -4000.
 west = 0.
@@ -36 +36 @@
                                         'bottom': [2],
                                         'right': [3]},
-                         maximum_triangle_area=100000,
+                         maximum_triangle_area=200000,
                          filename=meshname)
                          #interior_regions=interior_regions)
@@ -91 +91 @@
 def gauge_line(west,east,north,south):
     from Numeric import arange
-    gaugex = arange(east,west,-1000.)
+    gaugex = arange(west,east,1000.)
     gauges = []
     gaugey = []
@@ -121 +121 @@
     elevation = []
     depths = []
+    bed = topography(g[0],g[1])
     for i, t in enumerate(f.get_time()):
         w = f(t, point_id = k)[0]
@@ -128 +129 @@
         elevation.append(elev)
         depths.append(elev)
-        if g[0] < (west+east)/2. and depth <= 0 and count == 0:
-            count += 1
-            posx = g[0]
+        if w-bed <= 0:
+            count+=1
+            posx=g[0]
             loc = k
+        #if g[0] < (west+east)/2. and depth <= 0 and count == 0:
+        #    count += 1
+        #    posx = g[0]
+        #    loc = k
+        
     maxw.append(max(stage))
     minw.append(min(stage))
@@ -152 +158 @@
 
 figure(1)
-plot(gaugex[:loc],maxw[:loc],'g-')
+plot(gaugex,maxw,'g+',gaugex,topography(gaugex,(north+south)/2.),'r-')
 xlabel('x')
 ylabel('stage')