Changeset 14808 for NEMO

Timestamp:

2021-05-07T14:00:45+02:00 (4 years ago)

Author:

jchanut

Message:

#2605: reproducibility issue with RK3: hdiv

File:

• NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/DYN/divhor.F90 (modified) (1 diff)

NEMO/branches/2021/dev_r14318_RK3_stage1/src/OCE/DYN/divhor.F90

@@ -101 +101 @@
       IF( ln_isf )   CALL isf_hdiv( kt, Kmm, hdiv )            !==  + ice-shelf mass exchange ==!
       !
+      CALL lbc_lnk( 'divhor', hdiv, 'T', 1._wp )   !   (no sign change)
+      !
 !!gm Patch before suppression of hdiv from all modules that use it
-      DO_3D( 0, 0, 0, 0, 1, jpkm1 )                            !==  e3t * Horizontal divergence  ==!
-         pe3divUh(ji,jj,jk) = hdiv(ji,jj,jk) * e3t(ji,jj,jk,Kmm)
-      END_3D
+!      DO_3D( 0, 0, 0, 0, 1, jpkm1 )                            !==  e3t * Horizontal divergence  ==!
+!         pe3divUh(ji,jj,jk) = hdiv(ji,jj,jk) * e3t(ji,jj,jk,Kmm)
+!      END_3D
+!JC: over whole domain, and after lbclnk on hdiv to prevent from reproducibility issues
+      DO jk=1, jpkm1
+         pe3divUh(:,:,jk) = hdiv(:,:,jk) * e3t(:,:,jk,Kmm)
+      END DO
 !!gm end
       !
-      CALL lbc_lnk( 'divhor', hdiv, 'T', 1._wp )   !   (no sign change)
       !
       IF( ln_timing )   CALL timing_stop('div_hor_RK3')