From 1a147abbbca515bcd5a72cc29de4cf62d554039e Mon Sep 17 00:00:00 2001
From: Gaelle Tanguy <gaelle.tanguy@meteo.fr>
Date: Thu, 26 Nov 2015 14:48:37 +0000
Subject: [PATCH] D.St Martin 11/2015 : bug on ZFACOA for NOVLP>= 5
---
src/LIB/RAD/ECMWF_RAD/swclr.f90 | 5 +++--
1 file changed, 3 insertions(+), 2 deletions(-)
diff --git a/src/LIB/RAD/ECMWF_RAD/swclr.f90 b/src/LIB/RAD/ECMWF_RAD/swclr.f90
index 45cdc930a..c9c801de3 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)
--
GitLab