Changeset 3511 for inundation/examples

Timestamp:

Aug 21, 2006, 1:17:30 PM (19 years ago)

Author:

ole

Message:

Grid on and snapshots as per request from Rosh.

File:

• inundation/examples/runup_convergence.py (modified) (9 diffs)

inundation/examples/runup_convergence.py

@@ -22 +22 @@
 from pyvolution.shallow_water import Transmissive_Momentum_Set_Stage_boundary
 from pyvolution.util import file_function
-from pylab import plot, xlabel, ylabel, title, ion, close, savefig, figure, axis
-ion()
+from pylab import plot, xlabel, ylabel, title, ion, close, savefig, figure, axis, legend, grid, hold
+
 
 
@@ -96 +96 @@
 
 def waveform(t): 
-    return sea_level + amplitude/cosh(t-5)**2
+    return sea_level + amplitude/cosh(t-10)**2
 
 Bw = Time_boundary(domain=domain,     # Time dependent boundary  
@@ -114 +114 @@
 
 stagestep = []
-for t in domain.evolve(yieldstep = 1, finaltime = 100):
+for t in domain.evolve(yieldstep = 1, finaltime = 150):
     domain.write_time()
 
@@ -126 +126 @@
 def gauge_line(west,east,north,south):
     from Numeric import arange
-    x_vector = arange(west, (east-west)/2, 10) # Gauges every 1 meter from west to midpt
+    x_vector = arange(west, east, dx/4) # Gauges every dx/2 meter from west to east
     y = (north+south)/2.
 
@@ -142 +142 @@
                   interpolation_points = gauges,
                   verbose = True,
-                  use_cache = True)
+                  use_cache = False)
 
 
@@ -171 +171 @@
 
 # Print
-
 print 'wave height [m]:                    ', amplitude
 runup_height = topography(runup_point, (north+south)/2.)
@@ -182 +181 @@
 
 
+# Take snapshots and plot
+ion()
+figure(1)
+plot(x_vector, topography(x_vector,(north+south)/2.), 'r-')
+xlabel('x')
+ylabel('Elevation')
+#legend(('Max stage', 'Min stage', 'Elevation'), shadow=True, loc='upper right')
+title('Stage snapshots (t=0, 10, ...) for gauge line')
+grid()
+hold(True)
+
+for i, t in enumerate(f.get_time()):
+    if i % 10 == 0:
+        # Take only some timesteps to avoid clutter
+        stages = []    
+        for k, g in enumerate(gauges):
+            w = f(t, point_id = k)[0]        
+            stages.append(w)
+
+        plot(x_vector, stages, 'b-')
+         
+savefig('snapshots')
+        
+
+
 
 # Store
 filename = 'maxrunup'+str(amplitude)+'.csv'
 fid = open(filename,'w')    
-s = 'Depth at eastern boundary,Waveheight,Runup distance,Runup height\n'
+s = 'Waveheight,Runup distance,Runup height\n'
 fid.write(s)
 
@@ -194 +218 @@
 fid.close()
 
-# Plot
-figure(1)
+
+
+
+
+# Plot max runup etc
+ion()
+hold(False)
+figure(2)
 plot(x_vector, max_stage, 'g+',
      x_vector, min_stage, 'b+',     
@@ -201 +231 @@
 xlabel('x')
 ylabel('stage')
+legend(('Max stage', 'Min stage', 'Elevation'), shadow=True, loc='upper right')
 title('Maximum stage for gauge line')
+grid()
 #axis([33000, 47000, -1000, 3000])
 savefig('max_stage')