Changeset 8446 for trunk/anuga_work/development/gareth/tests/runup_sinusoid/runup_sinusoid.py

Timestamp:

Jun 14, 2012, 12:31:32 AM (13 years ago)

Author:

davies

Message:

bal_dev: Updates

File:

• trunk/anuga_work/development/gareth/tests/runup_sinusoid/runup_sinusoid.py (modified) (2 diffs)

trunk/anuga_work/development/gareth/tests/runup_sinusoid/runup_sinusoid.py

@@ -65 +65 @@
 Bt=anuga.Transmissive_boundary(domain)          # Continue all values of boundary -- not used in this example
 Bd=anuga.Dirichlet_boundary([-0.1*scale_me,0.,0.])       # Constant boundary values -- not used in this example
+def waveform(t):
+    return (0.0*sin(t*2*pi)-0.1)*exp(-t)-0.1
+Bt2=anuga.Transmissive_n_momentum_zero_t_momentum_set_stage_boundary(domain,waveform)
 #Bw=anuga.Time_boundary(domain=domain,
 #       f=lambda t: [(0.0*sin(t*2*pi)-0.1)*exp(-t)-0.1,0.0,0.0]) # Time varying boundary -- get rid of the 0.0 to do a runup.
@@ -82 +85 @@
     yy = domain.quantities['ymomentum'].centroid_values
     dd = domain.quantities['stage'].centroid_values - domain.quantities['elevation'].centroid_values 
-    dd = (dd)*(dd>1.0e-03)+1.0e-06
+    dd = (dd)*(dd>1.0e-06)+1.0e-03
     vv = ( (xx/dd)**2 + (yy/dd)**2 )**0.5
     vv = vv*(dd>1.0e-03)