diff --git a/src/OCE/LDF/ldfslp.F90 b/src/OCE/LDF/ldfslp.F90
index 4a85f6991d3fdfa7cb16b5cad538b18191e43d7d..d4804673554de52a04ebd0b8dd6de52161314407 100644
--- a/src/OCE/LDF/ldfslp.F90
+++ b/src/OCE/LDF/ldfslp.F90
@@ -684,6 +684,11 @@ CONTAINS
ALLOCATE( ah_wslp2(jpi,jpj,jpk) , akz(jpi,jpj,jpk) , STAT=ierr )
IF( ierr > 0 ) CALL ctl_stop( 'STOP', 'ldf_slp_init : unable to allocate ah_slp2 or akz' )
!
+ DO_3D( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpk )
+ akz (ji,jj,jk) = 0._wp
+ ah_wslp2(ji,jj,jk) = 0._wp
+ END_3D
+ !
IF( ln_traldf_triad ) THEN ! Griffies operator : triad of slopes
IF(lwp) WRITE(numout,*) ' ==>>> triad) operator (Griffies)'
ALLOCATE( triadi_g(jpi,jpj,jpk,0:1,0:1) , triadj_g(jpi,jpj,jpk,0:1,0:1) , &
diff --git a/src/OCE/TRA/traldf_iso.F90 b/src/OCE/TRA/traldf_iso.F90
index f93bf2dcda0500c21d02604f78db72b1e7d35209..8d84e76891bb958dec0bd0c1b561d7e8cb894f66 100644
--- a/src/OCE/TRA/traldf_iso.F90
+++ b/src/OCE/TRA/traldf_iso.F90
@@ -145,11 +145,6 @@ CONTAINS
IF(lwp) WRITE(numout,*) 'tra_ldf_iso : rotated laplacian diffusion operator on ', cdtype
IF(lwp) WRITE(numout,*) '~~~~~~~~~~~'
ENDIF
- !
- DO_3D_OVR( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpk )
- akz (ji,jj,jk) = 0._wp
- ah_wslp2(ji,jj,jk) = 0._wp
- END_3D
ENDIF
!
IF( .NOT. l_istiled .OR. ntile == 1 ) THEN ! Do only on the first tile
diff --git a/src/OCE/TRA/traldf_triad.F90 b/src/OCE/TRA/traldf_triad.F90
index 70eddec54f8399bc0e4531d1198843206aa64952..19039b88c12d763f483faede1d5fa33ec2c9a10e 100644
--- a/src/OCE/TRA/traldf_triad.F90
+++ b/src/OCE/TRA/traldf_triad.F90
@@ -147,11 +147,6 @@ CONTAINS
!!----------------------------------------------------------------------
!
IF( kpass == 1 ) THEN !== first pass only and whatever the tracer is ==!
- !
- DO_3D_OVR( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpk )
- akz (ji,jj,jk) = 0._wp
- ah_wslp2(ji,jj,jk) = 0._wp
- END_3D
!
DO kp = 0, 1 ! i-k triads
DO_3D_OVR( nn_hls-1, nn_hls-1, nn_hls-1, nn_hls-1, 1, jpkm1 )