Changeset 3508 for inundation/examples

Timestamp:

Aug 17, 2006, 6:30:41 PM (19 years ago)

Author:

ole

Message:

Demos for Rosh

Location:

inundation/examples

Files:

• runup.py (added)
• runup_convergence.py (modified) (4 diffs)

inundation/examples/runup_convergence.py ¶

@@ -31 +31 @@
 slope = -0.02     # 1:50 Slope
 highest_point = 3 # Highest elevation (m)
-sea_level = -2    # Mean sea level
-wave_height = 3   # Solitary wave height H
+sea_level = -1    # Mean sea level
+amplitude = 2     # Solitary wave height H
 simulation_name = 'runup_convergence'   
 
@@ -96 +96 @@
 
 def waveform(t): 
-    return wave_height/cosh(t-5)**2
+    return sea_level + amplitude/cosh(t-5)**2
 
 Bw = Time_boundary(domain=domain,     # Time dependent boundary  
@@ -172 +172 @@
 # Print
 
-print 'wave height [m]:                    ', wave_height
+print 'wave height [m]:                    ', amplitude
 runup_height = topography(runup_point, (north+south)/2.)
 print 'run up height [m]:                  ', runup_height 
@@ -184 +184 @@
 
 # Store
-filename = 'maxrunup'+str(wave_height)+'.csv'
+filename = 'maxrunup'+str(amplitude)+'.csv'
 fid = open(filename,'w')    
 s = 'Depth at eastern boundary,Waveheight,Runup distance,Runup height\n'
 fid.write(s)
 
-s = '%.2f,%.2f,%.2f\n' %(wave_height, runup_distance, runup_height)
+s = '%.2f,%.2f,%.2f\n' %(amplitude, runup_distance, runup_height)
 fid.write(s)