_RL ZSTAR(N),E1(N),E2(N),TEMP1(N)
ZSTAR(I) = 100. * Z(I) - 14.
IF( Z(I) .GT. 2. ) ZSTAR(I) = TEMP1(I)
IF( Z(I).GT.1.78e10 ) ZSTAR(I) = 384.9999
I1(I) = ZSTAR(I)
TEMP1(I) = ZSTAR(I) - I1(I)
ZSTAR(I) = -Z(I)
IF( Z(I) .LT. 0.15 ) PHIM(I) = 1. + ZSTAR(I)
2 *(0.25+ZSTAR(I)*(0.09375+ZSTAR(I)*
3 (0.03125+0.00732422 * ZSTAR(I))))
PHIH(I) =1.+ Z(I) * (0.5+ZSTAR(I)*(0.375+ZSTAR(I)*
1 (0.5+ZSTAR(I)*(0.8203125+ZSTAR(I)*
2 (1.5+2.93262*ZSTAR(I))))))