namespace rossby_front_2d_initialization

Overview

Initial conditions for the 2D Rossby front test. More…

namespace rossby_front_2d_initialization {

// global functions

subroutine, public rossby_front_initialize_thickness(
    h h,
    G G,
    GV GV,
    param_file param_file,
    just_read_params just_read_params
    );

subroutine, public rossby_front_initialize_temperature_salinity(
    T T,
    S S,
    h h,
    G G,
    GV GV,
    param_file param_file,
    eqn_of_state eqn_of_state,
    just_read_params just_read_params
    );

subroutine, public rossby_front_initialize_velocity(
    u u,
    v v,
    h h,
    G G,
    GV GV,
    US US,
    param_file param_file,
    just_read_params just_read_params
    );

} // namespace rossby_front_2d_initialization

Detailed Documentation

Initial conditions for the 2D Rossby front test.

Description of the 2d Rossby front initial conditions

Consistent with a linear equation of state, the system has a constant stratification below the mixed layer, stratified in temperature only. Isotherms are flat below the mixed layer and vertical within. Salinity is constant. The mixed layer has a half sine form so that there are no mixed layer or temperature gradients at the side walls.

Below the mixed layer the potential temperature, \(\theta(z)\), is given by

\[\theta(z) = \theta_0 - \Delta \theta \left( z + h_{ML} \right)\]

where \(\theta_0\) and \(\Delta \theta\) are external model parameters.

The depth of the mixed layer, \(H_{ML}\) is

\[h_{ML}(y) = h_{min} + \left( h_{max} - h_{min} \right) \cos{\pi y/L}\]

. The temperature in mixed layer is given by the reference temperature at \(z=h_{ML}\) so that

\[\theta(y,z) = \theta_0 - \Delta \theta \left( z + h_{ML} \right) & \forall & z < h_{ML}(y) T.B.D.\]

Global Functions

subroutine, public rossby_front_initialize_thickness(
    h h,
    G G,
    GV GV,
    param_file param_file,
    just_read_params just_read_params
    )

Initialization of thicknesses in 2D Rossby front test.

Parameters:

g	Grid structure
gv	Vertical grid structure
h	The thickness that is being initialized [H ~> m or kg m-2]
param_file	A structure indicating the open file to parse for model parameter values.
just_read_params	If present and true, this call will only read parameters without changing h.

subroutine, public rossby_front_initialize_temperature_salinity(
    T T,
    S S,
    h h,
    G G,
    GV GV,
    param_file param_file,
    eqn_of_state eqn_of_state,
    just_read_params just_read_params
    )

Initialization of temperature and salinity in the Rossby front test.

Parameters:

g	Grid structure
gv	The ocean’s vertical grid structure.
t	Potential temperature [degC]
s	Salinity [ppt]
h	Thickness [H ~> m or kg m-2]
param_file	Parameter file handle
eqn_of_state	Equation of state structure
just_read_params	If present and true, this call will only read parameters without changing T & S.

subroutine, public rossby_front_initialize_velocity(
    u u,
    v v,
    h h,
    G G,
    GV GV,
    US US,
    param_file param_file,
    just_read_params just_read_params
    )

Initialization of u and v in the Rossby front test.

Parameters:

g	Grid structure
gv	Vertical grid structure
u	i-component of velocity [L T-1 ~> m s-1]
v	j-component of velocity [L T-1 ~> m s-1]
h	Thickness [H ~> m or kg m-2]
us	A dimensional unit scaling type
param_file	A structure indicating the open file to parse for model parameter values.
just_read_params	If present and true, this call will only read parameters without setting u & v.