Changeset 8446 for trunk/anuga_work/development/gareth/tests/wave/run_wave.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/wave/run_wave.py (modified) (2 diffs)

trunk/anuga_work/development/gareth/tests/wave/run_wave.py

@@ -81 +81 @@
 Br = anuga.Reflective_boundary(domain)      # Solid reflective wall
 Bt = anuga.Transmissive_boundary(domain)    # Continue all values on boundary
-Bz = swb2_boundary_conditions.zero_mass_flux_zero_momentum_boundary(domain) # Strong reflections
-Bz1 = swb2_boundary_conditions.zero_mass_flux_zero_t_transmissive_n_momentum_boundary(domain) # Strong reflections
-Bz2 = swb2_boundary_conditions.zero_mass_flux_zero_n_transmissive_t_momentum_boundary(domain) # Strong reflections
+#Bz = swb2_boundary_conditions.zero_mass_flux_zero_momentum_boundary(domain) # Strong reflections
+#Bz1 = swb2_boundary_conditions.zero_mass_flux_zero_t_transmissive_n_momentum_boundary(domain) # Strong reflections
+#Bz2 = swb2_boundary_conditions.zero_mass_flux_zero_n_transmissive_t_momentum_boundary(domain) # Strong reflections
 Bs = anuga.Transmissive_stage_zero_momentum_boundary(domain) # Strong reflections
 Bd = anuga.Dirichlet_boundary([1,0.,0.]) # Constant boundary values
@@ -100 +100 @@
     return -amplitude*sin((1./wave_length)*t*2*pi)
 
+def waveform2(t):
+    return 0.
+
 Bw2 = anuga.shallow_water.boundaries.Transmissive_n_momentum_zero_t_momentum_set_stage_boundary(domain, waveform)
+Bf_in = swb2_boundary_conditions.Characteristic_boundary(domain, waveform)
 #Bw3 = swb2_boundary_conditions.Transmissive_momentum_nudge_stage_boundary(domain, waveform)
+Bf=swb2_boundary_conditions.Characteristic_boundary(domain,waveform2)
                     
 
 
 # Associate boundary tags with boundary objects
-domain.set_boundary({'left': Bw2, 'right': Bt, 'top': Br, 'bottom': Br})
+domain.set_boundary({'left': Bw2, 'right': Bf, 'top': Br, 'bottom': Br})