namespace mom_pressureforce_blk_afv

Overview

Analytically integrated finite volume pressure gradient. More…

namespace mom_pressureforce_blk_afv {

// global functions

subroutine, public pressureforce_blk_afv(
    h h,
    tv tv,
    PFu PFu,
    PFv PFv,
    G G,
    GV GV,
    US US,
    CS CS,
    ALE_CSp ALE_CSp,
    p_atm p_atm,
    pbce pbce,
    eta eta
    );

subroutine, public pressureforce_blk_afv_nonbouss(
    h h,
    tv tv,
    PFu PFu,
    PFv PFv,
    G G,
    GV GV,
    US US,
    CS CS,
    p_atm p_atm,
    pbce pbce,
    eta eta
    );

subroutine, public pressureforce_blk_afv_bouss(
    h h,
    tv tv,
    PFu PFu,
    PFv PFv,
    G G,
    GV GV,
    US US,
    CS CS,
    ALE_CSp ALE_CSp,
    p_atm p_atm,
    pbce pbce,
    eta eta
    );

subroutine, public pressureforce_blk_afv_init(
    Time Time,
    G G,
    GV GV,
    US US,
    param_file param_file,
    diag diag,
    CS CS,
    tides_CSp tides_CSp
    );

subroutine, public pressureforce_blk_afv_end(CS CS);

} // namespace mom_pressureforce_blk_afv

Detailed Documentation

Analytically integrated finite volume pressure gradient.

Global Functions

subroutine, public pressureforce_blk_afv(
    h h,
    tv tv,
    PFu PFu,
    PFv PFv,
    G G,
    GV GV,
    US US,
    CS CS,
    ALE_CSp ALE_CSp,
    p_atm p_atm,
    pbce pbce,
    eta eta
    )

Thin interface between the model and the Boussinesq and non-Boussinesq pressure force routines.

Parameters:

g

Ocean grid structure

gv

Vertical grid structure

us

A dimensional unit scaling type

h

Layer thickness [H ~> m or kg m-2]

tv

Thermodynamic variables

pfu

Zonal acceleration [L T-2 ~> m s-2]

pfv

Meridional acceleration [L T-2 ~> m s-2]

cs

Finite volume PGF control structure

ale_csp

ALE control structure

p_atm

The pressure at the ice-ocean or atmosphere-ocean interface [Pa].

pbce

The baroclinic pressure anomaly in each layer due to eta anomalies [m2 s-2 H-1 ~> m s-2 or m4 s-2 kg-1].

eta

The bottom mass used to calculate PFu and PFv [H ~> m or kg m-2], with any tidal contributions or compressibility compensation.

subroutine, public pressureforce_blk_afv_nonbouss(
    h h,
    tv tv,
    PFu PFu,
    PFv PFv,
    G G,
    GV GV,
    US US,
    CS CS,
    p_atm p_atm,
    pbce pbce,
    eta eta
    )

Non-Boussinesq analytically-integrated finite volume form of pressure gradient.

Determines the acceleration due to hydrostatic pressure forces, using the analytic finite volume form of the Pressure gradient, and does not make the Boussinesq approximation. This version uses code-blocking for threads.

To work, the following fields must be set outside of the usual (is:ie,js:je) range before this subroutine is called: h(isB:ie+1,jsB:je+1), T(isB:ie+1,jsB:je+1), and S(isB:ie+1,jsB:je+1).

Parameters:

g

Ocean grid structure

gv

Vertical grid structure

us

A dimensional unit scaling type

h

Layer thickness [H ~> kg m-2]

tv

Thermodynamic variables

pfu

Zonal acceleration [L T-2 ~> m s-2]

pfv

Meridional acceleration [L T-2 ~> m s-2]

cs

Finite volume PGF control structure

p_atm

The pressure at the ice-ocean or atmosphere-ocean interface [Pa].

pbce

The baroclinic pressure anomaly in each layer due to eta anomalies [m2 s-2 H-1 ~> m s-2 or m4 s-2 kg-1].

eta

The bottom mass used to calculate PFu and PFv [H ~> m or kg m-2], with any tidal contributions or compressibility compensation.

subroutine, public pressureforce_blk_afv_bouss(
    h h,
    tv tv,
    PFu PFu,
    PFv PFv,
    G G,
    GV GV,
    US US,
    CS CS,
    ALE_CSp ALE_CSp,
    p_atm p_atm,
    pbce pbce,
    eta eta
    )

Boussinesq analytically-integrated finite volume form of pressure gradient.

Determines the acceleration due to hydrostatic pressure forces, using the finite volume form of the terms and analytic integrals in depth, making the Boussinesq approximation. This version uses code-blocking for threads.

To work, the following fields must be set outside of the usual (is:ie,js:je) range before this subroutine is called: h(isB:ie+1,jsB:je+1), T(isB:ie+1,jsB:je+1), and S(isB:ie+1,jsB:je+1).

Parameters:

g

Ocean grid structure

gv

Vertical grid structure

us

A dimensional unit scaling type

h

Layer thickness [H ~> m or kg m-2]

tv

Thermodynamic variables

pfu

Zonal acceleration [L T-2 ~> m s-2]

pfv

Meridional acceleration [L T-2 ~> m s-2]

cs

Finite volume PGF control structure

ale_csp

ALE control structure

p_atm

The pressure at the ice-ocean or atmosphere-ocean interface [Pa].

pbce

The baroclinic pressure anomaly in each layer due to eta anomalies [m2 s-2 H-1 ~> m s-2 or m4 s-2 kg-1].

eta

The bottom mass used to calculate PFu and PFv [H ~> m or kg m-2], with any tidal contributions or compressibility compensation.

subroutine, public pressureforce_blk_afv_init(
    Time Time,
    G G,
    GV GV,
    US US,
    param_file param_file,
    diag diag,
    CS CS,
    tides_CSp tides_CSp
    )

Initializes the finite volume pressure gradient control structure.

Parameters:

time

Current model time

g

Ocean grid structure

gv

Vertical grid structure

us

A dimensional unit scaling type

param_file

Parameter file handles

diag

Diagnostics control structure

cs

Finite volume PGF control structure

tides_csp

Tides control structure

subroutine, public pressureforce_blk_afv_end(CS CS)

Deallocates the finite volume pressure gradient control structure.

Parameters:

cs

Blocked AFV pressure control structure that will be deallocated in this subroutine.