* installing *source* package 'GeoModels' ... ** using staged installation ** libs using C compiler: 'gcc.exe (GCC) 13.2.0' using Fortran compiler: 'GNU Fortran (GCC) 13.2.0' gcc -I"D:/RCompile/recent/R/include" -DNDEBUG -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include" -pedantic -Wstrict-prototypes -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c CompositeLikelihood2.c -o CompositeLikelihood2.o gcc -I"D:/RCompile/recent/R/include" -DNDEBUG -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include" -pedantic -Wstrict-prototypes -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c CompositeLikelihood2_ani.c -o CompositeLikelihood2_ani.o gcc -I"D:/RCompile/recent/R/include" -DNDEBUG -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include" -pedantic -Wstrict-prototypes -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c CompositeLikelihoodCond2.c -o CompositeLikelihoodCond2.o gcc -I"D:/RCompile/recent/R/include" -DNDEBUG -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include" -pedantic -Wstrict-prototypes -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c CompositeLikelihoodCond2_ani.c -o CompositeLikelihoodCond2_ani.o gcc -I"D:/RCompile/recent/R/include" -DNDEBUG -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include" -pedantic -Wstrict-prototypes -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c CorrelationFunction.c -o CorrelationFunction.o gcc -I"D:/RCompile/recent/R/include" -DNDEBUG -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include" -pedantic -Wstrict-prototypes -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c Distributions.c -o Distributions.o Distributions.c: In function 'corr_tukeygh': Distributions.c:1689:54: warning: 'rho1' may be used uninitialized [-Wmaybe-uninitialized] 1689 | double mu,rho2,a,eta2,tail2,u,A1,A2,A3,cova,vari,rho1; | ^~~~ Distributions.c: In function 'biv_two_piece_bimodal': Distributions.c:1521:22: warning: 'res' may be used uninitialized [-Wmaybe-uninitialized] 1521 | return(R_pow(alpha,2)*res/sill); | ~~~~~~~~~~~~~~^~~~ Distributions.c:1504:8: note: 'res' was declared here 1504 | double res; | ^~~ Distributions.c: In function 'biv_two_pieceGaussian': Distributions.c:3399:11: warning: 'res' may be used uninitialized [-Wmaybe-uninitialized] 3399 | return(res/sill); | ~~~~^~~~~~ Distributions.c:3379:8: note: 'res' was declared here 3379 | double res; | ^~~ Distributions.c: In function 'biv_two_pieceTukeyh': Distributions.c:3758:11: warning: 'res' may be used uninitialized [-Wmaybe-uninitialized] 3758 | return(res/sill); | ~~~~^~~~~~ Distributions.c:3733:8: note: 'res' was declared here 3733 | double res; | ^~~ Distributions.c: In function 'one_log_tukeyhh': Distributions.c:4863:9: warning: 'res' may be used uninitialized [-Wmaybe-uninitialized] 4863 | return(res); | ^ Distributions.c:4856:11: note: 'res' was declared here 4856 | double res; | ^~~ Distributions.c: In function 'one_log_two_pieceTukey': Distributions.c:4943:13: warning: 'res' may be used uninitialized [-Wmaybe-uninitialized] 4943 | return(res-log(sqrt(sill))); | ~~~~^~~~~~~~~~~~~~~~~ Distributions.c:4939:11: note: 'res' was declared here 4939 | double res; | ^~~ Distributions.c: In function 'one_log_two_pieceT': Distributions.c:4951:13: warning: 'res' may be used uninitialized [-Wmaybe-uninitialized] 4951 | return(res-log(sqrt(sill))); | ~~~~^~~~~~~~~~~~~~~~~ Distributions.c:4947:11: note: 'res' was declared here 4947 | double res; | ^~~ Distributions.c: In function 'one_log_two_pieceGauss': Distributions.c:4959:13: warning: 'res' may be used uninitialized [-Wmaybe-uninitialized] 4959 | return(res-log(sqrt(sill))); | ~~~~^~~~~~~~~~~~~~~~~ Distributions.c:4955:11: note: 'res' was declared here 4955 | double res; | ^~~ Distributions.c: In function 'one_log_bomidal': Distributions.c:4984:9: warning: 'res' may be used uninitialized [-Wmaybe-uninitialized] 4984 | return(res); | ^ Distributions.c:4974:11: note: 'res' was declared here 4974 | double res; | ^~~ gcc -I"D:/RCompile/recent/R/include" -DNDEBUG -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include" -pedantic -Wstrict-prototypes -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c GeoModels_init.c -o GeoModels_init.o gcc -I"D:/RCompile/recent/R/include" -DNDEBUG -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include" -pedantic -Wstrict-prototypes -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c Godambe.c -o Godambe.o In file included from D:/RCompile/recent/R/include/R.h:78, from header.h:8, from Godambe.c:1: Godambe.c: In function 'Vari_SubSamp': D:/RCompile/recent/R/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; | ^~~~~~~ D:/RCompile/recent/R/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': D:/RCompile/recent/R/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; | ^~~~~~~ D:/RCompile/recent/R/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 -I"D:/RCompile/recent/R/include" -DNDEBUG -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include" -pedantic -Wstrict-prototypes -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c Gradient.c -o Gradient.o gfortran -pedantic -O2 -mfpmath=sse -msse2 -mstackrealign -c HypergeometricFunctions.f -o HypergeometricFunctions.o HypergeometricFunctions.f:43:21: 43 | DOUBLE COMPLEX CHGF,A,B,Z,CONHYP | 1 Warning: GNU Extension: DOUBLE COMPLEX at (1) HypergeometricFunctions.f:168:21: 168 | DOUBLE COMPLEX A,B,Z,FINAL,CHGF | 1 Warning: GNU Extension: DOUBLE COMPLEX at (1) HypergeometricFunctions.f:675:21: 675 | DOUBLE COMPLEX C | 1 Warning: GNU Extension: DOUBLE COMPLEX at (1) HypergeometricFunctions.f:863:21: 863 | DOUBLE COMPLEX CN | 1 Warning: GNU Extension: DOUBLE COMPLEX at (1) HypergeometricFunctions.f:908:21: 908 | DOUBLE COMPLEX CN | 1 Warning: GNU Extension: DOUBLE COMPLEX at (1) HypergeometricFunctions.f:987:21: 987 | DOUBLE COMPLEX Z, A, B, CHF, CONHYP | 1 Warning: GNU Extension: DOUBLE COMPLEX at (1) gcc -I"D:/RCompile/recent/R/include" -DNDEBUG -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include" -pedantic -Wstrict-prototypes -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c TB.c -o TB.o TB.c: In function 'for_c': TB.c:373:9: warning: unused variable 'model_1' [-Wunused-variable] 373 | int model_1 = 14; | ^~~~~~~ gcc -I"D:/RCompile/recent/R/include" -DNDEBUG -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include" -pedantic -Wstrict-prototypes -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c Utility.c -o Utility.o gfortran -pedantic -O2 -mfpmath=sse -msse2 -mstackrealign -c mvn.f -o mvn.o gfortran -pedantic -O2 -mfpmath=sse -msse2 -mstackrealign -c mvndstpack.f -o mvndstpack.o mvndstpack.f:559:72: 559 | & 8*12, 7, 3*3, 3*7, 21*3/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:563:72: 563 | & 2*10, 8*15, 6, 2, 3, 2, 3, 12*2/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:566:72: 566 | & 16*14, 2*31, 3*5, 11, 13, 7*11, 2*13, 11, 13, 4*5, 14, 13, 8*5/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:569:72: 569 | & 17, 3*29, 17, 4*5, 16*21, 3*17, 6, 2*17, 6, 3, 2*6, 5*3/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:572:72: 572 | & 6*45, 19*66, 11, 9*66, 45, 11, 7, 3, 3*2, 27, 5, 2*3, 2*5, 7*2/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:576:72: 576 | & 51, 12, 51, 5, 2*3, 2*2, 5/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:579:72: 579 | & 6*116, 5*100, 5*138, 19*101, 8*38, 5*3/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:582:72: 582 | & 12*49, 34*171, 8*161, 17*14, 6*10, 103, 4*10, 5/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:585:72: 585 | & 2*231, 14*90, 4*48, 23*90, 10*243, 9*283, 16, 283, 16, 2*283/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:589:72: 589 | & 8*101, 7*193, 21*101, 17*122, 4*101/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:593:72: 593 | & 7*195, 11*132, 13*387/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:597:72: 597 | & 838, 3*563, 2*759, 564, 2*759, 4*801, 5*759, 8*563, 22*226/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:602:72: 602 | & 4*301, 2*86, 15/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:605:72: 605 | & 2*3117, 5*1101, 8*2503, 7*429, 3*1702, 5*184, 34*105, 13*784/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:609:72: 609 | & 13*2339, 18*3148, 3*1776, 3*3354, 925, 2*3354, 5*925, 8*2133/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:615:72: 615 | & 4*1097, 4*222, 754,11*1982/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:620:72: 620 | & 24*8204, 3*4688, 8*2831/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:625:72: 625 | & 4377, 7096, 5*4377, 2*5410, 32*4377, 2*440, 3*1199/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:631:72: 631 | & 6*4544, 4*8394/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:636:72: 636 | & 5*6545, 14*12138, 5*30483, 19*12138, 9305, 13*11107, 2*9305/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:642:72: 642 | & 13*26401, 12982, 6*40398, 3*3518, 9*37799, 4*4721, 4*7067/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:648:72: 648 | & 9*62060, 13*1803, 12*51108, 2*55315, 5*54140, 13134/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:655:72: 655 | & 7*40586, 4*38725, 5*88329/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:662:72: 662 | & 6*80822, 16*131661, 7114, 131661/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:669:72: 669 | & 7*245397, 94241, 2*66575, 19*217673, 10*210249, 15*94453/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:677:72: 677 | & 5*13942/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:684:72: 684 | & 19*542095, 3*277743, 12*457259/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement mvndstpack.f:692:72: 692 | & 4*178309, 8*74373, 3*214965/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement gfortran -pedantic -O2 -mfpmath=sse -msse2 -mstackrealign -c specfun.f -o specfun.o specfun.f:49:28: 49 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 at (1) specfun.f:102:29: 102 | IMPLICIT COMPLEX *16 (C,S,Z,D) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 at (1) 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:281:28: 281 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:1245:28: 1245 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:1582:20: 1582 | COMPLEX *16 Z,CER,CDER | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:2007:28: 2007 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:2216:29: 2216 | IMPLICIT COMPLEX *16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:2601:72: 2601 | & 1.7588273098916D+4,1.4950639538279D+5/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement 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:3099:28: 3099 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:3363:28: 3363 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:3584:72: 3584 | & -.3833534661393944D+13,.4854014686852901D+15/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement specfun.f:3590:72: 3590 | & .4218971570284096D+14,-.5827244631566907D+16/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement specfun.f:3596:72: 3596 | & .4011838599133198D+13,-.5060568503314727D+15/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement specfun.f:3602:72: 3602 | & -.4406481417852278D+14,.6065091351222699D+16/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement 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:3771:72: 3771 | & 9.515939374212D+4,9.0412425769041D+5/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement 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:4125:72: 4125 | & -.5725014209747314D+00,.6074042001273483D+01/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement specfun.f:4127:72: 4127 | & .1727727502584457D+01,-.2438052969955606D+02/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement specfun.f:4129:72: 4129 | & .6765925884246826D+00,-.6883914268109947D+01/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement specfun.f:4131:72: 4131 | & -.1993531733751297D+01,.2724882731126854D+02/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement 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:4567:29: 4567 | IMPLICIT COMPLEX *16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:4648:72: 4648 | & -.206D-13, -.54D-14, .14D-14, .1D-15/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement 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:4963:29: 4963 | IMPLICIT COMPLEX *16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:5123:28: 5123 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:5317:23: 5317 | DOUBLE COMPLEX Z, CEI | 1 Warning: GNU Extension: DOUBLE COMPLEX at (1) specfun.f:5384:28: 5384 | IMPLICIT COMPLEX*16 (C) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:5683:29: 5683 | IMPLICIT COMPLEX *16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:5811:29: 5811 | IMPLICIT COMPLEX *16 (Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:6150:72: 6150 | & .231369166433050D+12,.358622522796969D+13/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement 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:6294:28: 6294 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:6354:72: 6354 | & -.5725014209747314D+00,.6074042001273483D+01/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement specfun.f:6356:72: 6356 | & .1727727502584457D+01,-.2438052969955606D+02/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement specfun.f:6358:72: 6358 | & .6765925884246826D+00,-.6883914268109947D+01/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement specfun.f:6360:72: 6360 | & -.1993531733751297D+01,.2724882731126854D+02/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement 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:6660:28: 6660 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 at (1) specfun.f:6723:72: 6723 | & -.3833534661393944D+13,.4854014686852901D+15/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement specfun.f:6729:72: 6729 | & .4218971570284096D+14,-.5827244631566907D+16/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement specfun.f:6735:72: 6735 | & .4011838599133198D+13,-.5060568503314727D+15/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement specfun.f:6741:72: 6741 | & -.4406481417852278D+14,.6065091351222699D+16/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement 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:7266:28: 7266 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 at (1) specfun.f:7298:29: 7298 | IMPLICIT COMPLEX *16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:7552:72: 7552 | & 1.796443723688307D-01,-1.39243221690590D+00/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement 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:7929:72: 7929 | & 5.5133589612202D02,3.0380905109224D03/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement specfun.f:7935:72: 7935 | & -6.0384407670507D02,-3.3022722944809D03/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement 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:7967:72: 7967 | & 2.3347645606175D05,1.2312234987631D07/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement 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:8000:28: 8000 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:8521:29: 8521 | IMPLICIT COMPLEX *16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:8768:28: 8768 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:9425:28: 9425 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:9967:28: 9967 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:10608:28: 10608 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:10770:28: 10770 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:10928:29: 10928 | IMPLICIT COMPLEX *16 (C,S,Z,D) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:11570:28: 11570 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:11773:28: 11773 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:11882:28: 11882 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 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:12091:28: 12091 | IMPLICIT COMPLEX*16 (C,Z) | 1 Warning: GNU Extension: Nonstandard type declaration COMPLEX*16 at (1) specfun.f:12133:72: 12133 | & 5.5133589612202D02,3.0380905109224D03/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement specfun.f:12139:72: 12139 | & -6.0384407670507D02,-3.3022722944809D03/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement 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:12171:72: 12171 | & 8.401390346421D08,7.2031420482627D10/ | 1 Warning: Obsolescent feature: DATA statement at (1) after the first executable statement 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) gcc -I"D:/RCompile/recent/R/include" -DNDEBUG -I"d:/rtools44/x86_64-w64-mingw32.static.posix/include" -pedantic -Wstrict-prototypes -O2 -Wall -mfpmath=sse -msse2 -mstackrealign -c weightedleastsquare.c -o weightedleastsquare.o gcc -shared -s -static-libgcc -o GeoModels.dll tmp.def 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 -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib/x64 -Ld:/rtools44/x86_64-w64-mingw32.static.posix/lib -lgfortran -lm -lquadmath -LD:/RCompile/recent/R/bin/x64 -lR D:\rtools44\x86_64-w64-mingw32.static.posix\bin/ld.exe: TB.o:TB.c:(.text+0x426): undefined reference to `dgemm_' D:\rtools44\x86_64-w64-mingw32.static.posix\bin/ld.exe: TB.o:TB.c:(.text+0x542): undefined reference to `dgemm_' D:\rtools44\x86_64-w64-mingw32.static.posix\bin/ld.exe: TB.o:TB.c:(.text+0x933): undefined reference to `zgemm_' D:\rtools44\x86_64-w64-mingw32.static.posix\bin/ld.exe: TB.o:TB.c:(.text+0xa47): undefined reference to `zgemm_' D:\rtools44\x86_64-w64-mingw32.static.posix\bin/ld.exe: TB.o:TB.c:(.text+0xcd4): undefined reference to `dnrm2_' D:\rtools44\x86_64-w64-mingw32.static.posix\bin/ld.exe: TB.o:TB.c:(.text+0x218c): undefined reference to `dsyev_' collect2.exe: error: ld returned 1 exit status no DLL was created ERROR: compilation failed for package 'GeoModels' * removing 'd:/RCompile/CRANincoming/R-devel/lib/GeoModels'