ZSTAR(I) = -Z(I)
IF( Z(I) .LT. 0.15 ) PHIM(I) = 1. + ZSTAR(I)
IF( Z(I) .LT. 0.15 ) THEN
PHIH(I) =1.+ Z(I) * (0.5+ZSTAR(I)*(0.375+ZSTAR(I)*
_RL ZZ0(irun),Z(irun),Z2(irun),Z1(irun),Z0(irun)
Z(indx) = VZZ(I)
Z(I) = -18. * Z(I)
CALL PHI( Z,X1,Y1,IFLAG,IBIT )
IF ( Z(I) .LT. 0.013 ) ARG1(I) =
1 Z(I) * ( 0.25 - 0.09375 * Z(I) )
IF( Z(I) .LT. 0.0065 ) ARG1(I) =
1 Z(I) * ( 0.5 - 0.625 * Z(I) )
Z1(I) = Z(I)
IF ( Z(I) .GT. ZWM ) THEN
Z2(I) = Z(I)
IF ( (INTZ0(I).EQ.1) .OR. (Z(I).GT.ZWM) )
Y1DEN(I) = 1. + CM1 * ( X1NUM(I) * Z(I) )