Skip to content
Snippets Groups Projects
Commit 1a147abb authored by Gaelle Tanguy's avatar Gaelle Tanguy Committed by WAUTELET Philippe
Browse files

D.St Martin 11/2015 : bug on ZFACOA for NOVLP>= 5

parent adf7fa21
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
src/LIB/RAD/ECMWF_RAD/swclr.f90
+ 3
− 2
View file @ 1a147abb
...@@ -53,6 +53,7 @@ SUBROUTINE SWCLR & ...@@ -53,6 +53,7 @@ SUBROUTINE SWCLR &
! JJMorcrette 990128 : sunshine duration ! JJMorcrette 990128 : sunshine duration
! 99-05-25 JJMorcrette Revised aerosols ! 99-05-25 JJMorcrette Revised aerosols
! JJMorcrette 001218 : 6 spectral intervals ! JJMorcrette 001218 : 6 spectral intervals
! D.St Martin 11/2015 : bug on ZFACOA for NOVLP>= 5
! ------------------------------------------------------------------ ! ------------------------------------------------------------------
...@@ -273,7 +274,7 @@ JKLP1 = JKL + 1 ...@@ -273,7 +274,7 @@ JKLP1 = JKL + 1
DO JL = KIDIA,KFDIA DO JL = KIDIA,KFDIA
!++MODIF_MESONH !++MODIF_MESONH
IF (NOVLP.GE.5) THEN IF (NOVLP.GE.5) THEN
ZFACOA = PTAUAZ(JL,JK) ZFACOA = PTAUAZ(JL,JKL)
ZCORAE = ZFACOA * PSEC(JL) ZCORAE = ZFACOA * PSEC(JL)
ELSE ELSE
ZFACOA = _ONE_ - PPIZAZ(JL,JKL)*PCGAZ(JL,JKL)*PCGAZ(JL,JKL) ZFACOA = _ONE_ - PPIZAZ(JL,JKL)*PCGAZ(JL,JKL)*PCGAZ(JL,JKL)
...@@ -313,7 +314,7 @@ DO JK = 2 , KLEV ...@@ -313,7 +314,7 @@ DO JK = 2 , KLEV
DO JL = KIDIA,KFDIA DO JL = KIDIA,KFDIA
!++MODIF_MESONH !++MODIF_MESONH
IF (NOVLP.GE.5) THEN IF (NOVLP.GE.5) THEN
ZFACOA = PTAUAZ(JL,JK) ZFACOA = PTAUAZ(JL,JKL)
ZCORAE = ZFACOA * PSEC(JL) ZCORAE = ZFACOA * PSEC(JL)
ELSE ELSE
ZFACOA = _ONE_ - PPIZAZ(JL,JKL)*PCGAZ(JL,JKL)*PCGAZ(JL,JKL) ZFACOA = _ONE_ - PPIZAZ(JL,JKL)*PCGAZ(JL,JKL)*PCGAZ(JL,JKL)
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment

Documentation | Mentions légales | CGU | Accessibilité : non conforme