From b800565a06981f05870f843802393aba4124644b Mon Sep 17 00:00:00 2001 From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr> Date: Mon, 29 Apr 2019 14:32:11 +0200 Subject: [PATCH] Philippe 29/04/2019: gamma_inc(_low): bug correction in argument check (cherry picked from commit f234597fa739b139fa98e7a322c15b3e2097f38c) --- src/MNH/gamma_inc.f90 | 2 +- src/MNH/gamma_inc_low.f90 | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/src/MNH/gamma_inc.f90 b/src/MNH/gamma_inc.f90 index 083bf6d47..93e6247d0 100644 --- a/src/MNH/gamma_inc.f90 +++ b/src/MNH/gamma_inc.f90 @@ -87,7 +87,7 @@ REAL :: ZFPMIN=1.E-30 REAL :: ZAP,ZDEL,ZSUM REAL :: ZAN,ZB,ZC,ZD,ZH ! -IF( PX<0.0 .OR. PA>=0.0 ) call Print_msg(NVERB_FATAL,'GEN','GAMMA_INC','invalid arguments: PX<0.0 .OR. PA>=0.0') +IF( PX<0.0 .OR. PA<=0.0 ) call Print_msg(NVERB_FATAL,'GEN','GAMMA_INC','invalid arguments: PX<0.0 .OR. PA<=0.0') ! IF( (PX.LT.PA+1.0) ) THEN ZAP = PA diff --git a/src/MNH/gamma_inc_low.f90 b/src/MNH/gamma_inc_low.f90 index dabcd722a..bec06e6f0 100644 --- a/src/MNH/gamma_inc_low.f90 +++ b/src/MNH/gamma_inc_low.f90 @@ -109,7 +109,7 @@ ZS(5) = 2.9092306039 ! !* 1 Compute coefficients ! -IF( PX<0.0 .OR. PA>=0.0 ) call Print_msg(NVERB_FATAL,'GEN','GAMMA_INC_LOW','invalid arguments: PX<0.0 .OR. PA>=0.0') +IF( PX<0.0 .OR. PA<=0.0 ) call Print_msg(NVERB_FATAL,'GEN','GAMMA_INC_LOW','invalid arguments: PX<0.0 .OR. PA<=0.0') ! ZC(1) = 1.+ZP(1)*PA+ZP(2)*PA**2+ZP(3)*PA**3+ZP(4)*PA**4+ZP(5)*(EXP(-ZP(6)*PA)-1) ! -- GitLab