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