added few missing dummy values in case of zero absorption

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