* installing *source* package ‘GeoModels’ ...
** package ‘GeoModels’ successfully unpacked and MD5 sums checked
** using staged installation
** libs
using C compiler: ‘Apple clang version 14.0.3 (clang-1403.0.22.14.1)’
using Fortran compiler: ‘GNU Fortran (GCC) 12.2.0’
using SDK: ‘MacOSX11.3.sdk’
clang -arch x86_64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I/opt/R/x86_64/include    -fPIC  -falign-functions=64 -Wall -g -O2  -c CompositeLikelihood2.c -o CompositeLikelihood2.o
clang -arch x86_64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I/opt/R/x86_64/include    -fPIC  -falign-functions=64 -Wall -g -O2  -c CompositeLikelihood2_ani.c -o CompositeLikelihood2_ani.o
clang -arch x86_64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I/opt/R/x86_64/include    -fPIC  -falign-functions=64 -Wall -g -O2  -c CompositeLikelihoodCond2.c -o CompositeLikelihoodCond2.o
clang -arch x86_64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I/opt/R/x86_64/include    -fPIC  -falign-functions=64 -Wall -g -O2  -c CompositeLikelihoodCond2_ani.c -o CompositeLikelihoodCond2_ani.o
clang -arch x86_64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I/opt/R/x86_64/include    -fPIC  -falign-functions=64 -Wall -g -O2  -c CorrelationFunction.c -o CorrelationFunction.o
clang -arch x86_64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I/opt/R/x86_64/include    -fPIC  -falign-functions=64 -Wall -g -O2  -c Distributions.c -o Distributions.o
clang -arch x86_64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I/opt/R/x86_64/include    -fPIC  -falign-functions=64 -Wall -g -O2  -c GeoModels_init.c -o GeoModels_init.o
clang -arch x86_64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I/opt/R/x86_64/include    -fPIC  -falign-functions=64 -Wall -g -O2  -c Godambe.c -o Godambe.o
Godambe.c:653:9: warning: variable 'nsub1' set but not used [-Wunused-but-set-variable]
    int nsub1 =0;
        ^
