From 9850319d79599d7390efcc175c383ad8e2cb30c2 Mon Sep 17 00:00:00 2001 From: He Wang Date: Mon, 24 Feb 2025 10:20:15 -0500 Subject: [PATCH] Bugfix: reference height offset for SAL using bpa Recover reference height offset (Z_ref) in calculating bottom pressure for self-attraction and loading in Boussinesq mode. --- src/core/MOM_PressureForce_FV.F90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/core/MOM_PressureForce_FV.F90 b/src/core/MOM_PressureForce_FV.F90 index a3272efb99..c497ddcfea 100644 --- a/src/core/MOM_PressureForce_FV.F90 +++ b/src/core/MOM_PressureForce_FV.F90 @@ -1242,7 +1242,7 @@ subroutine PressureForce_FV_Bouss(h, tv, PFu, PFv, G, GV, US, CS, ALE_CSp, p_atm if (CS%sal_use_bpa) then !$OMP parallel do default(shared) do j=Jsq,Jeq+1 ; do i=Isq,Ieq+1 - pbot(i,j) = pa(i,j,nz+1) - GxRho_ref * e(i,j,nz+1) + pbot(i,j) = pa(i,j,nz+1) - GxRho_ref * (e(i,j,nz+1) - G%Z_ref) enddo ; enddo call calc_SAL(pbot, e_sal, G, CS%SAL_CSp, tmp_scale=US%Z_to_m) else