Coriolis

The Coriolis option determines whether the fluid experiences the effect of the Coriolis force, or rotation. Currently three options are available: no rotation, -plane, and -plane.

Coriolis vs. rotation

If you are wondering why this option is called "Coriolis" it is because rotational effects could include the Coriolis and centripetal forces, both of which arise in non-inertial reference frames. But here the model only considers the Coriolis force.

No rotation

By default there is no rotation. This can be made explicit by passing coriolis = nothing to a model constructor.

Traditional -plane

To set up an -plane with, for example, Coriolis parameter  

julia> coriolis = FPlane(f=1e-4)
FPlane{Float64}(f=0.0001)

An -plane can also be specified at some latitude on a spherical planet with a planetary rotation rate. For example, to specify an -plane at a latitude of    on Earth which has a rotation rate of   

julia> coriolis = FPlane(rotation_rate=7.292115e-5, latitude=45)
FPlane{Float64}(f=0.000103126)

in which case the value of is given by .

Coriolis term for constant rotation in a Cartesian coordinate system

One can use ConstantCartesianCoriolis to set up a Coriolis acceleration term where the Coriolis parameter is constant and the rotation axis is arbitrary. For example, with    ,

julia> coriolis = ConstantCartesianCoriolis(fx=0, fy=2e-4, fz=1e-4)
ConstantCartesianCoriolis{Float64}: fx = 0.00e+00, fy = 2.00e-04, fz = 1.00e-04

Or alternatively, the same result can be achieved by specifying the magnitude of the Coriolis frequency f and the rotation_axis. So another way to get a Coriolis acceleration with the same values is:

julia> rotation_axis = (0, 2e-4, 1e-4)./√(2e-4^2 + 1e-4^2) # rotation_axis has to be a unit vector
(0.0, 0.8944271909999159, 0.4472135954999579)

julia> coriolis = ConstantCartesianCoriolis(f=√(2e-4^2+1e-4^2), rotation_axis=rotation_axis)
ConstantCartesianCoriolis{Float64}: fx = 0.00e+00, fy = 2.00e-04, fz = 1.00e-04

An -plane with non-traditional Coriolis terms can also be specified at some latitude on a spherical planet with a planetary rotation rate. For example, to specify an -plane at a latitude of    on Earth which has a rotation rate of   

julia> coriolis = ConstantCartesianCoriolis(rotation_rate=7.292115e-5, latitude=45)
ConstantCartesianCoriolis{Float64}: fx = 0.00e+00, fy = 1.03e-04, fz = 1.03e-04

in which case   and  .

Traditional -plane {#Traditional-\beta-plane}

To set up a -plane the background rotation rate and the parameter must be specified. For example, a -plane with   and    can be set up with

julia> coriolis = BetaPlane(f₀=1e-4, β=1.5e-11)
BetaPlane{Float64}(f₀=0.0001, β=1.5e-11)

Alternatively, a -plane can also be set up at some latitude on a spherical planet with a planetary rotation rate and planetary radius. For example, to specify a -plane at a latitude of   on Earth which has a rotation rate of    and a radius of  

julia> coriolis = BetaPlane(rotation_rate=7.292115e-5, latitude=-10, radius=6371e3)
BetaPlane{Float64}(f₀=-2.53252e-5, β=2.25438e-11)

in which case   and  .

Non-traditional -plane {#Non-traditional-\beta-plane}

A non-traditional -plane requires either 5 parameters (by default Earth's radius and rotation rate are used):

julia> NonTraditionalBetaPlane(fz=1e-4, fy=2e-4, β=4e-11, γ=-8e-11)
NonTraditionalBetaPlane{Float64}(fz = 1.00e-04, fy = 2.00e-04, β = 4.00e-11, γ = -8.00e-11, R = 6.37e+06)

or the rotation rate, radius, and latitude:

julia> NonTraditionalBetaPlane(rotation_rate=5.31e-5, radius=252.1e3, latitude=10)
NonTraditionalBetaPlane{Float64}(fz = 1.84e-05, fy = 1.05e-04, β = 4.15e-10, γ = -1.46e-10, R = 2.52e+05)