Changeset 14978 for NEMO/trunk/src/OCE/TRA/traadv_qck.F90

Timestamp:

2021-06-11T15:21:08+02:00 (3 years ago)

Author:

hadcv

Message:

#2640: correct tra_adv_qck results on the north fold

File:

• NEMO/trunk/src/OCE/TRA/traadv_qck.F90 (modified) (1 diff)

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

@@ -246 +246 @@
          IF (nn_hls==1) CALL lbc_lnk( 'traadv_qck', zfc(:,:,:), 'T', 1.0_wp , zfd(:,:,:), 'T', 1.0_wp, ld4only= .TRUE. )   ! Lateral boundary conditions
 
+         ! Correct zfd on northfold after lbc_lnk; see #2640
+         IF( nn_hls == 1 .AND. l_IdoNFold .AND. ntej == Nje0 ) THEN
+            DO jk = 1, jpkm1
+               WHERE( tmask_i(ntsi:ntei,ntej:jpj) == 0._wp ) zfd(ntsi:ntei,ntej:jpj,jk) = zfc(ntsi:ntei,ntej:jpj,jk)
+            END DO
+         ENDIF
          !
          ! Horizontal advective fluxes