diff --git a/src/MNH/turb.f90 b/src/MNH/turb.f90
index 3cf06ade78de9d27a03e5a00081163200e2b9a6d..921cdd662446e728f5593fefd9b4ee70dd24a5dc 100644
--- a/src/MNH/turb.f90
+++ b/src/MNH/turb.f90
@@ -610,7 +610,7 @@ IF (HTURBLEN=='BL89' .OR. HTURBLEN=='RM17' .OR. ORMC01) THEN
   ZRM(:,:,:,:) = PRT(:,:,:,:)
 END IF
 !
-!
+ZTRH(:, :, : ) = XUNDEF
 !
 !----------------------------------------------------------------------------
 !
diff --git a/src/MNH/turb_ver_thermo_flux.f90 b/src/MNH/turb_ver_thermo_flux.f90
index b5b7355cea1e24394b3623be1c1258a5b18d18a7..770c4cf96d51798455a2e30723040c80be28559d 100644
--- a/src/MNH/turb_ver_thermo_flux.f90
+++ b/src/MNH/turb_ver_thermo_flux.f90
@@ -794,6 +794,7 @@ ZFLXZ(:,:,KKA) = ZFLXZ(:,:,IKB)
   PWTH(:,:,IKB)=0.5*(ZFLXZ(:,:,IKB)+ZFLXZ(:,:,IKB+KKL))
   PWTH(:,:,KKA)=0.5*(ZFLXZ(:,:,KKA)+ZFLXZ(:,:,KKA+KKL))
   PWTH(:,:,IKE)=PWTH(:,:,IKE-KKL)
+  PWTH(:, :, IKE + 1 : ) = XUNDEF
 !$acc end kernels
 !
 IF ( OTURB_FLX .AND. OCLOSE_OUT ) THEN
@@ -1193,6 +1194,7 @@ IF (KRR /= 0) THEN
   PWRC(:,:,IKB)=0.5*(ZFLXZ(:,:,IKB)+ZFLXZ(:,:,IKB+KKL))
   PWRC(:,:,KKA)=0.5*(ZFLXZ(:,:,KKA)+ZFLXZ(:,:,KKA+KKL))
   PWRC(:,:,IKE)=PWRC(:,:,IKE-KKL)
+  PWRC(:, :, IKE + 1 : ) = XUNDEF
 !$acc end kernels
   !
   !