1 SUBROUTINE zs1s2(ZRR, ZRI, S1R, S1I, S2R, S2I, NZ, ASCLE, ALIM,
17 DOUBLE PRECISION AA, ALIM, ALN, ASCLE, AS1, AS2, C1I, C1R, S1DI,
18 * s1dr, s1i, s1r, s2i, s2r, zeroi, zeror, zri, zrr, xzabs
20 DATA zeror,zeroi / 0.0d0 , 0.0d0 /
24 IF (s1r.EQ.0.0d0 .AND. s1i.EQ.0.0d0) go to 10
25 IF (as1.EQ.0.0d0) go to 10
26 aln = -zrr - zrr + dlog(as1)
32 IF (aln.LT.(-alim)) go to 10
33 CALL
xzlog(s1dr, s1di, c1r, c1i, idum)
36 CALL
xzexp(c1r, c1i, s1r, s1i)
41 IF (aa.GT.ascle)
RETURN
subroutine xzlog(AR, AI, BR, BI, IERR)
subroutine xzexp(AR, AI, BR, BI)
subroutine zs1s2(ZRR, ZRI, S1R, S1I, S2R, S2I, NZ, ASCLE, ALIM, IUF)