1      SUBROUTINE FIG53A ( RV,Z,R)
2C
3C*********          SUBROUTINE TO COMPUTE FIGURE 4.1.5.2-53-A
4C
5      DIMENSION A(5),B(5),C(5),D(5),E(5),ZZ(5)
6      DATA A/8.98425E-03,4.50351E-03,6.0128E-03,1.07637E-02,8.48758E-03/
7      DATA B/-.138262,-.064239,-.08858,-.167056,-.130342/
8      DATA C/0.718213,0.263365,0.410583,0.893775,0.67405/
9      DATA D/-1.32167,-9.04994E-02,-0.489542,-1.80535,-1.22853/
10      DATA E/0.493821,-0.735445,-0.317567,1.02609,0.471355/
11      DATA ZZ/0.0,1.0,2.0,4.0,10.0/
12      IF(RV.EQ.0.0) GO TO 1070
13      X=ALOG10(RV)
14      M=1
15 1000 TEMP= ABS(Z-ZZ(M))
16      IF(Z.GE.ZZ(M).AND.Z.LT.ZZ(M+1)) GO TO 1030
17 1010 M=M+1
18      IF(M.GT.4) GO TO 1060
19 1020 GO TO 1000
20 1030 CONTINUE
21      YM=A(M)*X**5+B(M)*X**4+C(M)*X**3+D(M)*X**2+E(M)*X
22      IF(TEMP.LE.0.001) GO TO 1050
23 1040 N=M+1
24      YN=A(N)*X**5+B(N)*X**4+C(N)*X**3+D(N)*X**2+E(N)*X
25      R=YM+((Z-ZZ(M))/(ZZ(N)-ZZ(M)))*(YN-YM)
26      IF(R.LT.0.)R=0.
27      RETURN
28 1050 R=YM
29      IF(R.LT.0.)R=0.
30      RETURN
31 1060 R=A(M)*X**5+B(M)*X**4+C(M)*X**3+D(M)*X**2+E(M)*X
32      IF(R.LT.0.)R=0.
33      RETURN
34 1070 R=0.0
35      RETURN
36      END
37