From e3e9c69cc69c2659c281994cf3d7709e9d04a66d Mon Sep 17 00:00:00 2001 From: Alan Wallcraft Date: Tue, 4 Mar 2025 12:05:30 -0500 Subject: [PATCH] Update btstep negative eta warning (#844) * Update btstep negative eta warning In Boussinesq mode the eta below bathyT warning includes the location, but in non-Boussinesq mode the corresponding negative eta warning does not. Added the location to the negative eta warning. No answers are changed. * 4 space indenting of continuation lines --- src/core/MOM_barotropic.F90 | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/src/core/MOM_barotropic.F90 b/src/core/MOM_barotropic.F90 index 0948580e9f..e7b2787bcf 100644 --- a/src/core/MOM_barotropic.F90 +++ b/src/core/MOM_barotropic.F90 @@ -2536,8 +2536,11 @@ subroutine btstep(U_in, V_in, eta_in, dt, bc_accel_u, bc_accel_v, forces, pbce, else do j=js,je ; do i=is,ie if (eta(i,j) < 0.0) then + write(mesg,'(" at ", ES12.4, ES12.4, i7, i7)') & + G%geoLonT(i,j), G%geoLatT(i,j), i + G%HI%idg_offset, j + G%HI%jdg_offset if (err_count < 2) & - call MOM_error(WARNING, "btstep: negative eta in a non-Boussinesq barotropic solver.", all_print=.true.) + call MOM_error(WARNING, "btstep: negative eta in a non-Boussinesq barotropic solver "//& + trim(mesg), all_print=.true.) err_count = err_count + 1 endif enddo ; enddo