diff --git a/src/MNH/mode_les_diachro.f90 b/src/MNH/mode_les_diachro.f90
index 536ab8fc28d797002c0a36d127cb808d7ae1b24d..e681a460c1de837c6ffda2a7e418be2ecd410723 100644
--- a/src/MNH/mode_les_diachro.f90
+++ b/src/MNH/mode_les_diachro.f90
@@ -412,7 +412,8 @@ SUBROUTINE LES_TIME_AVG(PWORK6,tpdates,KRESP)
 ! this routine computes time averaging
 !
 ! Modifications:
-!  03/2018     (P.Wautelet)   replace ADD_FORECAST_TO_DATE by DATETIME_CORRECTDATE
+!  P. Wautelet    03/2018: replace ADD_FORECAST_TO_DATE by DATETIME_CORRECTDATE
+!  P. Wautelet 28/08/2020: bugfix: allocate ZWORK6 at the correct dimension (3rd one)
 !
 use modd_time,      only: tdtseg
 USE MODD_LES
@@ -456,7 +457,7 @@ END IF
 deallocate( tpdates )
 !
 allocate( tpdates( iavg ) )
-ALLOCATE (ZWORK6(SIZE(PWORK6,1),SIZE(PWORK6,2),NLES_K,IAVG,SIZE(PWORK6,5),SIZE(PWORK6,6)))
+ALLOCATE (ZWORK6(SIZE(PWORK6,1),SIZE(PWORK6,2),SIZE(PWORK6,3),IAVG,SIZE(PWORK6,5),SIZE(PWORK6,6)))
 !
 ZWORK6(:,:,:,:,:,:) = 0.
 !
@@ -501,7 +502,7 @@ DO JAVG=1,IAVG
 END DO
 !
 DEALLOCATE(PWORK6)
-ALLOCATE(PWORK6(SIZE(ZWORK6,1),SIZE(ZWORK6,2),NLES_K,IAVG,SIZE(ZWORK6,5),SIZE(ZWORK6,6)))
+ALLOCATE(PWORK6(SIZE(ZWORK6,1),SIZE(ZWORK6,2),SIZE(ZWORK6,3),IAVG,SIZE(ZWORK6,5),SIZE(ZWORK6,6)))
 PWORK6 = ZWORK6
 DEALLOCATE(ZWORK6)