diff --git a/src/LIB/RAD/ECMWF_RAD/swclr.f90 b/src/LIB/RAD/ECMWF_RAD/swclr.f90
index 45cdc930a7f942f6c1379d33ce42bc6b8df6e270..c9c801de301734a1be9a30208d12125bddd7b5d7 100644
--- a/src/LIB/RAD/ECMWF_RAD/swclr.f90
+++ b/src/LIB/RAD/ECMWF_RAD/swclr.f90
@@ -53,6 +53,7 @@ SUBROUTINE SWCLR &
 !        JJMorcrette 990128 : sunshine duration
 !        99-05-25   JJMorcrette    Revised aerosols
 !        JJMorcrette 001218 : 6 spectral intervals
+!        D.St Martin 11/2015 : bug on ZFACOA for NOVLP>= 5
    
 !     ------------------------------------------------------------------
 
@@ -273,7 +274,7 @@ JKLP1 = JKL + 1
 DO JL = KIDIA,KFDIA
 !++MODIF_MESONH
   IF (NOVLP.GE.5) THEN
-   ZFACOA = PTAUAZ(JL,JK)
+   ZFACOA = PTAUAZ(JL,JKL)
    ZCORAE = ZFACOA *  PSEC(JL)
   ELSE
    ZFACOA = _ONE_ - PPIZAZ(JL,JKL)*PCGAZ(JL,JKL)*PCGAZ(JL,JKL)
@@ -313,7 +314,7 @@ DO JK = 2 , KLEV
   DO JL = KIDIA,KFDIA
 !++MODIF_MESONH
     IF (NOVLP.GE.5) THEN
-     ZFACOA = PTAUAZ(JL,JK)
+     ZFACOA = PTAUAZ(JL,JKL)
      ZCORAE = ZFACOA *  PSEC(JL)
     ELSE
      ZFACOA = _ONE_ - PPIZAZ(JL,JKL)*PCGAZ(JL,JKL)*PCGAZ(JL,JKL)