diff --git a/src/MNH/ini_lesn.f90 b/src/MNH/ini_lesn.f90
index 40fd4bb1d7e6eabfb499160f32fc0698c045cc35..34b5101e40877c22b46f6bedfdab70c8e156e62d 100644
--- a/src/MNH/ini_lesn.f90
+++ b/src/MNH/ini_lesn.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 2000-2020 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2000-2021 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
 !MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
@@ -35,6 +35,7 @@
 !!  Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O
 !!                     02/2019 (C. Lac) Add rain fraction as a LES diagnostic
 !  P. Wautelet 12/08/2020: bugfix: use NUNDEF instead of XUNDEF for integer variables
+!  P. Wautelet 04/01/2021: bugfix: nles_k was used instead of nspectra_k for a loop index
 !! --------------------------------------------------------------------------
 !
 !*      0. DECLARATIONS
@@ -427,7 +428,7 @@ IF (ANY(XSPECTRA_ALTITUDES(:)/=XUNDEF)) THEN
   ALLOCATE(NKLIN_SPEC   (SIZE(XZZ,1),SIZE(XZZ,2),NSPECTRA_K))
   !
   ALLOCATE(ZZ_SPEC      (SIZE(XZZ,1),SIZE(XZZ,2),NSPECTRA_K))
-  DO JK=1,NLES_K
+  DO JK=1,NSPECTRA_K
     DO JJ=1,SIZE(XZZ,2)
       DO JI=1,SIZE(XZZ,1)
         ZZ_SPEC(JI,JJ,JK) = XSPECTRA_ALTITUDES(JK)