GfsSourceCoriolis

From Gerris

As its name suggests, GfsSourceCoriolis adds a Coriolis source term to the (horizontal) components of the velocity or momentum. The object is further generalised and can be used to optionally add a linear friction term (which has a formally similar form to the Coriolis terms).

The syntax in parameter files is

[ GfsSourceVelocity ] [ GfsFunction ] ( [ GfsFunction ] )

where the second GfsFunction is optional. The first GfsFunction is the value of the Coriolis coefficient, the second GfsFunction is the value of the linear friction coefficient.

For example

SourceCoriolis 1e-4 3e-3

is formally equivalent to

Source U (-1e-4*V - 3e-3*U)
Source V ( 1e-4*U - 3e-3*V)

but uses a robust implicit discretisation.

Examples

• Lunar tides in Cook Strait, New Zealand

    SourceCoriolis -1e-4

• Geostrophic adjustment

  SourceCoriolis {} 1

• Geostrophic adjustment on a beta-plane

  SourceCoriolis {} (1. + 0.156246961595*(y + 0.5))

• Coastally-trapped waves

  SourceCoriolis {} 1.

• Coastally-trapped waves with adaptive refinement

  SourceCoriolis {} 1.

• Oscillations in a parabolic container

    SourceCoriolis 0 tau

• Circular dam break on a rotating sphere

    SourceCoriolis 10.*sin(y*M_PI/180.)