From 3b82d4f5fd0164396b914c0c2d27786a6c3c3ccc Mon Sep 17 00:00:00 2001 From: Robert Hallberg Date: Fri, 14 Dec 2018 19:19:38 -0500 Subject: [PATCH] Corrected uh_tot size in zonal_mass_flux Corrected the array size of hu_tot in zonal_mass_flux. This changes answers with symmetric memory, but not without symmetric memory, and it is an unambiguous bug that has been recently introduced. All answers now reproduce the reference solutions. --- src/SIS_continuity.F90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/SIS_continuity.F90 b/src/SIS_continuity.F90 index 4c1170d7..c6511cba 100644 --- a/src/SIS_continuity.F90 +++ b/src/SIS_continuity.F90 @@ -716,7 +716,7 @@ subroutine zonal_mass_flux(u, dt, G, IG, CS, LB, h_in, uh, htot_in, uh_tot) optional, intent(out) :: uh !< Category volume flux through zonal faces = u*h*dy, H m2 s-1. real, dimension(SZI_(G),SZJ_(G)), & optional, intent(in) :: htot_in !< Total thicknesses used to calculate the fluxes, in H. - real, dimension(SZI_(G),SZJB_(G)), & + real, dimension(SZIB_(G),SZJ_(G)), & optional, intent(out) :: uh_tot !< Total mass flux through zonal faces = u*htot*dy, H m2 s-1. ! This subroutine calculates the mass or volume fluxes through the zonal ! faces, and other related quantities.