diff --git a/src/common/turb/mode_turb_ver_dyn_flux.F90 b/src/common/turb/mode_turb_ver_dyn_flux.F90
index c4f3cc9cad53fd5532c952b42cf843fa93d72327..680129b12928fe1a5925d23659136687f3cefcfd 100644
--- a/src/common/turb/mode_turb_ver_dyn_flux.F90
+++ b/src/common/turb/mode_turb_ver_dyn_flux.F90
@@ -1060,7 +1060,7 @@ IF(HTURBDIM=='3DIM') THEN
   ZWORK31D(IIJB:IIJE  ) = - ZFLXZ(IIJB:IIJE,IKB+D%NKL) &
                                     * ( ZWORK1(IIJB:IIJE,IKB+D%NKL) - ZWORK51D(IIJB:IIJE  ) &
                                     *   PDZY(IIJB:IIJE,IKB+D%NKL) ) &
-                                    / (0.5*(PDYY(IIJB:IIJE,IKB+D%NKL)+PDYY(IIJB:IIJE,IKE)))
+                                    / (0.5*(PDYY(IIJB:IIJE,IKB+D%NKL)+PDYY(IIJB:IIJE,IKB)))
   !$mnh_end_expand_array(JIJ=IIJB:IIJE)
   CALL MYF2D_PHY(D,ZWORK31D,ZWORK41D)
   ZA(IIJB:IIJE,IKB) = ZWORK41D(IIJB:IIJE)