## Bessel monopole encounters a north-south ridge

*Details on the computations*

**Initialisation**
Monopole: Bessel monopole of radius 0.5 and strength 4
centred at (+3,-3)
Tracers: single tracer point at the centre of the vortex
Computation: domain: X = (-10, ..., 10), Y = (-10, ..., 10)
no. of grid points from wall to wall: (257, 257)
boundary conditions: stress-free along all walls
standard Arakawa scheme used
rotation effects: beta = 0.3, gamma = 0
and: f0 = 0, +3 or -3
bottom topography: smooth cosine-shaped ridge
along the y-axis; fluid depth:
1.d0 - ampl*dcos((xx-xc)*pi/wd) - ampl
for |xx| <= wd; elsewhere: 1.d0
with: ampl = 0.2 ==> max. height = 0.4
wd = 1.0
solver for streamfunction: NAG
viscosity: nu = 10^-3 (i.e. Re=1000 at T=0)
time step: dt = 0.05
final time: T = 50

Bessel monopole encounters a north-south ridge:

<=== with `f0`=0:
*equatorial case*

<=== with `f0`=+3:
*northern hemisphere*

<=== with `f0`=-3:
*southern hemisphere*

**
Jos van Geffen** --
Home |
Site Map |
Contact Me

###### created: 10 April 1998

last modified: 26 May 2001