Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

+Use RHO_PGF_REF for the pressure gradient forces #392

Merged
merged 3 commits into from
Jul 20, 2023
Merged

+Use RHO_PGF_REF for the pressure gradient forces #392

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
0bd4613
+Use RHO_PGF_REF for the pressure gradient forces
Hallberg-NOAA Jun 30, 2023
c0bdb6d
Merge branch 'dev/gfdl' into add_RHO_PGF_REF
adcroft Jul 18, 2023
09c99f7
Merge branch 'dev/gfdl' into add_RHO_PGF_REF
marshallward Jul 20, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
13 changes: 6 additions & 7 deletions src/core/MOM_PressureForce_FV.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -519,7 +519,7 @@ subroutine PressureForce_FV_Bouss(h, tv, PFu, PFv, G, GV, US, CS, ALE_CSp, p_atm
if (associated(ALE_CSp)) use_ALE = CS%reconstruct .and. use_EOS

h_neglect = GV%H_subroundoff
dz_neglect = GV%H_subroundoff * GV%H_to_Z
dz_neglect = GV%dZ_subroundoff
I_Rho0 = 1.0 / GV%Rho0
G_Rho0 = GV%g_Earth / GV%Rho0
rho_ref = CS%Rho0
Expand Down Expand Up @@ -827,12 +827,11 @@ subroutine PressureForce_FV_init(Time, G, GV, US, param_file, diag, CS, tides_CS

mdl = "MOM_PressureForce_FV"
call log_version(param_file, mdl, version, "")
call get_param(param_file, mdl, "RHO_0", CS%Rho0, &
"The mean ocean density used with BOUSSINESQ true to "//&
"calculate accelerations and the mass for conservation "//&
"properties, or with BOUSSINSEQ false to convert some "//&
"parameters from vertical units of m to kg m-2.", &
units="kg m-3", default=1035.0, scale=US%kg_m3_to_R)
call get_param(param_file, mdl, "RHO_PGF_REF", CS%Rho0, &
"The reference density that is subtracted off when calculating pressure "//&
"gradient forces. Its inverse is subtracted off of specific volumes when "//&
"in non-Boussinesq mode. The default is RHO_0.", &
units="kg m-3", default=GV%Rho0*US%R_to_kg_m3, scale=US%kg_m3_to_R)
call get_param(param_file, mdl, "TIDES", CS%tides, &
"If true, apply tidal momentum forcing.", default=.false.)
call get_param(param_file, "MOM", "USE_REGRIDDING", use_ALE, &
Expand Down