Changeset 12366

Timestamp:

2020-02-11T12:17:24+01:00 (5 years ago)

Author:

jchanut

Message:

#2317, revert to the correct sign for asselin term on e3t and runoffs

Location:

NEMO/trunk/src/OCE

Files:

• DYN/dynnxt.F90 (modified) (1 diff)
• TRA/tranxt.F90 (modified) (1 diff)

NEMO/trunk/src/OCE/DYN/dynnxt.F90

@@ -232 +232 @@
             IF ( ln_rnf ) THEN
                DO jk = 1, jpkm1
-                  e3t_b(:,:,jk) = e3t_b(:,:,jk) - zcoef * ( rnf_b(:,:) - rnf(:,:) ) * tmask(:,:,jk) & 
+                  e3t_b(:,:,jk) = e3t_b(:,:,jk) + zcoef * ( rnf_b(:,:) - rnf(:,:) ) * tmask(:,:,jk) & 
                                 &                       * e3t_n(:,:,jk) /  ( ht_n(:,:) + 1._wp - ssmask(:,:) )
                END DO

NEMO/trunk/src/OCE/TRA/tranxt.F90

@@ -314 +314 @@
                   zscale = zfact2 * e3t_n(ji,jj,jk) / ( ht_n(ji,jj) + 1._wp - ssmask(ji,jj) )
                   ze3t_f = ze3t_f - zscale * ( emp_b(ji,jj) - emp(ji,jj) )
-                  IF ( ll_rnf ) ze3t_f = ze3t_f - zscale * (    rnf_b(ji,jj) -    rnf(ji,jj) )
+                  IF ( ll_rnf ) ze3t_f = ze3t_f + zscale * (    rnf_b(ji,jj) -    rnf(ji,jj) )
                   IF ( ll_isf ) ze3t_f = ze3t_f - zscale * ( fwfisf_b(ji,jj) - fwfisf(ji,jj) )