mom_pressureforce_fv mom_pressureforce_fv::pressureforce_fv_cs subroutine, public subroutine, public mom_pressureforce_fv::pressureforce_fv_nonbouss (h, tv, PFu, PFv, G, GV, US, CS, ALE_CSp, p_atm, pbce, eta) pressureforce_fv_nonbouss 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 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.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). g Ocean grid structure gv Vertical grid structure us A dimensional unit scaling type h Layer thickness [H ~> kg/m2] 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 [R L2 T-2 ~> Pa]. pbce The baroclinic pressure anomaly in each layer due to eta anomalies [L2 T-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. mom_tidal_forcing::calc_tidal_forcing mom_error_handler::mom_error mom_pressureforce_mont::set_pbce_nonbouss mom_ale::ts_plm_edge_values mom_ale::ts_ppm_edge_values subroutine, public subroutine, public mom_pressureforce_fv::pressureforce_fv_bouss (h, tv, PFu, PFv, G, GV, US, CS, ALE_CSp, p_atm, pbce, eta) pressureforce_fv_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.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). g Ocean grid structure gv Vertical grid structure us A dimensional unit scaling type h Layer thickness [H ~> m] 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 [R L2 T-2 ~> Pa]. pbce The baroclinic pressure anomaly in each layer due to eta anomalies [L2 T-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. mom_tidal_forcing::calc_tidal_forcing mom_error_handler::mom_error mom_pressureforce_mont::set_pbce_bouss mom_ale::ts_plm_edge_values mom_ale::ts_ppm_edge_values subroutine, public subroutine, public mom_pressureforce_fv::pressureforce_fv_init (Time, G, GV, US, param_file, diag, CS, tides_CSp) pressureforce_fv_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. 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 mom_error_handler::mom_error subroutine, public subroutine, public mom_pressureforce_fv::pressureforce_fv_end (CS) pressureforce_fv_end CS CS Deallocates the finite volume pressure gradient control structure. cs Finite volume pressure control structure that will be deallocated in this subroutine. Finite volume pressure gradient (integrated by quadrature or analytically) Provides the Boussinesq and non-Boussinesq forms of horizontal accelerations due to pressure gradients using a vertically integrated finite volume form, as described by Adcroft et al., 2008. Integration in the vertical is made either by quadrature or analytically.This form eliminates the thermobaric instabilities that had been a problem with previous forms of the pressure gradient force calculation, as described by Hallberg, 2005.Adcroft, A., R. Hallberg, and M. Harrison, 2008: A finite volume discretization of the pressure gradient force using analytic integration. Ocean Modelling, 22, 106-113. http://doi.org/10.1016/j.ocemod.2008.02.001Hallberg, 2005: A thermobaric instability of Lagrangian vertical coordinate ocean models. Ocean Modelling, 8, 279-300. http://dx.doi.org/10.1016/j.ocemod.2004.01.001