From 94212dba99b4eb4b4402c75baa2e3d85690bc86f Mon Sep 17 00:00:00 2001 From: Lorenzo Figini Date: Fri, 9 Jan 2015 15:05:20 +0000 Subject: [PATCH] added few missing dummy values in case of zero absorption --- src/gray.f | 8 +++++--- 1 file changed, 5 insertions(+), 3 deletions(-) diff --git a/src/gray.f b/src/gray.f index dcacff4..c59513b 100644 --- a/src/gray.f +++ b/src/gray.f @@ -4631,7 +4631,7 @@ c common/pstab/rpstab common/eqnn/nr,nz,npp,nintp c - ip=int((nintp-1)*rpsi+1) + ip=int((nintp-1)*rpsi+1) c if(ip.eq.0) ip=1 c if(ip.eq.nintp) ip=nintp-1 ip=min(max(1,ip),nintp-1) @@ -6580,6 +6580,8 @@ c of gaussian profile else ajmxfi=0.0d0 + ajphip=0.0d0 + dpdvp=0.0d0 dpdvmx=0.0d0 rhotjfi=1.0d0 rhotjav=1.0d0 @@ -6590,8 +6592,8 @@ c of gaussian profile drhotjava=0.0d0 drhotp=0.0d0 drhotpav=0.0d0 - ratjamx=0.0d0 - ratjbmx=0.0d0 + ratjamx=1.0d0 + ratjbmx=1.0d0 taumn=0 taumx=0 stmx=stf