Changeset 3248

Timestamp:

Jun 27, 2006, 6:16:24 PM (18 years ago)

Author:

sexton

Message:

small change

File:

• documentation/user_manual/examples/harbour.py (modified) (4 diffs)

documentation/user_manual/examples/harbour.py

@@ -25 +25 @@
 harbour_length = 2000.0 
 shelf_length = 10000.0
+bottom_length = 10000.0
 harbour_width = 250.0
 shelf_depth = -100. 
@@ -33 +34 @@
 #------------------------------------------------------------------------------
 
-centre = 340000.
-west = centre-harbour_length_shelf_length-10000.
-east = centre+10000.
-south = 6265000.
-north = 6250000.
+west = 10000.+harbour_length+shelf_length+bottom_length
+east = 0.
+south = 0.
+north = 5000.
 polygon = [[east,north],[west,north],[west,south],[east,south]]
 meshname = 'harbour_test.msh'
@@ -113 +113 @@
 def gauge_line(west,east,north,south):
     from Numeric import arange
-    gaugex = arange(west,east,1000.)
+    gaugex = arange(east,west,-1000.)
     gauges = []
+    gaugey = []
     for x in gaugex:
         y = (north+south)/2.
+        gaugey.append(y)
         gauges.append([x,y])
-    return gauges
+    return gauges, gaugex, gaugey
 
-gauges = gauge_line(west,east,north,south)
+gauges, gaugex, gaugey = gauge_line(west,east,north,south)
 
 from pyvolution.util import file_function
@@ -130 +132 @@
                   use_cache = True)
 
-maxw = 0.0
-themaxw = maxw
-maxws = []
+from pylab import plot, xlabel, ylabel, title, ion, close, savefig, figure
+ion()
+
+maxw = []
+minw = []
+maxd = []
+count = 0
 for k, g in enumerate(gauges):
+    stage = []
+    elevation = []
+    depths = []
     for i, t in enumerate(f.get_time()):
         w = f(t, point_id = k)[0]
-        elevation = f(t, point_id = k)[1]
-        if w > maxw: maxw = w
-    maxws.append(maxw) 
-    if themaxw < maxw:
-        themaxw = maxw
-        gaugeloc = k
+        elev = f(t, point_id = k)[1]
+        depth = w-elev
+        stage.append(w)
+        elevation.append(elev)
+        depths.append(elev)
+    maxw.append(max(stage))
+    minw.append(min(stage))
+    maxd.append(max(depths))
 
-from pylab import plot, xlabel, ylabel, title, ion, close, savefig
-ion()
-plot(gaugex,maxws,'g-',gaugestar[0],gaugestar[1],'r*')
+filename = 'point'+str(waveheight)+'.csv'
+fid = open(filename,'w')    
+s = 'Waveheight,\n'
+fid.write(s)
+
+figure(1)
+plot(gaugex[:loc],maxw[:loc],'g-')
 xlabel('x')
 ylabel('stage')
 title('Maximum stage for gauge line')
-savefig('max_stage_harbour')
+savefig('max_stage')
+
 close('all')