clang -arch x86_64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I/opt/R/x86_64/include    -fPIC  -falign-functions=64 -Wall -g -O2  -c Gradient.c -o Gradient.o
/opt/gfortran/bin/gfortran -arch x86_64  -fPIC  -Wall -g -O2  -c HypergeometricFunctions.f -o HypergeometricFunctions.o
HypergeometricFunctions.f:233:53:

  233 |           IF (CDABS(CMPLX(AR,AI)*CMPLX(XR,XI)/(CMPLX(CR,CI)*CNT))
      |                                                     1
Warning: Conversion from REAL(8) to default-kind COMPLEX(4) at (1) might lose precision, consider using the KIND argument [-Wconversion]
HypergeometricFunctions.f:233:39:

  233 |           IF (CDABS(CMPLX(AR,AI)*CMPLX(XR,XI)/(CMPLX(CR,CI)*CNT))
      |                                       1
Warning: Conversion from REAL(8) to default-kind COMPLEX(4) at (1) might lose precision, consider using the KIND argument [-Wconversion]
HypergeometricFunctions.f:233:26:

  233 |           IF (CDABS(CMPLX(AR,AI)*CMPLX(XR,XI)/(CMPLX(CR,CI)*CNT))
      |                          1
Warning: Conversion from REAL(8) to default-kind COMPLEX(4) at (1) might lose precision, consider using the KIND argument [-Wconversion]
HypergeometricFunctions.f:302:12:

  302 |       EDIFF=DNINT(A(L+1)-B(L+1))
      |            1
Warning: Possible change of value in conversion from REAL(8) to INTEGER(4) at (1) [-Wconversion]
HypergeometricFunctions.f:989:16:

  989 |       A = CMPLX(ARE, AIM)
      |                1
Warning: Conversion from REAL(8) to default-kind COMPLEX(4) at (1) might lose precision, consider using the KIND argument [-Wconversion]
HypergeometricFunctions.f:990:16:

  990 |       B = CMPLX(BRE, BIM)
      |                1
Warning: Conversion from REAL(8) to default-kind COMPLEX(4) at (1) might lose precision, consider using the KIND argument [-Wconversion]
HypergeometricFunctions.f:992:18:

  992 |         Z = CMPLX(ZRE(I), ZIM(I))
      |                  1
Warning: Conversion from REAL(8) to default-kind COMPLEX(4) at (1) might lose precision, consider using the KIND argument [-Wconversion]
clang -arch x86_64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I/opt/R/x86_64/include    -fPIC  -falign-functions=64 -Wall -g -O2  -c TB.c -o TB.o
clang -arch x86_64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I/opt/R/x86_64/include    -fPIC  -falign-functions=64 -Wall -g -O2  -c Utility.c -o Utility.o
/opt/gfortran/bin/gfortran -arch x86_64  -fPIC  -Wall -g -O2  -c mvn.f -o mvn.o
/opt/gfortran/bin/gfortran -arch x86_64  -fPIC  -Wall -g -O2  -c mvndstpack.f -o mvndstpack.o
mvn.f:970:3:

  970 |  10      END DO
      |   1
Warning: Label 10 at (1) defined but not used [-Wunused-label]
mvn.f:72:15:

   72 |       INFORM = MVNNIT( N, CORREL, LOWER, UPPER, INFIN, INFIS, D, E )
      |               1
Warning: Possible change of value in conversion from REAL(8) to INTEGER(4) at (1) [-Wconversion]
mvn.f:209:38:

  209 |         IF ( MINCLS .LT. 0 ) SBRGNS = WORK(LENWRK)
      |                                      1
Warning: Possible change of value in conversion from REAL(8) to INTEGER(4) at (1) [-Wconversion]
mvn.f:267:21:

  267 |             RGNCLS = RGNCLS*(MESH(DIVAXN)+1)/MESH(DIVAXN)
      |                     1
Warning: Possible change of value in conversion from REAL(8) to INTEGER(4) at (1) [-Wconversion]
mvn.f:268:21:

  268 |             NWRGNS = NWRGNS*(MESH(DIVAXN)+1)/MESH(DIVAXN)
      |                     1
Warning: Possible change of value in conversion from REAL(8) to INTEGER(4) at (1) [-Wconversion]
mvn.f:1268:24:

 1268 |                POINTS = PONTRS(SUBTMP)
      |                        1
Warning: Possible change of value in conversion from REAL(8) to INTEGER(4) at (1) [-Wconversion]
mvn.f:1269:24:

 1269 |                POINTP = PONTRS(SUBTMP+1)
      |                        1
Warning: Possible change of value in conversion from REAL(8) to INTEGER(4) at (1) [-Wconversion]
mvn.f:1277:21:

 1277 |             POINTS = PONTRS(SUBTMP)
      |                     1
Warning: Possible change of value in conversion from REAL(8) to INTEGER(4) at (1) [-Wconversion]
mvn.f:1297:21:

 1297 |             POINTS = PONTRS(SUBTMP)
      |                     1
Warning: Possible change of value in conversion from REAL(8) to INTEGER(4) at (1) [-Wconversion]
mvn.f:316:15:

  316 |          TOP = PONTRS(1)
      |               1
Warning: Possible change of value in conversion from REAL(8) to INTEGER(4) at (1) [-Wconversion]
mvn.f:323:21:

  323 |             RGNCLS = RGNCLS*MESH(I)
      |                     1
Warning: Possible change of value in conversion from REAL(8) to INTEGER(4) at (1) [-Wconversion]
mvn.f:328:18:

  328 |          RGNCLS = RGNCLS*(MESH(DIVAXN)+1)/MESH(DIVAXN)
      |                  1
Warning: Possible change of value in conversion from REAL(8) to INTEGER(4) at (1) [-Wconversion]
mvn.f:1573:3:

 1573 |  10      END DO
      |   1
Warning: Label 10 at (1) defined but not used [-Wunused-label]
mvn.f:1418:15:

 1418 |       INFORM = MVTNIT( N, NU, CORREL, LOWER, UPPER, INFIN, INFIS, D, E )
      |               1
Warning: Possible change of value in conversion from REAL(8) to INTEGER(4) at (1) [-Wconversion]
mvndstpack.f:116:18:

  116 |          INFORM = MVNDNT(N, CORREL, LOWER, UPPER, INFIN, INFIS, D, E)
      |                  1
Warning: Possible change of value in conversion from REAL(8) to INTEGER(4) at (1) [-Wconversion]
mvndstpack.f:704:14:

  704 |          JP = J + MVNUNI()*( NK + 1 - J )
      |              1
Warning: Possible change of value in conversion from REAL(8) to INTEGER(4) at (1) [-Wconversion]
mvndstpack.f:940:72:

  940 |       END
      |                                                                        ^
Warning: ‘__result_bvnmvn’ may be used uninitialized [-Wmaybe-uninitialized]
mvndstpack.f:898:38:

  898 |       DOUBLE PRECISION FUNCTION BVNMVN( LOWER, UPPER, INFIN, CORREL )
      |                                      ^
note: ‘__result_bvnmvn’ was declared here
1 warning generated.
/opt/gfortran/bin/gfortran -arch x86_64  -fPIC  -Wall -g -O2  -c specfun.f -o specfun.o
specfun.f:121:72:

  121 | 10            WB0=WB
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:131:72:

  131 | 15            CF0=CF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:157:72:

  157 | 20            CF=CF+CR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:162:72:

  162 | 25            CG=CG+CR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:190:72:

  190 |            DO 10 J=0,N
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 10 at (1)
specfun.f:207:72:

  207 |            DO 15 J=2,N
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 15 at (1)
specfun.f:212:72:

  212 |            DO 20 I=2,M
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 20 at (1)
specfun.f:229:72:

  229 | 25            QF1=QF0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:231:72:

  231 | 30            QM(0,K)=Q0*QM(0,K)/QF0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:238:72:

  238 | 35            QF1=QF0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:241:72:

  241 | 40            QM(1,K)=Q10*QM(1,K)/QF0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:245:72:

  245 |               DO 45 I=0,M-2
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 45 at (1)
specfun.f:254:72:

  254 | 50         QD(0,J)=LS*J*(QM(0,J-1)-X*QM(0,J))/XS
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:256:72:

  256 |         DO 55 I=1,M
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 55 at (1)
specfun.f:285:72:

  285 |         DO 10 J=0,M
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 10 at (1)
specfun.f:287:72:

  287 | 10         CPD(J,I)=(0.0D0,0.0D0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:293:72:

  293 | 15            CPD(0,I)=0.5D0*I*(I+1)*X**(I+1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:295:72:

  295 |            DO 20 I=1,M
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 20 at (1)
specfun.f:320:72:

  320 | 25         CPM(I,I)=(2.0D0*I-1.0D0)*ZQ*CPM(I-1,I-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:323:72:

  323 | 30         CPM(I,I+1)=(2.0D0*I+1.0D0)*Z*CPM(I,I)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:325:72:

  325 |         DO 35 J=I+2,N
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 35 at (1)
specfun.f:333:72:

  333 | 40         CPD(0,J)=LS*J*(Z*CPM(0,J)-CPM(0,J-1))/ZS
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:335:72:

  335 |         DO 45 J=I,N
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 45 at (1)
specfun.f:488:72:

  488 | 40               P(K)=P1(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:506:72:

  506 | 45         L0=L2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:530:72:

  530 | 10         U(J)=C*C
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:532:72:

  532 | 15         V(J)=(2.0*J-1.0-IP)*(2.0*(J-M)-IP)+M*(M-1.0)-CV
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:534:72:

  534 | 20         W(J)=(2.0*J-IP)*(2.0*J+1.0-IP)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:544:72:

  544 | 25                  R1=R1*(I+M-J)/J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:560:72:

  560 | 45                  R1=R1*(I+M-J)/J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:574:72:

  574 | 65         BK(K)=(BK(K)-BK(K-1)*U(K))/T
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 65 at (1)
specfun.f:576:72:

  576 | 70         BK(K)=BK(K)-W(K)*BK(K+1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 70 at (1)
specfun.f:621:72:

  621 | 10        SW=SU0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:627:72:

  627 | 20        SW=SD0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:634:72:

  634 | 30           SD1=SD1+DN(K)*QD(J)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:642:72:

  642 | 35              R1=(M+J1)*R1
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:645:72:

  645 | 40              R2=J2*R2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:650:72:

  650 | 45              SF=SF+R3
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:662:72:

  662 | 50              SD=SD+R4
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:673:72:

  673 | 60           SW=SU2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:679:72:

  679 | 70           SW=SD2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 70 at (1)
specfun.f:714:72:

  714 | 20         BN(M)=R1*R2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:737:72:

  737 | 10               R=R*(J+M-K)/J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:738:72:

  738 | 20            S=S-R*BN(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:739:72:

  739 | 30         BN(M)=S
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:741:72:

  741 | 40         BN(M)=0.0D0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:765:72:

  765 | 10               SK=SK+CK(K+1)*CK(L-K+1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:766:72:

  766 | 15            S=S+SK*AP(I-L+1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:767:72:

  767 | 20      AP(I+1)=-R*S
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:772:72:

  772 | 25            R=R*(2.0D0*K+IP)*(2.0D0*K-1.0D0+IP)/(2.0D0*K)**2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:773:72:

  773 | 30         QS0=QS0+AP(M-L+1)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:1029:72:

 1029 |  10        F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:1068:72:

 1068 | 10         F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:1214:72:

 1214 | 55            G0=G0+R0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 55 at (1)
specfun.f:1219:72:

 1219 | 60            G1=G1+R1
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
clang -arch x86_64 -I"/Library/Frameworks/R.framework/Resources/include" -DNDEBUG   -I/opt/R/x86_64/include    -fPIC  -falign-functions=64 -Wall -g -O2  -c weightedleastsquare.c -o weightedleastsquare.o
specfun.f:1262:72:

 1262 | 10               CF(K)=CF(K)*CT2+A(I)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:1263:72:

 1263 | 15            CF(K)=CF(K)*CT**K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:1267:72:

 1267 | 20            CSJ=CSJ+CF(K)*VR**K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:1272:72:

 1272 | 25            CSY=CSY+(-1)**K*CF(K)*VR**K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:1310:72:

 1310 | 10         R0=R0*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:1318:72:

 1318 | 15         SW=SUC
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:1336:72:

 1336 | 50         SW=R2F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:1359:72:

 1359 | 60         SW=SUD
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:1389:72:

 1389 | 10            S=S+1.0D0/K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:1394:72:

 1394 | 20            S=S+1.0/(2.0D0*K-1.0D0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:1400:72:

 1400 | 30               S=S+1.0D0/(XA+K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:1527:72:

 1527 | 10         PD(K)=0.0D0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:1551:72:

 1551 | 20         PM0=PMK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:1559:72:

 1559 | 25         PM1=PM2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:1562:72:

 1562 | 30         PD(K)=(K*X*PM(K)-(K+M)*PM(K-1))/(X*X-1.0D0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:1565:72:

 1565 | 35         PD(K)=(-1)**M*PD(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:1596:72:

 1596 | 10            W=ER
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:1604:72:

 1604 | 20            ER=ER+R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:1622:72:

 1622 | 25            W1=ER2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:1631:72:

 1631 | 35            W2=EI2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:1748:72:

 1748 | 20         GR=GR*X+G(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:1809:72:

 1809 | 10               R=R*(2.0D0*M-2.0D0*K+J)/J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:1810:72:

 1810 | 20            S=S+R*EN(2*K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:1811:72:

 1811 | 30         EN(2*M)=-S
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:1847:72:

 1847 | 10         F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:1899:72:

 1899 | 25            XA0=XA
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:1902:72:

 1902 | 30            XS=XS+2.0D0*BJ(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:1905:72:

 1905 | 35            BJ(K)=BJ(K)/XS
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:1910:72:

 1910 | 40            XG1=XG1+BJ(K)*XR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:1915:72:

 1915 | 45            XG2=XG2+BJ(K)*XR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:1925:72:

 1925 | 50            XF=XF+XR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:1930:72:

 1930 | 55            XG=XG+XR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 55 at (1)
specfun.f:1980:72:

 1980 | 30            A1=AF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:1985:72:

 1985 | 35            TI=TI+A(K)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:2014:72:

 2014 | 10            CQD(K)=(1.0D+300,0.0D0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:2030:72:

 2030 | 15            CQF1=CQF2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:2043:72:

 2043 | 20            CQF1=CQF0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:2045:72:

 2045 | 25            CQN(K)=CQN(K)*CQ0/CQF0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:2049:72:

 2049 | 30         CQD(K)=(K*Z*CQN(K)-K*CQN(K-1))/(Z*Z-1.0D0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:2197:72:

 2197 | 20            ER=ER+R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:2297:72:

 2297 | 20         EN(M)=R1*R2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:2334:72:

 2334 | 10               CV(IC)=4.0D0*(IC-1.0D0)**2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:2337:72:

 2337 | 15               CV(IC)=(2.0D0*IC-1.0D0)**2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:2340:72:

 2340 | 20               CV(IC)=4.0D0*IC*IC
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:2351:72:

 2351 | 25               F(I)=Q*Q
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:2359:72:

 2359 | 30               F(I)=Q*Q
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:2365:72:

 2365 | 35               F(I)=Q*Q
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:2379:72:

 2379 | 45            H(I)=XB
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:2490:72:

 2490 | 10            QD(K)=1.0D+300
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:2506:72:

 2506 | 15            Q1=QF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:2536:72:

 2536 | 40            QF1=QF0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:2539:72:

 2539 | 45            QD(K)=K*(QN(K-1)-X*QN(K))/(1.0D0-X*X)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:2570:72:

 2570 | 10         G0=G
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:2572:72:

 2572 |            DO 15 J=1,K
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 15 at (1)
specfun.f:2621:72:

 2621 | 20            TTI=TTI+C(K)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:2644:72:

 2644 | 35            TTK=TTK+C(K)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
mvn.f:1121:31:

 1121 |       INTEGER S, IXCHNG, LXCHNG, I, L
      |                               ^
Warning: ‘lxchng’ may be used uninitialized [-Wmaybe-uninitialized]
specfun.f:2696:72:

 2696 | 25            DL(K)=-0.5D0*X/(VK+1.0D0)*UK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:2711:72:

 2711 | 30            PX=PX+RP
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:2717:72:

 2717 | 35            QX=QX+RQ
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:2752:72:

 2752 | 45            VL(K)=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:2767:72:

 2767 | 50            F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:2775:72:

 2775 | 55            VL(K)=CS*VL(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 55 at (1)
specfun.f:2782:72:

 2782 | 65         R0=2.0D0*(J+V0+1)/X*R0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 65 at (1)
specfun.f:2923:72:

 2923 | 10         W(I)=(K+1.0D0)*(K+2.0D0)*CS/((GK2+2.0D0)*GK3)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:2927:72:

 2927 | 15            T=V(M-L)-W(M-L+1)*U(M-L)/T
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:2928:72:

 2928 | 20         RK(K)=-U(K)/T
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:2932:72:

 2932 | 25         DN(K)=DF(1)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:2943:72:

 2943 | 35         DN(K)=RK(K)*DN(K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:2946:72:

 2946 | 40         R1=R1*(J+0.5D0*(N+M+IP))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:2950:72:

 2950 | 45         R=R*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:2957:72:

 2957 | 50         SW=SU0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:2961:72:

 2961 | 60         R2=2.0D0*C*R2*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:2964:72:

 2964 | 65         R3=R3*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 65 at (1)
specfun.f:2970:72:

 2970 | 75         R4=4.0D0*R4*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 75 at (1)
specfun.f:2973:72:

 2973 | 80         R5=R5*(J+M)/C
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 80 at (1)
specfun.f:3010:72:

 3010 | 15            P=P*(Z-X(I))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:3017:72:

 3017 | 20            F1=PF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:3124:72:

 3124 | 10            CDY(K)=(1.0D+300,0.0D0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:3163:72:

 3163 | 30               CPZ=CPZ+CRP
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:3169:72:

 3169 | 35               CQZ=CQZ+CRQ
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:3211:72:

 3211 | 60               CS0=CS0+CR0*W0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:3219:72:

 3219 | 65               CS1=CS1+CR1*(2.0D0*W1+1.0D0/(K+1.0D0))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 65 at (1)
specfun.f:3247:72:

 3247 | 70            CF1=CF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 70 at (1)
specfun.f:3261:72:

 3261 | 75            CF1=CF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 75 at (1)
specfun.f:3265:72:

 3265 | 80            CBJ(K)=CS*CBJ(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 80 at (1)
specfun.f:3269:72:

 3269 | 85         CDJ(K)=-(K+V0)/Z*CBJ(K)+CBJ(K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 85 at (1)
specfun.f:3294:72:

 3294 | 100        CH1=CH0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 100 at (1)
specfun.f:3302:72:

 3302 | 105        CH1=CH0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 105 at (1)
specfun.f:3319:72:

 3319 | 110        CYL1=CYLK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 110 at (1)
specfun.f:3326:72:

 3326 | 115        CYL2=CYLK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 115 at (1)
specfun.f:3334:72:

 3334 | 130        CDY(K)=CBY(K-1)-(K+V0)/Z*CBY(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 130 at (1)
specfun.f:3381:72:

 3381 | 10            CDY(K)=(1.0D+300,0.0D0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:3414:72:

 3414 | 25            CPZ=CPZ+CRP
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:3420:72:

 3420 | 30            CQZ=CQZ+CRQ
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:3451:72:

 3451 | 45               CS0=CS0+CR0*W0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:3468:72:

 3468 | 50         CF1=CF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:3471:72:

 3471 | 55         CBJ(K)=CS*CBJ(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 55 at (1)
specfun.f:3495:72:

 3495 | 70         CDJ(K)=-(K+V0)/Z*CBJ(K)+CBJ(K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 70 at (1)
specfun.f:3498:72:

 3498 | 75         CDY(K)=CBY(K-1)-(K+V0)/Z*CBY(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 75 at (1)
specfun.f:3611:72:

 3611 | 45            Q0=Q0+B(K)*X**(-2*K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:3620:72:

 3620 | 50            Q1=Q1+B1(K)*X**(-2*K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:3792:72:

 3792 | 20            TI=TI+A(K)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:3810:72:

 3810 | 25            TW=TK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:3817:72:

 3817 | 35            TK=TK+A(K)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:3858:72:

 3858 | 10            DY(K)=1.0D+300
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:3899:72:

 3899 | 30               PX=PX+RP
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:3905:72:

 3905 | 35               QX=QX+RQ
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:3931:72:

 3931 | 45            F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:3946:72:

 3946 | 50            F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:3953:72:

 3953 | 55            BJ(K)=CS*BJ(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 55 at (1)
specfun.f:3956:72:

 3956 | 60         DJ(K)=-(K+V0)/X*BJ(K)+BJ(K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:3988:72:

 3988 | 80               CS0=CS0+R0*W0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 80 at (1)
specfun.f:3996:72:

 3996 | 85               CS1=CS1+R1*(2.0D0*W1+1.0D0/(K+1.0D0))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 85 at (1)
specfun.f:4006:72:

 4006 | 90         BYV1=BYVK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 90 at (1)
specfun.f:4009:72:

 4009 | 95         DY(K)=-(K+V0)/X*BY(K)+BY(K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 95 at (1)
specfun.f:4041:72:

 4041 |  10           DY(K) = 1.0D+300
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:4046:72:

 4046 |  40           DJ(K)=BJ(K-1)-K/X*BJ(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:4049:72:

 4049 |  50           DY(K)=BY(K-1)-K*BY(K)/X
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:4080:72:

 4080 | 10            BY(K-NMIN)=-1.0D+300
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:4109:72:

 4109 | 15            F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:4112:72:

 4112 | 20            BJ(K-NMIN)=BJ(K-NMIN)/S0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:4137:72:

 4137 | 25            Q0=Q0+B(K)*X**(-2*K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:4148:72:

 4148 | 30            Q1=Q1+B1(K)*X**(-2*K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:4157:72:

 4157 | 35            BJ1=BJK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:4165:72:

 4165 | 45         BY1=BYK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:4193:72:

 4193 | 15            P=P*(Z-X(I))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:4201:72:

 4201 | 20            F1=PF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:4210:72:

 4210 | 35            Q=Q+WP
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:4217:72:

 4217 | 45         W(N+1-NR)=W(NR)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:4318:72:

 4318 | 10            DY(K)=1.0D+300
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:4338:72:

 4338 | 20            BJ1=BJK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:4353:72:

 4353 | 30            F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:4360:72:

 4360 | 40            BJ(K)=CS*BJ(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:4363:72:

 4363 | 50         DJ(K)=BJ(K-1)-K/X*BJ(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:4370:72:

 4370 | 60         F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:4372:72:

 4372 | 70         DY(K)=BY(K-1)-K*BY(K)/X
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 70 at (1)
specfun.f:4427:72:

 4427 | 15            PD1=PDF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:4445:72:

 4445 | 20               PD1=PD
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:4459:72:

 4459 | 25               F0=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:4472:72:

 4472 | 30               F0=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:4475:72:

 4475 | 35               DV(K)=S0*DV(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:4536:72:

 4536 | 30            A1=AF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:4541:72:

 4541 | 35            BF=BF+A(2*K)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:4546:72:

 4546 | 40            BG=BG+A(2*K+1)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:4582:72:

 4582 | 20            ZP=ZP*(Z-ZO(I))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:4591:72:

 4591 | 30            ZQ=ZQ+ZW
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:4597:72:

 4597 | 35         ZO(NR)=Z
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:4622:72:

 4622 | 10               GA=GA*K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:4632:72:

 4632 | 15               R=R*(Z-K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:4651:72:

 4651 | 20            GR=GR*Z+G(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:4762:72:

 4762 | 10            DL(K)=0.0D0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:4778:72:

 4778 | 25            IF (K.GE.1) DL(K-1)=-0.5D0*X/K*BK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:4805:72:

 4805 | 40         F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:4808:72:

 4808 | 45         BL(K)=BL(K)/BG
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:4812:72:

 4812 | 50         BL(K)=R0*BL(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:4815:72:

 4815 | 55         DL(K)=2.0D0*K/X*(BL(K-1)-BL(K))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 55 at (1)
specfun.f:4873:72:

 4873 | 10            DK(2*K)=K*(B-K)*X/(A+2.0D0*K-1.0D0)/(A+2.0D0*K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:4875:72:

 4875 | 15            DK(2*K+1)=-(A+K)*(A+B+K)*X/(A+2.D0*K)/(A+2.0*K+1.0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:4878:72:

 4878 | 20            T1=DK(K)/(1.0D0+T1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:4883:72:

 4883 | 25            FK(2*K)=K*(A-K)*(1.0D0-X)/(B+2.*K-1.0)/(B+2.0*K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:4886:72:

 4886 |      &                   (B+2.D0*K)/(B+2.D0*K+1.D0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:4889:72:

 4889 | 35            T2=FK(K)/(1.0D0+T2)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:4924:72:

 4924 | 10         T1=-Q*Q/((2.0D0*J+L)**2-B+T1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:4939:72:

 4939 | 15            T2=-Q*Q/((2.0D0*J-L-L0)**2-B+T2)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:4981:72:

 4981 | 10         CP1=CPF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:5003:72:

 5003 | 10         QD(K)=0.0D0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:5007:72:

 5007 | 15            QD(K)=1.0D+300
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:5025:72:

 5025 | 20         QF1=QM0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:5037:72:

 5037 | 25               QF1=QF2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:5044:72:

 5044 | 30            QG1=QM1
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:5055:72:

 5055 | 35               QH1=QH2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:5071:72:

 5071 | 40                  QF1=QMK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:5076:72:

 5076 | 45               QH1=Q1L
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:5090:72:

 5090 | 50            QF1=QF0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:5092:72:

 5092 | 55            QM(K)=QM(K)*QM0/QF0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 55 at (1)
specfun.f:5096:72:

 5096 | 60            QM(K)=(-1)**M*QM(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:5100:72:

 5100 | 65         QD(K)=(K*X*QM(K)-(K+M)*QM(K-1))/(X*X-1.0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 65 at (1)
specfun.f:5140:72:

 5140 | 10               CF(K)=CF(K)*CT2+A(I)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:5141:72:

 5141 | 15            CF(K)=CF(K)*CT**K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:5145:72:

 5145 | 20            CSI=CSI+CF(K)*VR**K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:5150:72:

 5150 | 25            CSK=CSK+(-1)**K*CF(K)*VR**K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:5262:72:

 5262 | 10         EL3=EL3+W(I)*(F1+F2)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:5300:72:

 5300 | 25            EI=EI+R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:5440:72:

 5440 | 20               SUM2=SUM2+R2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:5460:72:

 5460 | 35            A=A+1.0D0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:5525:72:

 5525 | 10            HF=HF+R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:5534:72:

 5534 | 15            HF=HF+R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:5563:72:

 5563 | 30               GM=GM*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:5566:72:

 5566 | 35               RM=RM*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:5577:72:

 5577 | 40                  F0=F0+R0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:5579:72:

 5579 | 45                  SP0=SP0+1.0D0/(A+K-1.0)+1.0/(B+K-1.0)-1.0/K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:5586:72:

 5586 |      &                    (B+J+K-1.0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:5591:72:

 5591 | 55                  HW=F1
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 55 at (1)
specfun.f:5599:72:

 5599 | 65                  F0=F0+R0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 65 at (1)
specfun.f:5601:72:

 5601 | 70                  SP0=SP0+1.0D0/K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 70 at (1)
specfun.f:5607:72:

 5607 | 75                     SM=SM+1.0D0/(J+K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 75 at (1)
specfun.f:5612:72:

 5612 | 80                  HW=F1
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 80 at (1)
specfun.f:5634:72:

 5634 | 90               HW=HF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 90 at (1)
specfun.f:5651:72:

 5651 | 100           HW=HF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 100 at (1)
specfun.f:5709:72:

 5709 | 10            CHG=CHG+CR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:5755:72:

 5755 | 20                  CS2=CS2+CR2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:5781:72:

 5781 | 35               A=A+1.0D0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:5858:72:

 5858 | 10            ZHF=ZHF+ZR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:5866:72:

 5866 | 15            ZHF=ZHF+ZR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:5882:72:

 5882 | 20               ZW=ZHF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:5898:72:

 5898 | 30                  GM=GM*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:5901:72:

 5901 | 35                  RM=RM*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:5912:72:

 5912 | 40                     ZF0=ZF0+ZR0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:5914:72:

 5914 | 45                     SP0=SP0+1.0D0/(A+K-1.0D0)+1.0/(B+K-1.0D0)-1.D0/K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:5929:72:

 5929 | 55                     ZW=ZF1
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 55 at (1)
specfun.f:5938:72:

 5938 | 65                     ZF0=ZF0+ZR0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 65 at (1)
specfun.f:5940:72:

 5940 | 70                     SP0=SP0+1.0D0/K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 70 at (1)
specfun.f:5947:72:

 5947 | 75                        SM=SM+1.0D0/(J+K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 75 at (1)
specfun.f:5952:72:

 5952 | 80                     ZW=ZF1
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 80 at (1)
specfun.f:5974:72:

 5974 | 90                  ZW=ZHF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 90 at (1)
specfun.f:5990:72:

 5990 | 100              ZW=ZHF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 100 at (1)
specfun.f:6014:72:

 6014 | 110              ZW=ZHF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 110 at (1)
specfun.f:6042:72:

 6042 | 120              ZF0=ZF0+ZR*G0/GCBK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 120 at (1)
specfun.f:6047:72:

 6047 | 125              SP=SP+1.0D0/J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 125 at (1)
specfun.f:6051:72:

 6051 | 130              SQ=SQ*(B+J-1.0D0)*(B-C+J)/J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 130 at (1)
specfun.f:6062:72:

 6062 | 135                 RK2=RK2*(B+J-1.0D0)*(B-C+J)/J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 135 at (1)
specfun.f:6067:72:

 6067 | 140                 SJ2=SJ2+1.0D0/J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 140 at (1)
specfun.f:6073:72:

 6073 | 145              W0=WS
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 145 at (1)
specfun.f:6161:72:

 6161 | 35               SU1=SU1+A(K)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:6166:72:

 6166 | 40               SU2=SU2+A(K)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:6174:72:

 6174 | 45               SU3=SU3+A(2*K)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:6179:72:

 6179 | 50               SU4=SU4+A(2*K+1)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:6217:72:

 6217 | 10            DK(K)=-1.0D+300
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:6239:72:

 6239 | 15         H1=H
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:6254:72:

 6254 | 20            F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:6257:72:

 6257 | 25            BI(K)=S0*BI(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:6265:72:

 6265 | 30         G1=G
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:6268:72:

 6268 | 40         DK(K)=-BK(K-1)-K/X*BK(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:6309:72:

 6309 | 10            CDY(K)=(1.0D+300,0.0D0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:6341:72:

 6341 | 15            CF1=CF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:6348:72:

 6348 | 20            CBJ(K)=CBJ(K)/CS0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:6364:72:

 6364 | 25            CP0=CP0+A(K)*Z**(-2*K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:6367:72:

 6367 | 30            CQ0=CQ0+B(K)*Z**(-2*K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:6376:72:

 6376 | 35            CP1=CP1+A1(K)*Z**(-2*K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:6379:72:

 6379 | 40            CQ1=CQ1+B1(K)*Z**(-2*K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:6388:72:

 6388 | 45            CBJ1=CBJK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:6392:72:

 6392 | 50         CDJ(K)=CBJ(K-1)-K/Z*CBJ(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:6406:72:

 6406 | 60         CDY(K)=CBY(K-1)-K/Z*CBY(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:6441:72:

 6441 | 10            DK(K)=-1.0D+300
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:6464:72:

 6464 | 15         F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:6467:72:

 6467 | 20         BI(K)=S0*BI(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:6483:72:

 6483 | 25               BKL=BKL+R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:6493:72:

 6493 | 35         G1=G
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:6498:72:

 6498 | 40         DK(K)=-BK(K-1)-K/X*BK(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:6524:72:

 6524 |         DO 10 J=0,M
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 10 at (1)
specfun.f:6526:72:

 6526 | 10         PD(J,I)=0.0D0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:6532:72:

 6532 | 15            PD(0,I)=0.5D0*I*(I+1.0D0)*X**(I+1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:6534:72:

 6534 |            DO 20 I=1,M
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 20 at (1)
specfun.f:6550:72:

 6550 | 30         PM(I,I)=-LS*(2.0D0*I-1.0D0)*XQ*PM(I-1,I-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:6552:72:

 6552 | 35         PM(I,I+1)=(2.0D0*I+1.0D0)*X*PM(I,I)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:6554:72:

 6554 |         DO 40 J=I+2,N
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 40 at (1)
specfun.f:6560:72:

 6560 | 45         PD(0,J)=LS*J*(PM(0,J-1)-X*PM(0,J))/XS
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:6562:72:

 6562 |         DO 50 J=I,N
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 50 at (1)
specfun.f:6748:72:

 6748 | 50            CP0=CP0+A(K)*Z1**(-2*K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:6751:72:

 6751 | 55            CQ0=CQ0+B(K)*Z1**(-2*K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 55 at (1)
specfun.f:6758:72:

 6758 | 60            CP1=CP1+A1(K)*Z1**(-2*K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:6761:72:

 6761 | 65            CQ1=CQ1+B1(K)*Z1**(-2*K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 65 at (1)
specfun.f:6861:72:

 6861 | 25               XF0=XF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:6871:72:

 6871 | 30               XF=XF+XR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:6876:72:

 6876 | 35               XG=XG+XR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:7057:72:

 7057 | 50            DK(K)=-(6.0D0*K+1.0D0)/(6.0D0*K-1.0D0)*CK(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:7065:72:

 7065 | 55               SAD=SAD+DK(K)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 55 at (1)
specfun.f:7072:72:

 7072 | 60               SBD=SBD+DK(K)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:7088:72:

 7088 | 65               SDA=SDA+DK(2*K)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 65 at (1)
specfun.f:7095:72:

 7095 | 70               SDB=SDB+DK(2*K+1)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 70 at (1)
specfun.f:7145:72:

 7145 | 5                   CK(K1)=CK(K1)*1.0D-100
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 5 at (1)
specfun.f:7168:72:

 7168 | 10                  F2=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:7176:72:

 7176 | 25         SU1=SU1+CK(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:7179:72:

 7179 | 30         SU2=SU2+CK(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:7182:72:

 7182 | 35         R1=R1*(J+0.5D0*(N+M+IP))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:7185:72:

 7185 | 40         R2=-R2*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:7192:72:

 7192 | 45         CK(K)=FL/FS*S0*CK(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:7194:72:

 7194 | 50         CK(K)=S0*CK(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:7231:72:

 7231 | 10            R=R*I
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:7234:72:

 7234 | 15            R=R*(I+0.5D0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:7243:72:

 7243 | 20            SW=SUM
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:7246:72:

 7246 | 30            R1=R1*I
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:7247:72:

 7247 | 35         CK(K+1)=FAC*SUM/R1
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:7320:72:

 7320 | 20            ZP=ZP*(Z-ZO(I))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:7329:72:

 7329 | 30            ZQ=ZQ+ZW
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:7335:72:

 7335 | 35         ZO(NR)=Z
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:7397:72:

 7397 | 10            RG=RG*(V*V-J*J)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:7401:72:

 7401 | 15            R0=.5D0*R0*XQ/J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:7411:72:

 7411 | 20            PMV=PMV+R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:7432:72:

 7432 | 35                  R2=R2*QR*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:7437:72:

 7437 | 40                  S0=S0+R1
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:7444:72:

 7444 | 45               S1=S1+(J*J+V*V)/(J*(J*J-V*V))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:7451:72:

 7451 | 50                  S=S+((K+J)**2+V*V)/((K+J)*((K+J)**2-V*V))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:7454:72:

 7454 | 55                  S2=S2+1.0D0/(J*(J*J-V*V))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 55 at (1)
specfun.f:7579:72:

 7579 | 10         GI=GI-A(K)*T*DSIN((2.0D0*K-1.0D0)*TH)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:7585:72:

 7585 | 15            GI1=GI1+DATAN(Y/(X+J))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:7648:72:

 7648 | 10         SW=SU1
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:7655:72:

 7655 | 20         SW=SU1
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:7704:72:

 7704 | 10         H0=HU
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:7787:72:

 7787 | 10         GL0=GL0*X2+A(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:7792:72:

 7792 | 15            X0=X0-1.0D0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:7830:72:

 7830 | 15            Q1=QF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:7943:72:

 7943 | 35            BI0=BI0+A(K)*XR**K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:7947:72:

 7947 | 40            BI1=BI1+B(K)*XR**K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:7961:72:

 7961 | 65            WW=BK0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 65 at (1)
specfun.f:7972:72:

 7972 | 75            BK0=BK0+A1(K)*XR2**K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 75 at (1)
specfun.f:8017:72:

 8017 | 10            CF1=CF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:8035:72:

 8035 | 15               CF1=CF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:8048:72:

 8048 | 20                  CFA=CF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:8057:72:

 8057 | 25                  CFB=CF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:8060:72:

 8060 | 30                  CPB(K)=CS0*CPB(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:8067:72:

 8067 | 35            CPD(K)=-0.5D0*Z*CPB(K)+K*CPB(K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:8070:72:

 8070 | 40            CPD(K)=0.5D0*Z*CPB(K)-CPB(K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:8205:72:

 8205 | 10         P1=PF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:8238:72:

 8238 | 5          FC(I)=0.0D0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 5 at (1)
specfun.f:8252:72:

 8252 |  7666            FC(I)=FNAN
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 7666 at (1)
specfun.f:8296:72:

 8296 |  7            FC(I)=FNAN
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 7 at (1)
specfun.f:8348:72:

 8348 | 30            FC(K)=S0*FC(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:8367:72:

 8367 | 40            FC(K)=S0*FC(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:8410:72:

 8410 | 65            FC(K)=S0*FC(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 65 at (1)
specfun.f:8438:72:

 8438 | 90            FC(J)=-FC(J)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 90 at (1)
specfun.f:8468:72:

 8468 | 10            DI(K)=0.0D0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:8490:72:

 8490 | 15            F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:8493:72:

 8493 | 20            SI(K)=CS*SI(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:8497:72:

 8497 | 25         DI(K)=SI(K-1)-(K+1.0D0)/X*SI(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:8547:72:

 8547 | 10         H1=HL
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:8574:72:

 8574 | 40         D2=DL
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:8626:72:

 8626 | 10         R0=R0*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:8634:72:

 8634 | 15         SW=SUC
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:8643:72:

 8643 | 25            SW1=SUM
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:8646:72:

 8646 | 35            R1=R1*(J+0.5D0*(N+M+IP))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:8649:72:

 8649 | 40            R2=2.0D0*C*R2*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:8652:72:

 8652 | 45            R3=R3*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:8682:72:

 8682 | 50         SW=R1F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:8699:72:

 8699 | 60         SW=SUD
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:8922:72:

 8922 | 15            HU=HU+R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:8952:72:

 8952 | 10         GW=GF0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:8963:72:

 8963 | 20         GW=GD0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:9020:72:

 9020 | 30            A1=AF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:9025:72:

 9025 | 35            BF=BF+A(2*K)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:9030:72:

 9030 | 40            BG=BG+A(2*K+1)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:9059:72:

 9059 | 10            DY(K)=1.0D+300
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:9073:72:

 9073 | 15         RF1=RF2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:9077:72:

 9077 | 25         DY(K)=-K*RY(K)/X+RY(K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:9120:72:

 9120 | 10            R=(0.5D0/J-1.0D0)*R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:9246:72:

 9246 | 85            QN0=QN0+FAC*RS
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 85 at (1)
specfun.f:9343:72:

 9343 | 15         H0=HM1
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:9352:72:

 9352 | 25         IF (B.LT.0.0) S0=S0+(1.0D0-A)/(M*(A+M-1.0D0))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:9362:72:

 9362 | 30               S1=S1-(M+2.0D0*A-2.0D0)/(M*(M+A-1.0D0))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:9364:72:

 9364 | 35               S2=S2+1.0D0/(K+M)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:9367:72:

 9367 | 40               S1=S1+(1.0D0-A)/(M*(M+A-1.0D0))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:9369:72:

 9369 | 45               S2=S2+1.0D0/M
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:9378:72:

 9378 | 50         H0=HM2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:9389:72:

 9389 | 60         HM3=HM3+R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:9451:72:

 9451 | 20            ZP=ZP*(Z-ZO(I))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:9637:72:

 9637 | 10            SW=SUM
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:9690:72:

 9690 | 20         F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:9692:72:

 9692 | 25         BJ(K+1)=BJ(K+1)/(BS-F)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:9697:72:

 9697 | 30         FJ(K+1)=(K*K/(X*X)-1.0D0)*BJ(K+1)-DJ(K+1)/X
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:9726:72:

 9726 | 10            DJ(K)=0.0D0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:9755:72:

 9755 | 15            F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:9760:72:

 9760 | 20            SJ(K)=CS*SJ(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:9763:72:

 9763 | 25         DJ(K)=SJ(K-1)-(K+1.0D0)*SJ(K)/X
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:9826:72:

 9826 | 10         DY1=DYN
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:9887:72:

 9887 | 15         RT=RT+4.44D0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:9977:72:

 9977 | 10            CHD1(K)=CDJ(K)+CI*CDY(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:9984:72:

 9984 | 15            CFAC=CFAC*CI
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:9993:72:

 9993 | 20            CFAC=CFAC*CF1
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:9997:72:

 9997 | 25            CHD2(K)=CDJ(K)-CI*CDY(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:10004:72:

10004 | 30            CHD2(K)=CDJ(K)-CI*CDY(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:10176:72:

10176 | 10            DK(K)=1.0D+300
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:10211:72:

10211 | 25            SUM=SUM+R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:10228:72:

10228 | 30         F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:10231:72:

10231 | 35         BI(K)=CS*BI(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:10234:72:

10234 | 40         DI(K)=-(K+V0)/X*BI(K)+BI(K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:10247:72:

10247 | 45               WW=WA
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:10263:72:

10263 | 55               WW=WA
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 55 at (1)
specfun.f:10272:72:

10272 | 65            SUM=SUM+R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 65 at (1)
specfun.f:10282:72:

10282 | 70         BK1=BK2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 70 at (1)
specfun.f:10285:72:

10285 | 80         DK(K)=-(K+V0)/X*BK(K)-BK(K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 80 at (1)
specfun.f:10316:72:

10316 | 5             DF(I)=0D0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 5 at (1)
specfun.f:10350:72:

10350 | 12                  DF(K1)=DF(K1)*1.0D-100
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 12 at (1)
specfun.f:10372:72:

10372 | 15                        DF(K1)=DF(K1)*1.0D-100
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:10377:72:

10377 | 20                  F2=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:10386:72:

10386 | 40         R1=R1*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:10390:72:

10390 | 45           SU1=SU1+R1*DF(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:10397:72:

10397 | 50         SW=SU2
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:10400:72:

10400 | 60         R3=R3*(J+0.5D0*(N+M+IP))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:10403:72:

10403 | 65         R4=-4.0D0*R4*J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 65 at (1)
specfun.f:10406:72:

10406 | 70         DF(K)=FL/FS*S0*DF(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 70 at (1)
specfun.f:10408:72:

10408 | 75         DF(K)=S0*DF(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 75 at (1)
specfun.f:10481:72:

10481 | 30               PX=PX+RP
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:10487:72:

10487 | 35               QX=QX+RQ
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:10545:72:

10545 | 75               VSL=VSL+R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 75 at (1)
specfun.f:10575:72:

10575 | 100              SUM=SUM+R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 100 at (1)
specfun.f:10626:72:

10626 | 10            CDK(K)=1.0D+300
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:10660:72:

10660 | 25            CS=CS+CR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:10675:72:

10675 | 30         CF1=CF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:10678:72:

10678 | 35         CBI(K)=CS*CBI(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:10716:72:

10716 | 60            CS=CS+CR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:10740:72:

10740 | 80         CDK(K)=-(K+V0)/Z*CBK(K)-CBK(K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 80 at (1)
specfun.f:10788:72:

10788 | 10            CDK(K)=1.0D+300
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:10822:72:

10822 | 25            CS=CS+CR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:10837:72:

10837 | 30         CF1=CF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:10840:72:

10840 | 35         CBI(K)=CS*CBI(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:10880:72:

10880 | 60            CS=CS+CR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:10892:72:

10892 | 65         CG1=CGK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 65 at (1)
specfun.f:10909:72:

10909 | 75         CDK(K)=-(K+V0)/Z*CBK(K)-CBK(K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 75 at (1)
specfun.f:10947:72:

10947 | 10            WA0=WA
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:10957:72:

10957 | 15            CF0=CF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:10978:72:

10978 | 20            CF=CF+CR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:10983:72:

10983 | 25            CG=CG+CR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:11047:72:

11047 | 25            F0=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:11056:72:

11056 | 30            F=F+R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:11061:72:

11061 | 35            G=G+R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:11096:72:

11096 | 10            DJ(K)=0.0D0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:11119:72:

11119 | 15            F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:11123:72:

11123 | 20            RJ(K)=CS*RJ(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:11127:72:

11127 | 25         DJ(K)=-K*RJ(K)/X+RJ(K-1)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:11166:72:

11166 | 15            F1=HF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:11169:72:

11169 | 20            P=P*(Z-X(I))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:11178:72:

11178 | 30            Q=Q+WP
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:11186:72:

11186 | 35            R=2.0D0*R*K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:11188:72:

11188 | 40         W(N+1-NR)=W(NR)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:11298:72:

11298 | 10            EN(K)=1.0D0/(K-1.0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:11306:72:

11306 | 15               RP=-RP*X/J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:11309:72:

11309 | 20               PS=PS+1.0D0/M
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:11316:72:

11316 | 25                  R=-R*X/J
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:11329:72:

11329 | 45               T0=(L+K-1.0D0)/(1.0D0+K/(X+T0))
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:11331:72:

11331 | 50            EN(L)=EXP(-X)*T
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:11356:72:

11356 | 10            DK(K)=-1.0D+300
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:11368:72:

11368 | 15         F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:11372:72:

11372 | 25         DK(K)=-SK(K-1)-(K+1.0D0)/X*SK(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:11396:72:

11396 | 10         E1=EK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:11418:72:

11418 | 10            GA=GA*K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:11423:72:

11423 | 15            GA=0.5D0*GA*(2.0D0*K-1.0D0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
mvn.f:796:72:

  796 |       END
      |                                                                        ^
Warning: ‘__result_bvn’ may be used uninitialized [-Wmaybe-uninitialized]
mvn.f:754:35:

  754 |       DOUBLE PRECISION FUNCTION BVN ( LOWER, UPPER, INFIN, CORREL )
      |                                   ^
note: ‘__result_bvn’ was declared here
specfun.f:11487:72:

11487 | 15            F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:11502:72:

11502 | 20               F0=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:11513:72:

11513 | 25               F0=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:11516:72:

11516 | 30               VV(K)=S0*VV(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:11533:72:

11533 | 35               F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:11576:72:

11576 |            DO 10 J=0,N
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 10 at (1)
specfun.f:11595:72:

11595 |            DO 15 J=2,N
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 15 at (1)
specfun.f:11600:72:

11600 |            DO 20 I=2,M
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 20 at (1)
specfun.f:11616:72:

11616 | 25            CQF1=CQF0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:11618:72:

11618 | 30            CQM(0,K)=CQ0*CQM(0,K)/CQF0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:11625:72:

11625 | 35            CQF1=CQF0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:11628:72:

11628 | 40            CQM(1,K)=CQ10*CQM(1,K)/CQF0
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:11632:72:

11632 |               DO 45 I=0,M-2
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 45 at (1)
specfun.f:11641:72:

11641 | 50         CQD(0,J)=LS*J*(CQM(0,J-1)-Z*CQM(0,J))/ZS
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:11643:72:

11643 |         DO 55 I=1,M
      |                                                                        1
Warning: Fortran 2018 deleted feature: Shared DO termination label 55 at (1)
specfun.f:11673:72:

11673 | 5             EG(I)=(I+M)*(I+M-1.0D0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 5 at (1)
specfun.f:11691:72:

11691 | 10            G(I)=K*(K-1.0)/((DK2-3.0)*(DK2-1.0))*CS
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:11694:72:

11694 | 15            F(K)=E(K)*E(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:11709:72:

11709 | 25            H(I)=XB
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:11785:72:

11785 | 10            CDK(K)=-(1.0D+300,0.0D0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:11810:72:

11810 | 15         CF1=CF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:11813:72:

11813 | 20         CBI(K)=CS0*CBI(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:11829:72:

11829 | 25               CBKL=CBKL+CR
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:11839:72:

11839 | 35         CG1=CG
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:11856:72:

11856 | 50         CDK(K)=-CBK(K-1)-K/Z*CBK(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:11892:72:

11892 | 10            CDK(K)=(1.0D+300,0.0D0)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:11919:72:

11919 | 45         CF1=CF
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 45 at (1)
specfun.f:11922:72:

11922 | 50         CBI(K)=CS*CBI(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:11929:72:

11929 | 60         CBK(K)=CKK
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 60 at (1)
specfun.f:11932:72:

11932 | 70         CDK(K)=-CBK(K-1)-K/Z*CBK(K)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 70 at (1)
specfun.f:12015:72:

12015 | 30         W1=F1R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:12032:72:

12032 | 40         W2=D1R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:12047:72:

12047 | 55         W1=F2R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 55 at (1)
specfun.f:12064:72:

12064 | 65         W2=D2R
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 65 at (1)
specfun.f:12147:72:

12147 | 30            CBI0=CBI0+A(K)*ZR**K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 30 at (1)
specfun.f:12151:72:

12151 | 35            CBI1=CBI1+B(K)*ZR**K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 35 at (1)
specfun.f:12164:72:

12164 | 40            CW=CS
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 40 at (1)
specfun.f:12176:72:

12176 | 50            CBK0=CBK0+A1(K)*ZR2**K
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 50 at (1)
specfun.f:12239:72:

12239 | 10            PSI=PSI-A(K)*Z2**(-K)*DSIN(2.0D0*K*TH)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:12245:72:

12245 | 20               RI=RI+Y/((X0-K)**2.0D0+Y*Y)
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:12282:72:

12282 | 10            DY(K)=1.0D+300
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:12298:72:

12298 | 15         F1=F
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 15 at (1)
specfun.f:12301:72:

12301 | 25         DY(K)=SY(K-1)-(K+1.0D0)*SY(K)/X
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 25 at (1)
specfun.f:12332:72:

12332 | 10         B0=B
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 10 at (1)
specfun.f:12339:72:

12339 | 20         DN=D
      |                                                                        1
Warning: Fortran 2018 deleted feature: DO termination statement which is not END DO or CONTINUE with label 20 at (1)
specfun.f:1276:72:

 1276 |         CDJV=-V/Z*CBJV+CFJ
      |                                                                        ^
Warning: ‘IMAGPART_EXPR <cfj>’ may be used uninitialized [-Wmaybe-uninitialized]
specfun.f:1269:26:

 1269 |            IF (L.EQ.1) CFJ=CBJV
      |                          ^
note: ‘IMAGPART_EXPR <cfj>’ was declared here
specfun.f:1276:72:

 1276 |         CDJV=-V/Z*CBJV+CFJ
      |                                                                        ^
Warning: ‘REALPART_EXPR <cfj>’ may be used uninitialized [-Wmaybe-uninitialized]
specfun.f:1269:26:

 1269 |            IF (L.EQ.1) CFJ=CBJV
      |                          ^
note: ‘REALPART_EXPR <cfj>’ was declared here
specfun.f:5154:72:

 5154 |         CDIV=CFI-V/Z*CBIV
      |                                                                        ^
Warning: ‘IMAGPART_EXPR <cfi>’ may be used uninitialized [-Wmaybe-uninitialized]
specfun.f:5147:26:

 5147 |            IF (L.EQ.1) CFI=CBIV
      |                          ^
note: ‘IMAGPART_EXPR <cfi>’ was declared here
specfun.f:5154:72:

 5154 |         CDIV=CFI-V/Z*CBIV
      |                                                                        ^
Warning: ‘REALPART_EXPR <cfi>’ may be used uninitialized [-Wmaybe-uninitialized]
specfun.f:5147:26:

 5147 |            IF (L.EQ.1) CFI=CBIV
      |                          ^
note: ‘REALPART_EXPR <cfi>’ was declared here
clang -arch x86_64 -dynamiclib -Wl,-headerpad_max_install_names -undefined dynamic_lookup -L/Library/Frameworks/R.framework/Resources/lib -L/opt/R/x86_64/lib -o GeoModels.so CompositeLikelihood2.o CompositeLikelihood2_ani.o CompositeLikelihoodCond2.o CompositeLikelihoodCond2_ani.o CorrelationFunction.o Distributions.o GeoModels_init.o Godambe.o Gradient.o HypergeometricFunctions.o TB.o Utility.o mvn.o mvndstpack.o specfun.o weightedleastsquare.o -L/Library/Frameworks/R.framework/Resources/lib -lRlapack -L/Library/Frameworks/R.framework/Resources/lib -lRblas -L/opt/gfortran/lib/gcc/x86_64-apple-darwin20.0/12.2.0 -L/opt/gfortran/lib -lgfortran -lquadmath -L/opt/gfortran/lib/gcc/x86_64-apple-darwin20.0/12.2.0 -L/opt/gfortran/lib -lgfortran -lquadmath -F/Library/Frameworks/R.framework/.. -framework R -Wl,-framework -Wl,CoreFoundation
ld: warning: could not create compact unwind for _dkbvrc_: stack subq instruction is too different from dwarf stack size
ld: warning: could not create compact unwind for _sadmvn_: stack subq instruction is too different from dwarf stack size
ld: warning: could not create compact unwind for _master.0.mvndfn_: stack subq instruction is too different from dwarf stack size
ld: warning: could not create compact unwind for _sadmvt_: stack subq instruction is too different from dwarf stack size
ld: warning: could not create compact unwind for _rmn2so_.part.0: stack subq instruction is too different from dwarf stack size
ld: warning: could not create compact unwind for _rmn2sp_: stack subq instruction is too different from dwarf stack size
ld: warning: could not create compact unwind for _ch12n_: stack subq instruction is too different from dwarf stack size
ld: warning: could not create compact unwind for _mtu12_: stack subq instruction is too different from dwarf stack size
installing to /Volumes/Builds/packages/big-sur-x86_64/results/4.3/GeoModels.Rcheck/00LOCK-GeoModels/00new/GeoModels/libs
** generating debug symbols (dSYM)
** R
** data
** byte-compile and prepare package for lazy loading
** help
*** installing help indices
** building package indices
** testing if installed package can be loaded from temporary location
** checking absolute paths in shared objects and dynamic libraries
** testing if installed package can be loaded from final location
** testing if installed package keeps a record of temporary installation path
* DONE (GeoModels)