* installing to library ‘/home/hornik/tmp/R.check/r-patched-gcc/Work/build/Packages’
* installing *source* package ‘GeoModels’ ...
** package ‘GeoModels’ successfully unpacked and MD5 sums checked
** using staged installation
** libs
using C compiler: ‘gcc-13 (Debian 13.2.0-24) 13.2.0’
using Fortran compiler: ‘GNU Fortran (Debian 13.2.0-24) 13.2.0’
make[1]: Entering directory '/tmp/Rtmp1mOIgS/R.INSTALL17ff5b48950a5a/GeoModels/src'
gcc-13 -I"/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include" -DNDEBUG -I/usr/local/include -D_FORTIFY_SOURCE=3 -fpic -g -O2 -Wall -Wstrict-prototypes -pedantic -mtune=native -c CompositeLikelihood2.c -o CompositeLikelihood2.o
gcc-13 -I"/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include" -DNDEBUG -I/usr/local/include -D_FORTIFY_SOURCE=3 -fpic -g -O2 -Wall -Wstrict-prototypes -pedantic -mtune=native -c CompositeLikelihood2_ani.c -o CompositeLikelihood2_ani.o
gcc-13 -I"/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include" -DNDEBUG -I/usr/local/include -D_FORTIFY_SOURCE=3 -fpic -g -O2 -Wall -Wstrict-prototypes -pedantic -mtune=native -c CompositeLikelihoodCond2.c -o CompositeLikelihoodCond2.o
gcc-13 -I"/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include" -DNDEBUG -I/usr/local/include -D_FORTIFY_SOURCE=3 -fpic -g -O2 -Wall -Wstrict-prototypes -pedantic -mtune=native -c CompositeLikelihoodCond2_ani.c -o CompositeLikelihoodCond2_ani.o
gcc-13 -I"/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include" -DNDEBUG -I/usr/local/include -D_FORTIFY_SOURCE=3 -fpic -g -O2 -Wall -Wstrict-prototypes -pedantic -mtune=native -c CorrelationFunction.c -o CorrelationFunction.o
gcc-13 -I"/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include" -DNDEBUG -I/usr/local/include -D_FORTIFY_SOURCE=3 -fpic -g -O2 -Wall -Wstrict-prototypes -pedantic -mtune=native -c Distributions.c -o Distributions.o
gcc-13 -I"/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include" -DNDEBUG -I/usr/local/include -D_FORTIFY_SOURCE=3 -fpic -g -O2 -Wall -Wstrict-prototypes -pedantic -mtune=native -c GeoModels_init.c -o GeoModels_init.o
gcc-13 -I"/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include" -DNDEBUG -I/usr/local/include -D_FORTIFY_SOURCE=3 -fpic -g -O2 -Wall -Wstrict-prototypes -pedantic -mtune=native -c Godambe.c -o Godambe.o
In file included from /home/hornik/tmp/R.check/r-patched-gcc/Work/build/include/R_ext/BLAS.h:38,
from header.h:7,
from Godambe.c:1:
Godambe.c: In function ‘Vari_SubSamp’:
/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include/R_ext/RS.h:55:25: warning: ‘ecoordx’ may be used uninitialized [-Wmaybe-uninitialized]
55 | #define Free(p) (R_chk_free( (void *)(p) ), (p) = NULL)
| ^~~~~~~~~~~~~~~~~~~~~~~~~
Godambe.c:490:15: note: in expansion of macro ‘Free’
490 | if(*grid){Free(ecoordx);Free(ecoordy);}
| ^~~~
Godambe.c:181:13: note: ‘ecoordx’ was declared here
181 | double *ecoordx,*ecoordy,*scoordx,*scoordy,*sdata,*sumgrad,*subvari,*xgrid,*ygrid;
| ^~~~~~~
/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include/R_ext/RS.h:55:25: warning: ‘ecoordy’ may be used uninitialized [-Wmaybe-uninitialized]
55 | #define Free(p) (R_chk_free( (void *)(p) ), (p) = NULL)
| ^~~~~~~~~~~~~~~~~~~~~~~~~
Godambe.c:490:29: note: in expansion of macro ‘Free’
490 | if(*grid){Free(ecoordx);Free(ecoordy);}
| ^~~~
Godambe.c:181:22: note: ‘ecoordy’ was declared here
181 | double *ecoordx,*ecoordy,*scoordx,*scoordy,*sdata,*sumgrad,*subvari,*xgrid,*ygrid;
| ^~~~~~~
Godambe.c: In function ‘Vari_SubSamp_biv’:
/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include/R_ext/RS.h:55:25: warning: ‘ecoordx’ may be used uninitialized [-Wmaybe-uninitialized]
55 | #define Free(p) (R_chk_free( (void *)(p) ), (p) = NULL)
| ^~~~~~~~~~~~~~~~~~~~~~~~~
Godambe.c:1592:9: note: in expansion of macro ‘Free’
1592 | Free(ecoordx);
| ^~~~
Godambe.c:1140:58: note: ‘ecoordx’ was declared here
1140 | double *gradcortt,*gradcortv,*gradcorvt,*gradcorvv,*ecoordx,*ecoordy,*scoordx,*scoordy,*sdata,*score,*subvari,*xgrid,*ygrid;
| ^~~~~~~
/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include/R_ext/RS.h:55:25: warning: ‘ecoordy’ may be used uninitialized [-Wmaybe-uninitialized]
55 | #define Free(p) (R_chk_free( (void *)(p) ), (p) = NULL)
| ^~~~~~~~~~~~~~~~~~~~~~~~~
Godambe.c:1593:9: note: in expansion of macro ‘Free’
1593 | Free(ecoordy);
| ^~~~
Godambe.c:1140:67: note: ‘ecoordy’ was declared here
1140 | double *gradcortt,*gradcortv,*gradcorvt,*gradcorvv,*ecoordx,*ecoordy,*scoordx,*scoordy,*sdata,*score,*subvari,*xgrid,*ygrid;
| ^~~~~~~
gcc-13 -I"/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include" -DNDEBUG -I/usr/local/include -D_FORTIFY_SOURCE=3 -fpic -g -O2 -Wall -Wstrict-prototypes -pedantic -mtune=native -c Gradient.c -o Gradient.o
gfortran-13 -fpic -g -O2 -Wall -mtune=native -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]
gcc-13 -I"/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include" -DNDEBUG -I/usr/local/include -D_FORTIFY_SOURCE=3 -fpic -g -O2 -Wall -Wstrict-prototypes -pedantic -mtune=native -c TB.c -o TB.o
gcc-13 -I"/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include" -DNDEBUG -I/usr/local/include -D_FORTIFY_SOURCE=3 -fpic -g -O2 -Wall -Wstrict-prototypes -pedantic -mtune=native -c Utility.c -o Utility.o
gfortran-13 -fpic -g -O2 -Wall -mtune=native -c mvn.f -o mvn.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]
mvn.f:1154:72:
1154 | G(I) = G(IXCHNG)
| ^
Warning: ‘lxchng’ may be used uninitialized [-Wmaybe-uninitialized]
mvn.f:1121:31:
1121 | INTEGER S, IXCHNG, LXCHNG, I, L
| ^
note: ‘lxchng’ was declared here
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
gfortran-13 -fpic -g -O2 -Wall -mtune=native -c mvndstpack.f -o mvndstpack.o
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
gfortran-13 -fpic -g -O2 -Wall -mtune=native -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)
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)
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)
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:569:72:
569 | W(1)=W(1)/V(1)
| ^
Warning: ‘v’ may be used uninitialized [-Wmaybe-uninitialized]
specfun.f:522:42:
522 | DIMENSION BK(200),CK(200),U(200),V(200),W(200)
| ^
note: ‘v’ declared here
specfun.f:4463:72:
4463 | DV(0)=PD0
| ^
Warning: ‘pd0’ may be used uninitialized [-Wmaybe-uninitialized]
specfun.f:4411:17:
4411 | PD0=EP
| ^
note: ‘pd0’ declared here
specfun.f:11468:72:
11468 | VV(0)=PV0
| ^
Warning: ‘pv0’ may be used uninitialized [-Wmaybe-uninitialized]
specfun.f:11464:47:
11464 | IF (XA.LE.7.5) CALL VVSA(V0,X,PV0)
| ^
note: ‘pv0’ declared here
gcc-13 -I"/home/hornik/tmp/R.check/r-patched-gcc/Work/build/include" -DNDEBUG -I/usr/local/include -D_FORTIFY_SOURCE=3 -fpic -g -O2 -Wall -Wstrict-prototypes -pedantic -mtune=native -c weightedleastsquare.c -o weightedleastsquare.o
gcc-13 -shared -L/home/hornik/tmp/R.check/r-patched-gcc/Work/build/lib -Wl,-O1 -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 -llapack -L/home/hornik/tmp/R.check/r-patched-gcc/Work/build/lib -lRblas -lgfortran -lm -lquadmath -lgfortran -lm -lquadmath -L/home/hornik/tmp/R.check/r-patched-gcc/Work/build/lib -lR
make[1]: Leaving directory '/tmp/Rtmp1mOIgS/R.INSTALL17ff5b48950a5a/GeoModels/src'
make[1]: Entering directory '/tmp/Rtmp1mOIgS/R.INSTALL17ff5b48950a5a/GeoModels/src'
make[1]: Leaving directory '/tmp/Rtmp1mOIgS/R.INSTALL17ff5b48950a5a/GeoModels/src'
installing to /home/hornik/tmp/R.check/r-patched-gcc/Work/build/Packages/00LOCK-GeoModels/00new/GeoModels/libs
** R
** data
** byte-compile and prepare package for lazy loading
** help
*** installing help indices
converting help for package ‘GeoModels’
finding HTML links ... done
Anomalies html
CheckBiv html
CheckDistance html
CheckST html
CheckSph html
CkCorrModel html
CkInput html
CkLikelihood html
CkModel html
CkType html
CkVarType html
CompIndLik2 html
CompLik html
CompLik2 html
CorrParam html
CorrelationPar html
GeoAniso html
GeoCV html
GeoCorrFct html
GeoCorrFct_Cop html
GeoCovDisplay html
GeoCovariogram html
GeoCovmatrix html
GeoFit html
GeoFit2 html
GeoKrig html
GeoKrigloc html
GeoNA html
GeoNeighIndex html
GeoNeighborhood html
GeoNosymindices html
GeoOutlier html
GeoPit html
GeoQQ html
GeoResiduals html
GeoScatterplot html
GeoSim html
GeoSimCopula html
GeoSimapprox html
GeoTests html
GeoVarestbootstrap html
GeoVariogram html
GeoWls html
Lik html
MatDecomp html
Matrixfun html
NuisParam html
NuisParam2 html
Prscores html
SimCE html
StartParam html
WlsStart html
austemp html
plot.GeoVariogram html
sp2Geo html
spanish_wind html
winds html
winds.coords html
** 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
* creating tarball
packaged installation of ‘GeoModels’ as ‘GeoModels_2.0.1_R_x86_64-pc-linux-gnu.tar.gz’
* DONE (GeoModels)