mom_continuity_ppm mom_continuity_ppm::continuity_ppm_cs mom_continuity_ppm::loop_bounds_type integer integer mom_continuity_ppm::id_clock_update id_clock_update CPU time clock IDs. continuity_ppm continuity_ppm_init meridional_mass_flux zonal_mass_flux integer integer mom_continuity_ppm::id_clock_correct id_clock_correct CPU time clock IDs. continuity_ppm_init meridional_mass_flux zonal_mass_flux subroutine, public subroutine, public mom_continuity_ppm::continuity_ppm (u, v, hin, h, uh, vh, dt, G, GV, US, CS, uhbt, vhbt, OBC, visc_rem_u, visc_rem_v, u_cor, v_cor, BT_cont) continuity_ppm u u v v hin hin h h uh uh vh vh dt dt G G GV GV US US CS CS uhbt uhbt vhbt vhbt OBC OBC visc_rem_u visc_rem_u visc_rem_v visc_rem_v u_cor u_cor v_cor v_cor BT_cont BT_cont Time steps the layer thicknesses, using a monotonically limit, directionally split PPM scheme, based on Lin (1994). g The ocean's grid structure. cs Module's control structure. u Zonal velocity [L T-1 ~> m s-1]. v Meridional velocity [L T-1 ~> m s-1]. hin Initial layer thickness [H ~> m or kg m-2]. h Final layer thickness [H ~> m or kg m-2]. uh Zonal volume flux, u*h*dy [H L2 T-1 ~> m3 s-1 or kg s-1]. vh Meridional volume flux, v*h*dx [H L2 T-1 ~> m3 s-1 or kg s-1]. dt Time increment [T ~> s]. gv Vertical grid structure. us A dimensional unit scaling type uhbt The summed volume flux through zonal faces vhbt The summed volume flux through meridional faces obc Open boundaries control structure. visc_rem_u The fraction of zonal momentum originally in a layer that remains after a time-step of viscosity, and the fraction of a time-step's worth of a barotropic acceleration that a layer experiences after viscosity is applied. Non-dimensional between 0 (at the bottom) and 1 (far above the bottom). visc_rem_v The fraction of meridional momentum originally in a layer that remains after a time-step of viscosity, and the fraction of a time-step's worth of a barotropic acceleration that a layer experiences after viscosity is applied. Non-dimensional between 0 (at the bottom) and 1 (far above the bottom). u_cor The zonal velocities that give uhbt as the depth-integrated transport [L T-1 ~> m s-1]. v_cor The meridional velocities that give vhbt as the depth-integrated transport [L T-1 ~> m s-1]. bt_cont A structure with elements that describe the effective open face areas as a function of barotropic flow. id_clock_update meridional_mass_flux mom_error_handler::mom_error zonal_mass_flux subroutine subroutine mom_continuity_ppm::zonal_mass_flux (u, h_in, uh, dt, G, GV, US, CS, LB, uhbt, OBC, visc_rem_u, u_cor, BT_cont) zonal_mass_flux u u h_in h_in uh uh dt dt G G GV GV US US CS CS LB LB uhbt uhbt OBC OBC visc_rem_u visc_rem_u u_cor u_cor BT_cont BT_cont Calculates the mass or volume fluxes through the zonal faces, and other related quantities. g Ocean's grid structure. gv Ocean's vertical grid structure. u Zonal velocity [L T-1 ~> m s-1]. h_in Layer thickness used to calculate fluxes [H ~> m or kg m-2]. uh Volume flux through zonal faces = u*h*dy dt Time increment [T ~> s]. us A dimensional unit scaling type cs This module's control structure. lb Loop bounds structure. obc Open boundaries control structure. visc_rem_u The fraction of zonal momentum originally in a layer that remains after a time-step of viscosity, and the fraction of a time-step's worth of a barotropic acceleration that a layer experiences after viscosity is applied. Non-dimensional between 0 (at the bottom) and 1 (far above the bottom). uhbt The summed volume flux through zonal faces u_cor The zonal velocitiess (u with a barotropic correction) that give uhbt as the depth-integrated transport, m s-1. bt_cont A structure with elements that describe the effective open face areas as a function of barotropic flow. id_clock_correct id_clock_update mom_open_boundary::obc_direction_e ppm_reconstruction_x ratio_max set_zonal_bt_cont zonal_face_thickness zonal_flux_adjust zonal_flux_layer continuity_ppm subroutine subroutine mom_continuity_ppm::zonal_flux_layer (u, h, h_L, h_R, uh, duhdu, visc_rem, dt, G, US, j, ish, ieh, do_I, vol_CFL, OBC) zonal_flux_layer u u h h h_L h_L h_R h_R uh uh duhdu duhdu visc_rem visc_rem dt dt G G US US j j ish ish ieh ieh do_I do_I vol_CFL vol_CFL OBC OBC Evaluates the zonal mass or volume fluxes in a layer. g Ocean's grid structure. u Zonal velocity [L T-1 ~> m s-1]. visc_rem Both the fraction of the momentum originally in a layer that remains after a time-step of viscosity, and the fraction of a time-step's worth of a barotropic acceleration that a layer experiences after viscosity is applied. Non-dimensional between 0 (at the bottom) and 1 (far above the bottom). h Layer thickness [H ~> m or kg m-2]. h_l Left thickness [H ~> m or kg m-2]. h_r Right thickness [H ~> m or kg m-2]. uh Zonal mass or volume transport [H L2 T-1 ~> m3 s-1 or kg s-1]. duhdu Partial derivative of uh with u [H L ~> m2 or kg m-1]. dt Time increment [T ~> s]. us A dimensional unit scaling type j Spatial index. ish Start of index range. ieh End of index range. do_i Which i values to work on. vol_cfl If true, rescale the ratio of face areas to the cell areas when estimating the CFL number. obc Open boundaries control structure. mom_open_boundary::obc_direction_e set_zonal_bt_cont zonal_flux_adjust zonal_mass_flux subroutine subroutine mom_continuity_ppm::zonal_face_thickness (u, h, h_L, h_R, h_u, dt, G, US, LB, vol_CFL, marginal, visc_rem_u, OBC) zonal_face_thickness u u h h h_L h_L h_R h_R h_u h_u dt dt G G US US LB LB vol_CFL vol_CFL marginal marginal visc_rem_u visc_rem_u OBC OBC Sets the effective interface thickness at each zonal velocity point. g Ocean's grid structure. u Zonal velocity [L T-1 ~> m s-1]. h Layer thickness used to calculate fluxes [H ~> m or kg m-2]. h_l Left thickness in the reconstruction [H ~> m or kg m-2]. h_r Right thickness in the reconstruction [H ~> m or kg m-2]. h_u Thickness at zonal faces [H ~> m or kg m-2]. dt Time increment [T ~> s]. us A dimensional unit scaling type lb Loop bounds structure. vol_cfl If true, rescale the ratio of face areas to the cell areas when estimating the CFL number. marginal If true, report the marginal face thicknesses; otherwise report transport-averaged thicknesses. visc_rem_u Both the fraction of the momentum originally in a layer that remains after a time-step of viscosity, and the fraction of a time-step's worth of a barotropic acceleration that a layer experiences after viscosity is applied. Non-dimensional between 0 (at the bottom) and 1 (far above the bottom). obc Open boundaries control structure. mom_open_boundary::obc_direction_e zonal_mass_flux subroutine subroutine mom_continuity_ppm::zonal_flux_adjust (u, h_in, h_L, h_R, uhbt, uh_tot_0, duhdu_tot_0, du, du_max_CFL, du_min_CFL, dt, G, US, CS, visc_rem, j, ish, ieh, do_I_in, full_precision, uh_3d, OBC) zonal_flux_adjust u u h_in h_in h_L h_L h_R h_R uhbt uhbt uh_tot_0 uh_tot_0 duhdu_tot_0 duhdu_tot_0 du du du_max_CFL du_max_CFL du_min_CFL du_min_CFL dt dt G G US US CS CS visc_rem visc_rem j j ish ish ieh ieh do_I_in do_I_in full_precision full_precision uh_3d uh_3d OBC OBC Returns the barotropic velocity adjustment that gives the desired barotropic (layer-summed) transport. g Ocean's grid structure. u Zonal velocity [L T-1 ~> m s-1]. h_in Layer thickness used to calculate fluxes [H ~> m or kg m-2]. h_l Left thickness in the reconstruction [H ~> m or kg m-2]. h_r Right thickness in the reconstruction [H ~> m or kg m-2]. visc_rem Both the fraction of the momentum originally in a layer that remains after a time-step of viscosity, and the fraction of a time-step's worth of a barotropic acceleration that a layer experiences after viscosity is applied. Non-dimensional between 0 (at the bottom) and 1 (far above the bottom). uhbt The summed volume flux through zonal faces [H L2 T-1 ~> m3 s-1 or kg s-1]. du_max_cfl Maximum acceptable value of du [L T-1 ~> m s-1]. du_min_cfl Minimum acceptable value of du [L T-1 ~> m s-1]. uh_tot_0 The summed transport with 0 adjustment [H L2 T-1 ~> m3 s-1 or kg s-1]. duhdu_tot_0 The partial derivative of du_err with du at 0 adjustment [H L ~> m2 or kg m-1]. du The barotropic velocity adjustment [L T-1 ~> m s-1]. dt Time increment [T ~> s]. us A dimensional unit scaling type cs This module's control structure. j Spatial index. ish Start of index range. ieh End of index range. do_i_in A logical flag indicating which I values to work on. full_precision A flag indicating how carefully to iterate. The default is .true. (more accurate). uh_3d Volume flux through zonal faces = u*h*dy [H L2 T-1 ~> m3 s-1 or kg s-1]. obc Open boundaries control structure. zonal_flux_layer set_zonal_bt_cont zonal_mass_flux subroutine subroutine mom_continuity_ppm::set_zonal_bt_cont (u, h_in, h_L, h_R, BT_cont, uh_tot_0, duhdu_tot_0, du_max_CFL, du_min_CFL, dt, G, US, CS, visc_rem, visc_rem_max, j, ish, ieh, do_I) set_zonal_bt_cont u u h_in h_in h_L h_L h_R h_R BT_cont BT_cont uh_tot_0 uh_tot_0 duhdu_tot_0 duhdu_tot_0 du_max_CFL du_max_CFL du_min_CFL du_min_CFL dt dt G G US US CS CS visc_rem visc_rem visc_rem_max visc_rem_max j j ish ish ieh ieh do_I do_I Sets a structure that describes the zonal barotropic volume or mass fluxes as a function of barotropic flow to agree closely with the sum of the layer's transports. g Ocean's grid structure. u Zonal velocity [L T-1 ~> m s-1]. h_in Layer thickness used to calculate fluxes [H ~> m or kg m-2]. h_l Left thickness in the reconstruction [H ~> m or kg m-2]. h_r Right thickness in the reconstruction [H ~> m or kg m-2]. bt_cont A structure with elements that describe the effective open face areas as a function of barotropic flow. uh_tot_0 The summed transport with 0 adjustment [H L2 T-1 ~> m3 s-1 or kg s-1]. duhdu_tot_0 The partial derivative of du_err with du at 0 adjustment [H L ~> m2 or kg m-1]. du_max_cfl Maximum acceptable value of du [L T-1 ~> m s-1]. du_min_cfl Minimum acceptable value of du [L T-1 ~> m s-1]. dt Time increment [T ~> s]. us A dimensional unit scaling type cs This module's control structure. visc_rem Both the fraction of the momentum originally in a layer that remains after a time-step of viscosity, and the fraction of a time-step's worth of a barotropic acceleration that a layer experiences after viscosity is applied. Non-dimensional between 0 (at the bottom) and 1 (far above the bottom). visc_rem_max Maximum allowable visc_rem. j Spatial index. ish Start of index range. ieh End of index range. do_i A logical flag indicating which I values to work on. zonal_flux_adjust zonal_flux_layer zonal_mass_flux subroutine subroutine mom_continuity_ppm::meridional_mass_flux (v, h_in, vh, dt, G, GV, US, CS, LB, vhbt, OBC, visc_rem_v, v_cor, BT_cont) meridional_mass_flux v v h_in h_in vh vh dt dt G G GV GV US US CS CS LB LB vhbt vhbt OBC OBC visc_rem_v visc_rem_v v_cor v_cor BT_cont BT_cont Calculates the mass or volume fluxes through the meridional faces, and other related quantities. g Ocean's grid structure. gv Ocean's vertical grid structure. v Meridional velocity [L T-1 ~> m s-1]. h_in Layer thickness used to calculate fluxes [H ~> m or kg m-2]. vh Volume flux through meridional faces = v*h*dx [H m2 s-1 ~> m3 s-1 or kg s-1]. dt Time increment [T ~> s]. us A dimensional unit scaling type cs This module's control structure.G lb Loop bounds structure. obc Open boundary condition type specifies whether, where, and what open boundary conditions are used. visc_rem_v Both the fraction of the momentum vhbt The summed volume flux through meridional faces [H L2 T-1 ~> m3 s-1 or kg s-1]. v_cor The meridional velocitiess (v with a barotropic correction) that give vhbt as the depth-integrated transport [L T-1 ~> m s-1]. bt_cont A structure with elements that describe the effective open face areas as a function of barotropic flow. id_clock_correct id_clock_update merid_face_thickness merid_flux_layer meridional_flux_adjust mom_open_boundary::obc_direction_n ppm_reconstruction_y ratio_max set_merid_bt_cont continuity_ppm subroutine subroutine mom_continuity_ppm::merid_flux_layer (v, h, h_L, h_R, vh, dvhdv, visc_rem, dt, G, US, J, ish, ieh, do_I, vol_CFL, OBC) merid_flux_layer v v h h h_L h_L h_R h_R vh vh dvhdv dvhdv visc_rem visc_rem dt dt G G US US J J ish ish ieh ieh do_I do_I vol_CFL vol_CFL OBC OBC Evaluates the meridional mass or volume fluxes in a layer. g Ocean's grid structure. v Meridional velocity [L T-1 ~> m s-1]. visc_rem Both the fraction of the momentum originally in a layer that remains after a time-step of viscosity, and the fraction of a time-step's worth of a barotropic acceleration that a layer experiences after viscosity is applied. Non-dimensional between 0 (at the bottom) and 1 (far above the bottom). h Layer thickness used to calculate fluxes, [H ~> m or kg m-2]. h_l Left thickness in the reconstruction [H ~> m or kg m-2]. h_r Right thickness in the reconstruction [H ~> m or kg m-2]. vh Meridional mass or volume transport [H L2 T-1 ~> m3 s-1 or kg s-1]. dvhdv Partial derivative of vh with v [H L ~> m2 or kg m-1]. dt Time increment [T ~> s]. us A dimensional unit scaling type j Spatial index. ish Start of index range. ieh End of index range. do_i Which i values to work on. vol_cfl If true, rescale the ratio of face areas to the cell areas when estimating the CFL number. obc Open boundaries control structure. mom_open_boundary::obc_direction_n meridional_flux_adjust meridional_mass_flux set_merid_bt_cont subroutine subroutine mom_continuity_ppm::merid_face_thickness (v, h, h_L, h_R, h_v, dt, G, US, LB, vol_CFL, marginal, visc_rem_v, OBC) merid_face_thickness v v h h h_L h_L h_R h_R h_v h_v dt dt G G US US LB LB vol_CFL vol_CFL marginal marginal visc_rem_v visc_rem_v OBC OBC Sets the effective interface thickness at each meridional velocity point. g Ocean's grid structure. v Meridional velocity [L T-1 ~> m s-1]. h Layer thickness used to calculate fluxes, [H ~> m or kg m-2]. h_l Left thickness in the reconstruction, [H ~> m or kg m-2]. h_r Right thickness in the reconstruction, [H ~> m or kg m-2]. h_v Thickness at meridional faces, [H ~> m or kg m-2]. dt Time increment [T ~> s]. lb Loop bounds structure. us A dimensional unit scaling type vol_cfl If true, rescale the ratio of face areas to the cell areas when estimating the CFL number. marginal If true, report the marginal face thicknesses; otherwise report transport-averaged thicknesses. visc_rem_v Both the fraction of the momentum originally in a layer that remains after a time-step of viscosity, and the fraction of a time-step's worth of a barotropic acceleration that a layer experiences after viscosity is applied. Non-dimensional between 0 (at the bottom) and 1 (far above the bottom). obc Open boundaries control structure. mom_open_boundary::obc_direction_n meridional_mass_flux subroutine subroutine mom_continuity_ppm::meridional_flux_adjust (v, h_in, h_L, h_R, vhbt, vh_tot_0, dvhdv_tot_0, dv, dv_max_CFL, dv_min_CFL, dt, G, US, CS, visc_rem, j, ish, ieh, do_I_in, full_precision, vh_3d, OBC) meridional_flux_adjust v v h_in h_in h_L h_L h_R h_R vhbt vhbt vh_tot_0 vh_tot_0 dvhdv_tot_0 dvhdv_tot_0 dv dv dv_max_CFL dv_max_CFL dv_min_CFL dv_min_CFL dt dt G G US US CS CS visc_rem visc_rem j j ish ish ieh ieh do_I_in do_I_in full_precision full_precision vh_3d vh_3d OBC OBC Returns the barotropic velocity adjustment that gives the desired barotropic (layer-summed) transport. g Ocean's grid structure. v Meridional velocity [L T-1 ~> m s-1]. h_in Layer thickness used to calculate fluxes [H ~> m or kg m-2]. h_l Left thickness in the reconstruction [H ~> m or kg m-2]. h_r Right thickness in the reconstruction [H ~> m or kg m-2]. visc_rem Both the fraction of the momentum originally in a layer that remains after a time-step of viscosity, and the fraction of a time-step's worth of a barotropic acceleration that a layer experiences after viscosity is applied. Non-dimensional between 0 (at the bottom) and 1 (far above the bottom). vhbt The summed volume flux through meridional faces dv_max_cfl Maximum acceptable value of dv [L T-1 ~> m s-1]. dv_min_cfl Minimum acceptable value of dv [L T-1 ~> m s-1]. vh_tot_0 The summed transport with 0 adjustment [H L2 T-1 ~> m3 s-1 or kg s-1]. dvhdv_tot_0 The partial derivative of dv_err with dv at 0 adjustment [H L ~> m2 or kg m-1]. dv The barotropic velocity adjustment [L T-1 ~> m s-1]. dt Time increment [T ~> s]. us A dimensional unit scaling type cs This module's control structure. j Spatial index. ish Start of index range. ieh End of index range. do_i_in A flag indicating which I values to work on. full_precision A flag indicating how carefully to iterate. The default is .true. (more accurate). vh_3d Volume flux through meridional obc Open boundaries control structure. merid_flux_layer meridional_mass_flux set_merid_bt_cont subroutine subroutine mom_continuity_ppm::set_merid_bt_cont (v, h_in, h_L, h_R, BT_cont, vh_tot_0, dvhdv_tot_0, dv_max_CFL, dv_min_CFL, dt, G, US, CS, visc_rem, visc_rem_max, j, ish, ieh, do_I) set_merid_bt_cont v v h_in h_in h_L h_L h_R h_R BT_cont BT_cont vh_tot_0 vh_tot_0 dvhdv_tot_0 dvhdv_tot_0 dv_max_CFL dv_max_CFL dv_min_CFL dv_min_CFL dt dt G G US US CS CS visc_rem visc_rem visc_rem_max visc_rem_max j j ish ish ieh ieh do_I do_I Sets of a structure that describes the meridional barotropic volume or mass fluxes as a function of barotropic flow to agree closely with the sum of the layer's transports. g Ocean's grid structure. v Meridional velocity [L T-1 ~> m s-1]. h_in Layer thickness used to calculate fluxes, [H ~> m or kg m-2]. h_l Left thickness in the reconstruction, [H ~> m or kg m-2]. h_r Right thickness in the reconstruction, [H ~> m or kg m-2]. bt_cont A structure with elements that describe the effective open face areas as a function of barotropic flow. vh_tot_0 The summed transport with 0 adjustment [H L2 T-1 ~> m3 s-1 or kg s-1]. dvhdv_tot_0 The partial derivative of du_err with dv at 0 adjustment [H L ~> m2 or kg m-1]. dv_max_cfl Maximum acceptable value of dv [L T-1 ~> m s-1]. dv_min_cfl Minimum acceptable value of dv [L T-1 ~> m s-1]. dt Time increment [T ~> s]. us A dimensional unit scaling type cs This module's control structure. visc_rem Both the fraction of the momentum originally in a layer that remains after a time-step of viscosity, and the fraction of a time-step's worth of a barotropic acceleration that a layer experiences after viscosity is applied. Non-dimensional between 0 (at the bottom) and 1 (far above the bottom). visc_rem_max Maximum allowable visc_rem. j Spatial index. ish Start of index range. ieh End of index range. do_i A logical flag indicating which I values to work on. merid_flux_layer meridional_flux_adjust meridional_mass_flux subroutine subroutine mom_continuity_ppm::ppm_reconstruction_x (h_in, h_L, h_R, G, LB, h_min, monotonic, simple_2nd, OBC) ppm_reconstruction_x h_in h_in h_L h_L h_R h_R G G LB LB h_min h_min monotonic monotonic simple_2nd simple_2nd OBC OBC Calculates left/right edge values for PPM reconstruction. g Ocean's grid structure. h_in Layer thickness [H ~> m or kg m-2]. h_l Left thickness in the reconstruction, [H ~> m or kg m-2]. h_r Right thickness in the reconstruction, [H ~> m or kg m-2]. lb Active loop bounds structure. h_min The minimum thickness that can be obtained by a concave parabolic fit. monotonic If true, use the Colella & Woodward monotonic limiter. Otherwise use a simple positive-definite limiter. simple_2nd If true, use the arithmetic mean thicknesses as the default edge values for a simple 2nd order scheme. obc Open boundaries control structure. mom_error_handler::mom_error mom_open_boundary::obc_direction_e mom_open_boundary::obc_direction_w ppm_limit_cw84 ppm_limit_pos zonal_mass_flux subroutine subroutine mom_continuity_ppm::ppm_reconstruction_y (h_in, h_L, h_R, G, LB, h_min, monotonic, simple_2nd, OBC) ppm_reconstruction_y h_in h_in h_L h_L h_R h_R G G LB LB h_min h_min monotonic monotonic simple_2nd simple_2nd OBC OBC Calculates left/right edge values for PPM reconstruction. g Ocean's grid structure. h_in Layer thickness [H ~> m or kg m-2]. h_l Left thickness in the reconstruction, [H ~> m or kg m-2]. h_r Right thickness in the reconstruction, [H ~> m or kg m-2]. lb Active loop bounds structure. h_min The minimum thickness that can be obtained by a concave parabolic fit. monotonic If true, use the Colella & Woodward monotonic limiter. Otherwise use a simple positive-definite limiter. simple_2nd If true, use the arithmetic mean thicknesses as the default edge values for a simple 2nd order scheme. obc Open boundaries control structure. mom_error_handler::mom_error mom_open_boundary::obc_direction_n mom_open_boundary::obc_direction_s ppm_limit_cw84 ppm_limit_pos meridional_mass_flux subroutine subroutine mom_continuity_ppm::ppm_limit_pos (h_in, h_L, h_R, h_min, G, iis, iie, jis, jie) ppm_limit_pos h_in h_in h_L h_L h_R h_R h_min h_min G G iis iis iie iie jis jis jie jie This subroutine limits the left/right edge values of the PPM reconstruction to give a reconstruction that is positive-definite. Here this is reinterpreted as giving a constant thickness if the mean thickness is less than h_min, with a minimum of h_min otherwise. g Ocean's grid structure. h_in Layer thickness [H ~> m or kg m-2]. h_l Left thickness in the reconstruction [H ~> m or kg m-2]. h_r Right thickness in the reconstruction [H ~> m or kg m-2]. h_min The minimum thickness that can be obtained by a concave parabolic fit. iis Start of i index range. iie End of i index range. jis Start of j index range. jie End of j index range. ppm_reconstruction_x ppm_reconstruction_y subroutine subroutine mom_continuity_ppm::ppm_limit_cw84 (h_in, h_L, h_R, G, iis, iie, jis, jie) ppm_limit_cw84 h_in h_in h_L h_L h_R h_R G G iis iis iie iie jis jis jie jie This subroutine limits the left/right edge values of the PPM reconstruction according to the monotonic prescription of Colella and Woodward, 1984. g Ocean's grid structure. h_in Layer thickness [H ~> m or kg m-2]. h_l Left thickness in the reconstruction, [H ~> m or kg m-2]. h_r Right thickness in the reconstruction, [H ~> m or kg m-2]. iis Start of i index range. iie End of i index range. jis Start of j index range. jie End of j index range. ppm_reconstruction_x ppm_reconstruction_y real function real function mom_continuity_ppm::ratio_max (a, b, maxrat) ratio_max a a b b maxrat maxrat Return the maximum ratio of a/b or maxrat. a Numerator b Denominator maxrat Maximum value of ratio. Return value. meridional_mass_flux zonal_mass_flux subroutine, public subroutine, public mom_continuity_ppm::continuity_ppm_init (Time, G, GV, US, param_file, diag, CS) continuity_ppm_init Time Time G G GV GV US US param_file param_file diag diag CS CS Initializes continuity_ppm_cs. time The current model time. g The ocean's grid structure. gv Vertical grid structure. us A dimensional unit scaling type param_file A structure indicating the open file to parse for model parameter values. diag A structure that is used to regulate diagnostic output. cs Module's control structure. This include declares and sets the variable "version". mom_cpu_clock::cpu_clock_id id_clock_correct id_clock_update mom_error_handler::mom_error integer function, public integer function, public mom_continuity_ppm::continuity_ppm_stencil (CS) continuity_ppm_stencil CS CS continuity_PPM_stencil returns the continuity solver stencil size cs Module's control structure. The continuity solver stencil size with the current settings. subroutine, public subroutine, public mom_continuity_ppm::continuity_ppm_end (CS) continuity_ppm_end CS CS Destructor for continuity_ppm_cs. cs Module's control structure. mom_continuity::continuity_end Solve the layer continuity equation using the PPM method for layer fluxes. This module contains the subroutines that advect layer thickness. The scheme here uses a Piecewise-Parabolic method with a positive definite limiter.