PUBLIC INTERFACE ~ PUBLIC ROUTINES ~ NAMELIST

Module ocean_velocity_advect_mod

Contact:  S.M. Griffies
Reviewers: 
Change History: WebCVS Log

OVERVIEW

Velocity advective transport.

This module computes advection of velocity using one of the following advection schemes: 1/ second order centered 2/ first order upwind

OTHER MODULES USED

    diag_manager_mod
             fms_mod
             mpp_mod
   ocean_domains_mod
       ocean_obc_mod
 ocean_operators_mod
ocean_parameters_mod
     ocean_types_mod
 ocean_workspace_mod
      ocean_util_mod

PUBLIC INTERFACE

ocean_velocity_advect_init:
horz_advection_of_velocity:
horz_advection_centered:
horz_advection_upwind:
vert_advection_of_velocity:
vert_advection_centered:
vert_advection_upwind:

PUBLIC ROUTINES

  1. ocean_velocity_advect_init

    DESCRIPTION
    Initialize the velocity advection module.


  2. horz_advection_of_velocity

    DESCRIPTION
    Compute thickness weighted and density weighted acceleration (kg/m^3)*(m^2/s^2) due to horizontal (constant k-level) advection of velocity.


  3. horz_advection_centered

    DESCRIPTION
    Compute thickness weighted and density weighted acceleration (kg/m^3)*(m^2/s^2) due to horizontal (constant k-level) advection of velocity. Use second order centered method for either Bgrid or Cgrid.


  4. horz_advection_upwind

    DESCRIPTION
    Compute thickness weighted and density weighted acceleration (kg/m^3)*(m^2/s^2) due to horizontal (constant k-level) advection of velocity. Use first order upwind method. Not coded for Cgrid, since first order upwind method is rarely used.


  5. vert_advection_of_velocity

    DESCRIPTION
    Compute thickness weighted and density weighted acceleration (kg/m^3)*(m^2/s^2) due to vertical advection of velocity. Include vertical advection due to fresh water entering surface cells.


  6. vert_advection_centered

    DESCRIPTION
    Compute thickness weighted and density weighted acceleration (kg/m^3)*(m^2/s^2) due to vertical advection of velocity. Include vertical advection due to fresh water entering surface cells. Assume upme and uriver are the most fundamental fields, passed through the FMS coupler on the Bgrid velocity point. So use these to average onto the C-grid velocity points. Likewise, pme and river are fundamental, and sit on the T-grid. Use second order centered method here.


  7. vert_advection_upwind

    DESCRIPTION
    Compute thickness weighted and density weighted acceleration (kg/m^3)*(m^2/s^2) due to vertical advection of velocity. Include vertical advection due to fresh water entering surface cells. Use first order upwind method here. Not coded for Cgrid, since first order upwind method is rarely used.


NAMELIST

&ocean_velocity_advect_nml

debug_this_module
For debugging
[logical]
zero_velocity_advect_horz
For debugging, it is often useful to remove horizontal advection of velocity.
[logical]
zero_velocity_advect_vert
For debugging, it is often useful to remove vertical advection of velocity.
[logical]
velocity_advect_centered
For using the standard second order centered method for computing the advection of linear momentum. This is the default: velocity_advect_centered=.true.
[logical]
velocity_advect_upwind
For using the first order upwind method for computing the advection of linear momentum. Default: velocity_advect_upwind=.false.
[logical]

REFERENCES

  1. R.C. Pacanowski and S.M. Griffies The MOM3 Manual (1999)
  2. S.M. Griffies, M.J. Harrison, R.C. Pacanowski, and A. Rosati A Technical Guide to MOM4 (2004)
  3. R.C. Pacanowski and S.M. Griffies The MOM3 Manual (1999)
  4. S.M. Griffies: Elements of MOM (2012)
  5. Hundsdorder and Trompert (1994), "Method of lines and direct discretization: a comparison for linear advection", Applied Numerical Mathematics, pages 469--490.

top