* installing *source* package ‘GeoModels’ ...
** package ‘GeoModels’ successfully unpacked and MD5 sums checked
** using staged installation
** libs
using C compiler: ‘clang version 18.1.4’
using Fortran compiler: ‘flang-new version 18.1.4’
make[1]: Entering directory '/data/gannet/ripley/R/packages/tests-clang/GeoModels/src'
/usr/local/clang18/bin/clang -I"/data/gannet/ripley/R/R-clang/include" -DNDEBUG   -isystem /usr/local/clang18/include -I/usr/local/clang/include    -fpic  -O3 -Wall -pedantic -Wp,-D_FORTIFY_SOURCE=3  -c CompositeLikelihood2.c -o CompositeLikelihood2.o
/usr/local/clang18/bin/clang -I"/data/gannet/ripley/R/R-clang/include" -DNDEBUG   -isystem /usr/local/clang18/include -I/usr/local/clang/include    -fpic  -O3 -Wall -pedantic -Wp,-D_FORTIFY_SOURCE=3  -c CompositeLikelihood2_ani.c -o CompositeLikelihood2_ani.o
/usr/local/clang18/bin/clang -I"/data/gannet/ripley/R/R-clang/include" -DNDEBUG   -isystem /usr/local/clang18/include -I/usr/local/clang/include    -fpic  -O3 -Wall -pedantic -Wp,-D_FORTIFY_SOURCE=3  -c CompositeLikelihoodCond2.c -o CompositeLikelihoodCond2.o
/usr/local/clang18/bin/clang -I"/data/gannet/ripley/R/R-clang/include" -DNDEBUG   -isystem /usr/local/clang18/include -I/usr/local/clang/include    -fpic  -O3 -Wall -pedantic -Wp,-D_FORTIFY_SOURCE=3  -c CompositeLikelihoodCond2_ani.c -o CompositeLikelihoodCond2_ani.o
/usr/local/clang18/bin/clang -I"/data/gannet/ripley/R/R-clang/include" -DNDEBUG   -isystem /usr/local/clang18/include -I/usr/local/clang/include    -fpic  -O3 -Wall -pedantic -Wp,-D_FORTIFY_SOURCE=3  -c CorrelationFunction.c -o CorrelationFunction.o
/usr/local/clang18/bin/clang -I"/data/gannet/ripley/R/R-clang/include" -DNDEBUG   -isystem /usr/local/clang18/include -I/usr/local/clang/include    -fpic  -O3 -Wall -pedantic -Wp,-D_FORTIFY_SOURCE=3  -c Distributions.c -o Distributions.o
/usr/local/clang18/bin/clang -I"/data/gannet/ripley/R/R-clang/include" -DNDEBUG   -isystem /usr/local/clang18/include -I/usr/local/clang/include    -fpic  -O3 -Wall -pedantic -Wp,-D_FORTIFY_SOURCE=3  -c GeoModels_init.c -o GeoModels_init.o
/usr/local/clang18/bin/clang -I"/data/gannet/ripley/R/R-clang/include" -DNDEBUG   -isystem /usr/local/clang18/include -I/usr/local/clang/include    -fpic  -O3 -Wall -pedantic -Wp,-D_FORTIFY_SOURCE=3  -c Godambe.c -o Godambe.o
Godambe.c:653:9: warning: variable 'nsub1' set but not used [-Wunused-but-set-variable]
  653 |     int nsub1 =0;
      |         ^
1 warning generated.
/usr/local/clang18/bin/clang -I"/data/gannet/ripley/R/R-clang/include" -DNDEBUG   -isystem /usr/local/clang18/include -I/usr/local/clang/include    -fpic  -O3 -Wall -pedantic -Wp,-D_FORTIFY_SOURCE=3  -c Gradient.c -o Gradient.o
/usr/local/clang18/bin/flang-new  -fpic  -O2 -pedantic  -c HypergeometricFunctions.f -o HypergeometricFunctions.o
./HypergeometricFunctions.f:43:7: portability: nonstandard usage: DOUBLE COMPLEX
        DOUBLE COMPLEX CHGF,A,B,Z,CONHYP 
        ^^^^^^^^^^^^^^
./HypergeometricFunctions.f:43:7: in the context: intrinsic type spec
        DOUBLE COMPLEX CHGF,A,B,Z,CONHYP 
        ^
./HypergeometricFunctions.f:43:7: in the context: declaration type spec
        DOUBLE COMPLEX CHGF,A,B,Z,CONHYP 
        ^
./HypergeometricFunctions.f:43:7: in the context: specification construct
        DOUBLE COMPLEX CHGF,A,B,Z,CONHYP 
        ^
./HypergeometricFunctions.f:43:7: in the context: declaration construct
        DOUBLE COMPLEX CHGF,A,B,Z,CONHYP 
        ^
./HypergeometricFunctions.f:42:7: in the context: specification part
        INTEGER LNCHF,I,BITS,IP 
        ^
./HypergeometricFunctions.f:40:7: in the context: FUNCTION subprogram
        FUNCTION CONHYP (A,B,Z,LNCHF,IP)                          
        ^
./HypergeometricFunctions.f:168:7: portability: nonstandard usage: DOUBLE COMPLEX
        DOUBLE COMPLEX A,B,Z,FINAL,CHGF                                       
        ^^^^^^^^^^^^^^
./HypergeometricFunctions.f:168:7: in the context: intrinsic type spec
        DOUBLE COMPLEX A,B,Z,FINAL,CHGF                                       
        ^
./HypergeometricFunctions.f:168:7: in the context: declaration type spec
        DOUBLE COMPLEX A,B,Z,FINAL,CHGF                                       
        ^
./HypergeometricFunctions.f:168:7: in the context: specification construct
        DOUBLE COMPLEX A,B,Z,FINAL,CHGF                                       
        ^
./HypergeometricFunctions.f:168:7: in the context: declaration construct
        DOUBLE COMPLEX A,B,Z,FINAL,CHGF                                       
        ^
./HypergeometricFunctions.f:166:7: in the context: specification part
        PARAMETER (LENGTH=777)
        ^
./HypergeometricFunctions.f:164:7: in the context: FUNCTION subprogram
        FUNCTION CHGF (A,B,Z,L,LNCHF)                          
        ^
./HypergeometricFunctions.f:675:7: portability: nonstandard usage: DOUBLE COMPLEX
        DOUBLE COMPLEX C                                                      
        ^^^^^^^^^^^^^^
./HypergeometricFunctions.f:675:7: in the context: intrinsic type spec
        DOUBLE COMPLEX C                                                      
        ^
./HypergeometricFunctions.f:675:7: in the context: declaration type spec
        DOUBLE COMPLEX C                                                      
        ^
./HypergeometricFunctions.f:675:7: in the context: specification construct
        DOUBLE COMPLEX C                                                      
        ^
./HypergeometricFunctions.f:675:7: in the context: declaration construct
        DOUBLE COMPLEX C                                                      
        ^
./HypergeometricFunctions.f:674:7: in the context: specification part
        INTEGER L,BIT,REXP,IR10,II10,LNCHF
        ^
./HypergeometricFunctions.f:672:7: in the context: SUBROUTINE subprogram
        SUBROUTINE ARYDIV(AR,AI,BR,BI,C,L,LNCHF,RMAX,BIT)                 
        ^
./HypergeometricFunctions.f:863:7: portability: nonstandard usage: DOUBLE COMPLEX
        DOUBLE COMPLEX CN                                                     
        ^^^^^^^^^^^^^^
./HypergeometricFunctions.f:863:7: in the context: intrinsic type spec
        DOUBLE COMPLEX CN                                                     
        ^
./HypergeometricFunctions.f:863:7: in the context: declaration type spec
        DOUBLE COMPLEX CN                                                     
        ^
./HypergeometricFunctions.f:863:7: in the context: specification construct
        DOUBLE COMPLEX CN                                                     
        ^
./HypergeometricFunctions.f:863:7: in the context: declaration construct
        DOUBLE COMPLEX CN                                                     
        ^
./HypergeometricFunctions.f:863:7: in the context: specification part
        DOUBLE COMPLEX CN                                                     
        ^
./HypergeometricFunctions.f:861:7: in the context: SUBROUTINE subprogram
        SUBROUTINE CONV12(CN,CAE)                                         
        ^
./HypergeometricFunctions.f:908:7: portability: nonstandard usage: DOUBLE COMPLEX
        DOUBLE COMPLEX CN                                                     
        ^^^^^^^^^^^^^^
./HypergeometricFunctions.f:908:7: in the context: intrinsic type spec
        DOUBLE COMPLEX CN                                                     
        ^
./HypergeometricFunctions.f:908:7: in the context: declaration type spec
        DOUBLE COMPLEX CN                                                     
        ^
./HypergeometricFunctions.f:908:7: in the context: specification construct
        DOUBLE COMPLEX CN                                                     
        ^
./HypergeometricFunctions.f:908:7: in the context: declaration construct
        DOUBLE COMPLEX CN                                                     
        ^
./HypergeometricFunctions.f:907:7: in the context: specification part
        DOUBLE PRECISION CAE                                             
        ^
./HypergeometricFunctions.f:905:7: in the context: SUBROUTINE subprogram
        SUBROUTINE CONV21(CAE,CN)                                         
        ^
./HypergeometricFunctions.f:987:7: portability: nonstandard usage: DOUBLE COMPLEX
        DOUBLE COMPLEX Z, A, B, CHF, CONHYP
        ^^^^^^^^^^^^^^
./HypergeometricFunctions.f:987:7: in the context: intrinsic type spec
        DOUBLE COMPLEX Z, A, B, CHF, CONHYP
        ^
./HypergeometricFunctions.f:987:7: in the context: declaration type spec
        DOUBLE COMPLEX Z, A, B, CHF, CONHYP
        ^
./HypergeometricFunctions.f:987:7: in the context: specification construct
        DOUBLE COMPLEX Z, A, B, CHF, CONHYP
        ^
./HypergeometricFunctions.f:987:7: in the context: declaration construct
        DOUBLE COMPLEX Z, A, B, CHF, CONHYP
        ^
./HypergeometricFunctions.f:985:7: in the context: specification part
        DOUBLE PRECISION ZRE(N), ZIM(N), CRE(N), CIM(N)
        ^
./HypergeometricFunctions.f:983:7: in the context: SUBROUTINE subprogram
        SUBROUTINE CHFM(ZRE,ZIM,ARE,AIM,BRE,BIM,CRE,CIM,N,LNCHF,IP)
        ^
/usr/local/clang18/bin/clang -I"/data/gannet/ripley/R/R-clang/include" -DNDEBUG   -isystem /usr/local/clang18/include -I/usr/local/clang/include    -fpic  -O3 -Wall -pedantic -Wp,-D_FORTIFY_SOURCE=3  -c TB.c -o TB.o
/usr/local/clang18/bin/clang -I"/data/gannet/ripley/R/R-clang/include" -DNDEBUG   -isystem /usr/local/clang18/include -I/usr/local/clang/include    -fpic  -O3 -Wall -pedantic -Wp,-D_FORTIFY_SOURCE=3  -c Utility.c -o Utility.o
/usr/local/clang18/bin/flang-new  -fpic  -O2 -pedantic  -c mvn.f -o mvn.o
/usr/local/clang18/bin/flang-new  -fpic  -O2 -pedantic  -c mvndstpack.f -o mvndstpack.o
/usr/local/clang18/bin/flang-new  -fpic  -O2 -pedantic  -c specfun.f -o specfun.o
./specfun.f:49:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:49:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:49:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:49:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:48:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A-B,D-H,O-Y)
        ^
./specfun.f:48:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A-B,D-H,O-Y)
        ^
./specfun.f:37:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CPDSA(N,Z,CDN)
        ^
./specfun.f:102:26: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX *16 (C,S,Z,D)
                           ^^^
./specfun.f:102:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX *16 (C,S,Z,D)
                   ^
./specfun.f:102:18: in the context: declaration type spec
          IMPLICIT COMPLEX *16 (C,S,Z,D)
                   ^
./specfun.f:102:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX *16 (C,S,Z,D)
          ^
./specfun.f:101:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (E,P,W)
        ^
./specfun.f:101:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (E,P,W)
        ^
./specfun.f:92:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CFS(Z,ZF,ZD)
        ^
./specfun.f:281:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:281:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:281:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:281:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:280:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (D,X,Y)
        ^
./specfun.f:280:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (D,X,Y)
        ^
./specfun.f:264:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CLPMN(MM,M,N,X,Y,NTYPE,CPM,CPD)
        ^
./specfun.f:1245:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:1245:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:1245:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:1245:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:1244:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:1244:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:1228:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CJYLV(V,Z,CBJV,CDJV,CBYV,CDYV)
        ^
./specfun.f:1582:17: portability: nonstandard usage: TYPE*KIND syntax
          COMPLEX *16 Z,CER,CDER
                  ^^^
./specfun.f:1582:9: in the context: intrinsic type spec
          COMPLEX *16 Z,CER,CDER
          ^
./specfun.f:1582:9: in the context: declaration type spec
          COMPLEX *16 Z,CER,CDER
          ^
./specfun.f:1582:7: in the context: specification construct
          COMPLEX *16 Z,CER,CDER
        ^
./specfun.f:1582:7: in the context: declaration construct
          COMPLEX *16 Z,CER,CDER
        ^
./specfun.f:1581:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A-H,O-Z)
        ^
./specfun.f:1571:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CERF(Z,CER,CDER)
        ^
./specfun.f:2007:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:2007:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:2007:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:2007:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:2006:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (X,Y)
        ^
./specfun.f:2006:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (X,Y)
        ^
./specfun.f:1993:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CLQN(N,X,Y,CQN,CQD)
        ^
./specfun.f:2216:26: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX *16 (C,Z)
                           ^^^
./specfun.f:2216:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX *16 (C,Z)
                   ^
./specfun.f:2216:18: in the context: declaration type spec
          IMPLICIT COMPLEX *16 (C,Z)
                   ^
./specfun.f:2216:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX *16 (C,Z)
          ^
./specfun.f:2216:7: in the context: implicit part
          IMPLICIT COMPLEX *16 (C,Z)
        ^
./specfun.f:2216:7: in the context: specification part
          IMPLICIT COMPLEX *16 (C,Z)
        ^
./specfun.f:2207:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CERROR(Z,CER)
        ^
./specfun.f:3099:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:3099:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:3099:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:3099:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:3098:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A,B,G,O-Y)
        ^
./specfun.f:3098:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A,B,G,O-Y)
        ^
./specfun.f:3079:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CJYVA(V,Z,VM,CBJ,CDJ,CBY,CDY)
        ^
./specfun.f:3363:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:3363:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:3363:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:3363:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:3362:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A,B,G,O-Y)
        ^
./specfun.f:3362:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A,B,G,O-Y)
        ^
./specfun.f:3343:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CJYVB(V,Z,VM,CBJ,CDJ,CBY,CDY)
        ^
./specfun.f:4567:26: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX *16 (C,Z)
                           ^^^
./specfun.f:4567:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX *16 (C,Z)
                   ^
./specfun.f:4567:18: in the context: declaration type spec
          IMPLICIT COMPLEX *16 (C,Z)
                   ^
./specfun.f:4567:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX *16 (C,Z)
          ^
./specfun.f:4566:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (E,P,W)
        ^
./specfun.f:4566:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (E,P,W)
        ^
./specfun.f:4556:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CERZO(NT,ZO)
        ^
./specfun.f:4963:26: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX *16 (C,Z)
                           ^^^
./specfun.f:4963:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX *16 (C,Z)
                   ^
./specfun.f:4963:18: in the context: declaration type spec
          IMPLICIT COMPLEX *16 (C,Z)
                   ^
./specfun.f:4963:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX *16 (C,Z)
          ^
./specfun.f:4962:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (X,Y)
        ^
./specfun.f:4962:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (X,Y)
        ^
./specfun.f:4949:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CLPN(N,X,Y,CPN,CPD)
        ^
./specfun.f:5123:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:5123:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:5123:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:5123:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:5122:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:5122:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:5106:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CIKLV(V,Z,CBIV,CDIV,CBKV,CDKV)
        ^
./specfun.f:5317:9: portability: nonstandard usage: DOUBLE COMPLEX
          DOUBLE COMPLEX Z, CEI
          ^^^^^^^^^^^^^^
./specfun.f:5317:9: in the context: intrinsic type spec
          DOUBLE COMPLEX Z, CEI
          ^
./specfun.f:5317:9: in the context: declaration type spec
          DOUBLE COMPLEX Z, CEI
          ^
./specfun.f:5317:7: in the context: specification construct
          DOUBLE COMPLEX Z, CEI
        ^
./specfun.f:5317:7: in the context: declaration construct
          DOUBLE COMPLEX Z, CEI
        ^
./specfun.f:5316:7: in the context: specification part
          IMPLICIT NONE
        ^
./specfun.f:5308:7: in the context: SUBROUTINE subprogram
          SUBROUTINE EIXZ(Z,CEI)
        ^
./specfun.f:5384:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C)
                          ^^^
./specfun.f:5384:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C)
                   ^
./specfun.f:5384:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C)
                   ^
./specfun.f:5384:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C)
          ^
./specfun.f:5383:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A-B,D-H,O-Z)
        ^
./specfun.f:5383:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A-B,D-H,O-Z)
        ^
./specfun.f:5371:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CHGM(A,B,X,HG)
        ^
./specfun.f:5683:26: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX *16 (C,Z)
                           ^^^
./specfun.f:5683:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX *16 (C,Z)
                   ^
./specfun.f:5683:18: in the context: declaration type spec
          IMPLICIT COMPLEX *16 (C,Z)
                   ^
./specfun.f:5683:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX *16 (C,Z)
          ^
./specfun.f:5682:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:5682:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:5669:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CCHG(A,B,Z,CHG)
        ^
./specfun.f:5811:26: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX *16 (Z)
                           ^^^
./specfun.f:5811:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX *16 (Z)
                   ^
./specfun.f:5811:18: in the context: declaration type spec
          IMPLICIT COMPLEX *16 (Z)
                   ^
./specfun.f:5811:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX *16 (Z)
          ^
./specfun.f:5810:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A-H,O-Y)
        ^
./specfun.f:5810:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A-H,O-Y)
        ^
./specfun.f:5794:7: in the context: SUBROUTINE subprogram
          SUBROUTINE HYGFZ(A,B,C,Z,ZHF,ISFER)
        ^
./specfun.f:6294:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:6294:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:6294:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:6294:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:6293:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:6293:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:6276:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CJYNB(N,Z,NM,CBJ,CDJ,CBY,CDY)
        ^
./specfun.f:6660:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:6660:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:6660:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:6660:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:6659:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A,B,E,P,R,W)
        ^
./specfun.f:6659:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A,B,E,P,R,W)
        ^
./specfun.f:6645:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CY01(KF,Z,ZF,ZD)
        ^
./specfun.f:7266:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:7266:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:7266:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:7266:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:7265:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A-B,D-H,O-Y)
        ^
./specfun.f:7265:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A-B,D-H,O-Y)
        ^
./specfun.f:7255:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CPDLA(N,Z,CDN)
        ^
./specfun.f:7298:26: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX *16 (C,Z)
                           ^^^
./specfun.f:7298:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX *16 (C,Z)
                   ^
./specfun.f:7298:18: in the context: declaration type spec
          IMPLICIT COMPLEX *16 (C,Z)
                   ^
./specfun.f:7298:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX *16 (C,Z)
          ^
./specfun.f:7297:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (E,P,W)
        ^
./specfun.f:7297:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (E,P,W)
        ^
./specfun.f:7283:7: in the context: SUBROUTINE subprogram
          SUBROUTINE FCSZO(KF,NT,ZO)
        ^
./specfun.f:8000:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:8000:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:8000:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:8000:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:7999:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A-B,D-H,O-Y)
        ^
./specfun.f:7999:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A-B,D-H,O-Y)
        ^
./specfun.f:7985:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CPBDN(N,Z,CPB,CPD)
        ^
./specfun.f:8521:26: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX *16 (C,Z)
                           ^^^
./specfun.f:8521:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX *16 (C,Z)
                   ^
./specfun.f:8521:18: in the context: declaration type spec
          IMPLICIT COMPLEX *16 (C,Z)
                   ^
./specfun.f:8521:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX *16 (C,Z)
          ^
./specfun.f:8520:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:8520:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:8505:7: in the context: SUBROUTINE subprogram
          SUBROUTINE PBWA(A,X,W1F,W1D,W2F,W2D)
        ^
./specfun.f:8768:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:8768:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:8768:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:8768:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:8768:7: in the context: implicit part
          IMPLICIT COMPLEX*16 (C,Z)
        ^
./specfun.f:8768:7: in the context: specification part
          IMPLICIT COMPLEX*16 (C,Z)
        ^
./specfun.f:8760:7: in the context: SUBROUTINE subprogram
          SUBROUTINE E1Z(Z,CE1)
        ^
./specfun.f:9425:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:9425:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:9425:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:9425:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:9424:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (H,O-Y)
        ^
./specfun.f:9424:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (H,O-Y)
        ^
./specfun.f:9404:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CYZO(NT,KF,KC,ZO,ZV)
        ^
./specfun.f:9967:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:9967:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:9967:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:9967:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:9966:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:9966:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:9948:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CH12N(N,Z,NM,CHF1,CHD1,CHF2,CHD2)
        ^
./specfun.f:10608:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:10608:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:10608:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:10608:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:10607:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A,D-H,O-Y)
        ^
./specfun.f:10607:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A,D-H,O-Y)
        ^
./specfun.f:10587:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CIKVB(V,Z,VM,CBI,CDI,CBK,CDK)
        ^
./specfun.f:10770:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:10770:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:10770:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:10770:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:10769:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A,G,P,R,V,W)
        ^
./specfun.f:10769:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A,G,P,R,V,W)
        ^
./specfun.f:10749:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CIKVA(V,Z,VM,CBI,CDI,CBK,CDK)
        ^
./specfun.f:10928:26: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX *16 (C,S,Z,D)
                           ^^^
./specfun.f:10928:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX *16 (C,S,Z,D)
                   ^
./specfun.f:10928:18: in the context: declaration type spec
          IMPLICIT COMPLEX *16 (C,S,Z,D)
                   ^
./specfun.f:10928:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX *16 (C,S,Z,D)
          ^
./specfun.f:10927:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (E,P,W)
        ^
./specfun.f:10927:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (E,P,W)
        ^
./specfun.f:10918:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CFC(Z,ZF,ZD)
        ^
./specfun.f:11570:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:11570:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:11570:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:11570:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:11569:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (X,Y)
        ^
./specfun.f:11569:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (X,Y)
        ^
./specfun.f:11554:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CLQMN(MM,M,N,X,Y,CQM,CQD)
        ^
./specfun.f:11773:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:11773:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:11773:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:11773:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:11772:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:11772:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:11755:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CIKNB(N,Z,NM,CBI,CDI,CBK,CDK)
        ^
./specfun.f:11882:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:11882:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:11882:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:11882:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:11881:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A,B,P,W,X,Y)
        ^
./specfun.f:11881:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A,B,P,W,X,Y)
        ^
./specfun.f:11863:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CIKNA(N,Z,NM,CBI,CDI,CBK,CDK)
        ^
./specfun.f:12091:25: portability: nonstandard usage: TYPE*KIND syntax
          IMPLICIT COMPLEX*16 (C,Z)
                          ^^^
./specfun.f:12091:18: in the context: intrinsic type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:12091:18: in the context: declaration type spec
          IMPLICIT COMPLEX*16 (C,Z)
                   ^
./specfun.f:12091:9: in the context: IMPLICIT statement
          IMPLICIT COMPLEX*16 (C,Z)
          ^
./specfun.f:12090:7: in the context: implicit part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:12090:7: in the context: specification part
          IMPLICIT DOUBLE PRECISION (A,B,D-H,O-Y)
        ^
./specfun.f:12073:7: in the context: SUBROUTINE subprogram
          SUBROUTINE CIK01(Z,CBI0,CDI0,CBI1,CDI1,CBK0,CDK0,CBK1,CDK1)
        ^
/usr/local/clang18/bin/clang -I"/data/gannet/ripley/R/R-clang/include" -DNDEBUG   -isystem /usr/local/clang18/include -I/usr/local/clang/include    -fpic  -O3 -Wall -pedantic -Wp,-D_FORTIFY_SOURCE=3  -c weightedleastsquare.c -o weightedleastsquare.o
./specfun.f:115:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=1,80
             ^^^^^^^^^^^^
./specfun.f:121:1: DO loop currently ends at statement:
  10            WB0=WB
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:127:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=M,0,-1
             ^^^^^^^^^^^^^^
./specfun.f:131:1: DO loop currently ends at statement:
  15            CF0=CF
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:155:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=1,20
             ^^^^^^^^^^^^
./specfun.f:157:1: DO loop currently ends at statement:
  20            CF=CF+CR
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:160:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=1,12
             ^^^^^^^^^^^^
./specfun.f:162:1: DO loop currently ends at statement:
  25            CG=CG+CR
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:225:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=KM,0,-1
             ^^^^^^^^^^^^^^^
./specfun.f:229:1: DO loop currently ends at statement:
  25            QF1=QF0
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:230:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=0,N
             ^^^^^^^^^^^
./specfun.f:231:1: DO loop currently ends at statement:
  30            QM(0,K)=Q0*QM(0,K)/QF0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:234:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=KM,0,-1
             ^^^^^^^^^^^^^^^
./specfun.f:238:1: DO loop currently ends at statement:
  35            QF1=QF0
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:240:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 40 K=0,N
             ^^^^^^^^^^^
./specfun.f:241:1: DO loop currently ends at statement:
  40            QM(1,K)=Q10*QM(1,K)/QF0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:253:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 50 J=1,N
          ^^^^^^^^^^^
./specfun.f:254:1: DO loop currently ends at statement:
  50         QD(0,J)=LS*J*(QM(0,J-1)-X*QM(0,J))/XS
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:284:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 I=0,N
          ^^^^^^^^^^^
./specfun.f:287:1: DO loop currently ends at statement:
  10         CPD(J,I)=(0.0D0,0.0D0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:285:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 J=0,M
          ^^^^^^^^^^^
./specfun.f:287:1: DO loop currently ends at statement:
  10         CPD(J,I)=(0.0D0,0.0D0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:291:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 I=1,N
             ^^^^^^^^^^^
./specfun.f:293:1: DO loop currently ends at statement:
  15            CPD(0,I)=0.5D0*I*(I+1)*X**(I+1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:318:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 25 I=1,M
          ^^^^^^^^^^^
./specfun.f:320:1: DO loop currently ends at statement:
  25         CPM(I,I)=(2.0D0*I-1.0D0)*ZQ*CPM(I-1,I-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:321:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 30 I=0,MIN(M,N-1)
          ^^^^^^^^^^^^^^^^^^^^
./specfun.f:323:1: DO loop currently ends at statement:
  30         CPM(I,I+1)=(2.0D0*I+1.0D0)*Z*CPM(I,I)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:331:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 40 J=1,N
          ^^^^^^^^^^^
./specfun.f:333:1: DO loop currently ends at statement:
  40         CPD(0,J)=LS*J*(Z*CPM(0,J)-CPM(0,J-1))/ZS
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:441:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 45 I=1,NM
          ^^^^^^^^^^^^
./specfun.f:506:1: DO loop currently ends at statement:
  45         L0=L2
  ^^^^^^^^^^^^^^^^
./specfun.f:484:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 40 K=1,L
                ^^^^^^^^^^^
./specfun.f:488:1: DO loop currently ends at statement:
  40               P(K)=P1(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:529:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 J=2,N2
          ^^^^^^^^^^^^
./specfun.f:530:1: DO loop currently ends at statement:
  10         U(J)=C*C
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:531:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 15 J=1,N2
          ^^^^^^^^^^^^
./specfun.f:532:1: DO loop currently ends at statement:
  15         V(J)=(2.0*J-1.0-IP)*(2.0*(J-M)-IP)+M*(M-1.0)-CV
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:533:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 20 J=1,NM-1
          ^^^^^^^^^^^^^^
./specfun.f:534:1: DO loop currently ends at statement:
  20         W(J)=(2.0*J-IP)*(2.0*J+1.0-IP)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:543:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 25 J=1,K
                   ^^^^^^^^^^^
./specfun.f:544:1: DO loop currently ends at statement:
  25                  R1=R1*(I+M-J)/J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:559:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 45 J=1,K
                   ^^^^^^^^^^^
./specfun.f:560:1: DO loop currently ends at statement:
  45                  R1=R1*(I+M-J)/J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:571:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 65 K=2,N2
          ^^^^^^^^^^^^
./specfun.f:574:1: DO loop currently ends at statement:
  65         BK(K)=(BK(K)-BK(K-1)*U(K))/T
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:575:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 70 K=N2-1,1,-1
          ^^^^^^^^^^^^^^^^^
./specfun.f:576:1: DO loop currently ends at statement:
  70         BK(K)=BK(K)-W(K)*BK(K+1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:617:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:621:1: DO loop currently ends at statement:
  10        SW=SU0
  ^^^^^^^^^^^^^^^^
./specfun.f:623:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 20 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:627:1: DO loop currently ends at statement:
  20        SW=SD0
  ^^^^^^^^^^^^^^^^
./specfun.f:630:11: portability: A DO loop should terminate with an END DO or CONTINUE
            DO 30 K=1,M
            ^^^^^^^^^^^
./specfun.f:634:1: DO loop currently ends at statement:
  30           SD1=SD1+DN(K)*QD(J)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:641:14: portability: A DO loop should terminate with an END DO or CONTINUE
               DO 35 J1=1,J
               ^^^^^^^^^^^^
./specfun.f:642:1: DO loop currently ends at statement:
  35              R1=(M+J1)*R1
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:644:14: portability: A DO loop should terminate with an END DO or CONTINUE
               DO 40 J2=1,M-J-2
               ^^^^^^^^^^^^^^^^
./specfun.f:645:1: DO loop currently ends at statement:
  40              R2=J2*R2
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:648:14: portability: A DO loop should terminate with an END DO or CONTINUE
               DO 45 L1=1,J
               ^^^^^^^^^^^^
./specfun.f:650:1: DO loop currently ends at statement:
  45              SF=SF+R3
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:659:14: portability: A DO loop should terminate with an END DO or CONTINUE
               DO 50 L1=1,J-1
               ^^^^^^^^^^^^^^
./specfun.f:662:1: DO loop currently ends at statement:
  50              SD=SD+R4
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:669:11: portability: A DO loop should terminate with an END DO or CONTINUE
            DO 60 K=KI,NM3
            ^^^^^^^^^^^^^^
./specfun.f:673:1: DO loop currently ends at statement:
  60           SW=SU2
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:675:11: portability: A DO loop should terminate with an END DO or CONTINUE
            DO 70 K=KI,NM3
            ^^^^^^^^^^^^^^
./specfun.f:679:1: DO loop currently ends at statement:
  70           SW=SD2
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:706:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 20 M=4,N,2
          ^^^^^^^^^^^^^
./specfun.f:714:1: DO loop currently ends at statement:
  20         BN(M)=R1*R2
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:732:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 30 M=2,N
          ^^^^^^^^^^^
./specfun.f:739:1: DO loop currently ends at statement:
  30         BN(M)=S
  ^^^^^^^^^^^^^^^^^^
./specfun.f:734:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=2,M-1
             ^^^^^^^^^^^^^
./specfun.f:738:1: DO loop currently ends at statement:
  20            S=S-R*BN(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:736:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 10 J=2,K
                ^^^^^^^^^^^
./specfun.f:737:1: DO loop currently ends at statement:
  10               R=R*(J+M-K)/J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:740:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 40 M=3,N,2
          ^^^^^^^^^^^^^
./specfun.f:741:1: DO loop currently ends at statement:
  40         BN(M)=0.0D0
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:760:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 20 I=1,M
          ^^^^^^^^^^^
./specfun.f:767:1: DO loop currently ends at statement:
  20      AP(I+1)=-R*S
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:762:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 L=1,I
             ^^^^^^^^^^^
./specfun.f:766:1: DO loop currently ends at statement:
  15            S=S+SK*AP(I-L+1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:764:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 10 K=0,L
                ^^^^^^^^^^^
./specfun.f:765:1: DO loop currently ends at statement:
  10               SK=SK+CK(K+1)*CK(L-K+1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:769:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 30 L=1,M
          ^^^^^^^^^^^
./specfun.f:773:1: DO loop currently ends at statement:
  30         QS0=QS0+AP(M-L+1)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:771:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=1,L
             ^^^^^^^^^^^
./specfun.f:772:1: DO loop currently ends at statement:
  25            R=R*(2.0D0*K+IP)*(2.0D0*K-1.0D0+IP)/(2.0D0*K)**2
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1022:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 IT=1,20
          ^^^^^^^^^^^^^
./specfun.f:1029:2: DO loop currently ends at statement:
   10        F1=F
   ^^^^^^^^^^^^^^
./specfun.f:1061:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 IT=1,20
          ^^^^^^^^^^^^^
./specfun.f:1068:1: DO loop currently ends at statement:
  10         F1=F
  ^^^^^^^^^^^^^^^
./specfun.f:1212:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 55 K=1,10
             ^^^^^^^^^^^^
./specfun.f:1214:1: DO loop currently ends at statement:
  55            G0=G0+R0
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1217:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 60 K=1,10
             ^^^^^^^^^^^^
./specfun.f:1219:1: DO loop currently ends at statement:
  60            G1=G1+R1
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1257:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=1,KM
             ^^^^^^^^^^^^
./specfun.f:1263:1: DO loop currently ends at statement:
  15            CF(K)=CF(K)*CT**K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1261:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 10 I=LF-1,L0,-1
                ^^^^^^^^^^^^^^^^^^
./specfun.f:1262:1: DO loop currently ends at statement:
  10               CF(K)=CF(K)*CT2+A(I)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1266:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=1,KM
             ^^^^^^^^^^^^
./specfun.f:1267:1: DO loop currently ends at statement:
  20            CSJ=CSJ+CF(K)*VR**K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1271:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=1,KM
             ^^^^^^^^^^^^
./specfun.f:1272:1: DO loop currently ends at statement:
  25            CSY=CSY+(-1)**K*CF(K)*VR**K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1309:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 J=1,2*M+IP
          ^^^^^^^^^^^^^^^^
./specfun.f:1310:1: DO loop currently ends at statement:
  10         R0=R0*J
  ^^^^^^^^^^^^^^^^^^
./specfun.f:1314:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 15 K=2,NM
          ^^^^^^^^^^^^
./specfun.f:1318:1: DO loop currently ends at statement:
  15         SW=SUC
  ^^^^^^^^^^^^^^^^^
./specfun.f:1323:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 50 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:1336:1: DO loop currently ends at statement:
  50         SW=R2F
  ^^^^^^^^^^^^^^^^^
./specfun.f:1346:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 60 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:1359:1: DO loop currently ends at statement:
  60         SW=SUD
  ^^^^^^^^^^^^^^^^^
./specfun.f:1388:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=1 ,N-1
             ^^^^^^^^^^^^^^
./specfun.f:1389:1: DO loop currently ends at statement:
  10            S=S+1.0D0/K
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1393:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=1,N
             ^^^^^^^^^^^
./specfun.f:1394:1: DO loop currently ends at statement:
  20            S=S+1.0/(2.0D0*K-1.0D0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1399:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 30 K=0,N-1
                ^^^^^^^^^^^^^
./specfun.f:1400:1: DO loop currently ends at statement:
  30               S=S+1.0D0/(XA+K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1525:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 K=0,N
          ^^^^^^^^^^^
./specfun.f:1527:1: DO loop currently ends at statement:
  10         PD(K)=0.0D0
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1549:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 20 K=1,M
          ^^^^^^^^^^^
./specfun.f:1551:1: DO loop currently ends at statement:
  20         PM0=PMK
  ^^^^^^^^^^^^^^^^^^
./specfun.f:1555:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 25 K=M+2,N
          ^^^^^^^^^^^^^
./specfun.f:1559:1: DO loop currently ends at statement:
  25         PM1=PM2
  ^^^^^^^^^^^^^^^^^^
./specfun.f:1561:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 30 K=1,N
          ^^^^^^^^^^^
./specfun.f:1562:1: DO loop currently ends at statement:
  30         PD(K)=(K*X*PM(K)-(K+M)*PM(K-1))/(X*X-1.0D0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1563:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 35 K=1,N
          ^^^^^^^^^^^
./specfun.f:1565:1: DO loop currently ends at statement:
  35         PD(K)=(-1)**M*PD(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1592:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=1,100
             ^^^^^^^^^^^^^
./specfun.f:1596:1: DO loop currently ends at statement:
  10            W=ER
  ^^^^^^^^^^^^^^^^^^
./specfun.f:1602:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=1,12
             ^^^^^^^^^^^^
./specfun.f:1604:1: DO loop currently ends at statement:
  20            ER=ER+R
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1618:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 N=1,100
             ^^^^^^^^^^^^^
./specfun.f:1622:1: DO loop currently ends at statement:
  25            W1=ER2
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:1627:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 N=1,100
             ^^^^^^^^^^^^^
./specfun.f:1631:1: DO loop currently ends at statement:
  35            W2=EI2
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:1747:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 20 K=24,1,-1
          ^^^^^^^^^^^^^^^
./specfun.f:1748:1: DO loop currently ends at statement:
  20         GR=GR*X+G(K)
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1804:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 30 M=1,N/2
          ^^^^^^^^^^^^^
./specfun.f:1811:1: DO loop currently ends at statement:
  30         EN(2*M)=-S
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1806:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=1,M-1
             ^^^^^^^^^^^^^
./specfun.f:1810:1: DO loop currently ends at statement:
  20            S=S+R*EN(2*K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1808:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 10 J=1,2*K
                ^^^^^^^^^^^^^
./specfun.f:1809:1: DO loop currently ends at statement:
  10               R=R*(2.0D0*M-2.0D0*K+J)/J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1839:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 IT=1,100
          ^^^^^^^^^^^^^^
./specfun.f:1847:1: DO loop currently ends at statement:
  10         F1=F
  ^^^^^^^^^^^^^^^
./specfun.f:1895:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=M,1,-1
             ^^^^^^^^^^^^^^
./specfun.f:1899:1: DO loop currently ends at statement:
  25            XA0=XA
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:1901:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=3,M,2
             ^^^^^^^^^^^^^
./specfun.f:1902:1: DO loop currently ends at statement:
  30            XS=XS+2.0D0*BJ(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1904:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=2,M
             ^^^^^^^^^^^
./specfun.f:1905:1: DO loop currently ends at statement:
  35            BJ(K)=BJ(K)/XS
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1908:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 40 K=2,M
             ^^^^^^^^^^^
./specfun.f:1910:1: DO loop currently ends at statement:
  40            XG1=XG1+BJ(K)*XR
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1913:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 45 K=2,M
             ^^^^^^^^^^^
./specfun.f:1915:1: DO loop currently ends at statement:
  45            XG2=XG2+BJ(K)*XR
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1923:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 50 K=1,9
             ^^^^^^^^^^^
./specfun.f:1925:1: DO loop currently ends at statement:
  50            XF=XF+XR
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1928:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 55 K=1,8
             ^^^^^^^^^^^
./specfun.f:1930:1: DO loop currently ends at statement:
  55            XG=XG+XR
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:1975:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=1,10
             ^^^^^^^^^^^^
./specfun.f:1980:1: DO loop currently ends at statement:
  30            A1=AF
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:1983:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=1,11
             ^^^^^^^^^^^^
./specfun.f:1985:1: DO loop currently ends at statement:
  35            TI=TI+A(K)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2012:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:2014:1: DO loop currently ends at statement:
  10            CQD(K)=(1.0D+300,0.0D0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2026:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=2,N
             ^^^^^^^^^^^
./specfun.f:2030:1: DO loop currently ends at statement:
  15            CQF1=CQF2
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2039:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=KM,0,-1
             ^^^^^^^^^^^^^^^
./specfun.f:2043:1: DO loop currently ends at statement:
  20            CQF1=CQF0
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2044:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=0,N
             ^^^^^^^^^^^
./specfun.f:2045:1: DO loop currently ends at statement:
  25            CQN(K)=CQN(K)*CQ0/CQF0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2048:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 30 K=1,N
          ^^^^^^^^^^^
./specfun.f:2049:1: DO loop currently ends at statement:
  30         CQD(K)=(K*Z*CQN(K)-K*CQN(K-1))/(Z*Z-1.0D0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2195:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=1,12
             ^^^^^^^^^^^^
./specfun.f:2197:1: DO loop currently ends at statement:
  20            ER=ER+R
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2287:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 20 M=4,N,2
          ^^^^^^^^^^^^^
./specfun.f:2297:1: DO loop currently ends at statement:
  20         EN(M)=R1*R2
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2333:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 10 IC=1,ICM
                ^^^^^^^^^^^^^^
./specfun.f:2334:1: DO loop currently ends at statement:
  10               CV(IC)=4.0D0*(IC-1.0D0)**2
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2336:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 15 IC=1,ICM
                ^^^^^^^^^^^^^^
./specfun.f:2337:1: DO loop currently ends at statement:
  15               CV(IC)=(2.0D0*IC-1.0D0)**2
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2339:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 20 IC=1,ICM
                ^^^^^^^^^^^^^^
./specfun.f:2340:1: DO loop currently ends at statement:
  20               CV(IC)=4.0D0*IC*IC
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2348:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 25 I=2,NM
                ^^^^^^^^^^^^
./specfun.f:2351:1: DO loop currently ends at statement:
  25               F(I)=Q*Q
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2356:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 30 I=2,NM
                ^^^^^^^^^^^^
./specfun.f:2359:1: DO loop currently ends at statement:
  30               F(I)=Q*Q
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2362:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 35 I=2,NM
                ^^^^^^^^^^^^
./specfun.f:2365:1: DO loop currently ends at statement:
  35               F(I)=Q*Q
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2377:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 45 I=1,ICM
             ^^^^^^^^^^^^^
./specfun.f:2379:1: DO loop currently ends at statement:
  45            H(I)=XB
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2488:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:2490:1: DO loop currently ends at statement:
  10            QD(K)=1.0D+300
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2501:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=2,N
             ^^^^^^^^^^^
./specfun.f:2506:1: DO loop currently ends at statement:
  15            Q1=QF
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:2532:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 40 K=N,2,-1
             ^^^^^^^^^^^^^^
./specfun.f:2536:1: DO loop currently ends at statement:
  40            QF1=QF0
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2538:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 45 K=1,N
             ^^^^^^^^^^^
./specfun.f:2539:1: DO loop currently ends at statement:
  45            QD(K)=K*(QN(K-1)-X*QN(K))/(1.0D0-X*X)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2562:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 K=0,KM-1
          ^^^^^^^^^^^^^^
./specfun.f:2570:1: DO loop currently ends at statement:
  10         G0=G
  ^^^^^^^^^^^^^^^
./specfun.f:2619:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=1,8
             ^^^^^^^^^^^
./specfun.f:2621:1: DO loop currently ends at statement:
  20            TTI=TTI+C(K)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2642:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=1,8
             ^^^^^^^^^^^
./specfun.f:2644:1: DO loop currently ends at statement:
  35            TTK=TTK+C(K)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2679:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=0,N
             ^^^^^^^^^^^
./specfun.f:2696:1: DO loop currently ends at statement:
  25            DL(K)=-0.5D0*X/(VK+1.0D0)*UK
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2708:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:2711:1: DO loop currently ends at statement:
  30            PX=PX+RP
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2714:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:2717:1: DO loop currently ends at statement:
  35            QX=QX+RQ
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2748:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 45 K=2,N
             ^^^^^^^^^^^
./specfun.f:2752:1: DO loop currently ends at statement:
  45            VL(K)=F
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2763:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 50 K=M,0,-1
             ^^^^^^^^^^^^^^
./specfun.f:2767:1: DO loop currently ends at statement:
  50            F1=F
  ^^^^^^^^^^^^^^^^^^
./specfun.f:2774:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 55 K=0,N
             ^^^^^^^^^^^
./specfun.f:2775:1: DO loop currently ends at statement:
  55            VL(K)=CS*VL(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2778:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 65 J=1,N
          ^^^^^^^^^^^
./specfun.f:2782:1: DO loop currently ends at statement:
  65         R0=2.0D0*(J+V0+1)/X*R0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2914:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 I=1,NN+3
          ^^^^^^^^^^^^^^
./specfun.f:2923:1: DO loop currently ends at statement:
  10         W(I)=(K+1.0D0)*(K+2.0D0)*CS/((GK2+2.0D0)*GK3)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2924:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 20 K=1,M
          ^^^^^^^^^^^
./specfun.f:2928:1: DO loop currently ends at statement:
  20         RK(K)=-U(K)/T
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2926:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 L=0,M-K-1
             ^^^^^^^^^^^^^^^
./specfun.f:2927:1: DO loop currently ends at statement:
  15            T=V(M-L)-W(M-L+1)*U(M-L)/T
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2930:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 25 K=1,M
          ^^^^^^^^^^^
./specfun.f:2932:1: DO loop currently ends at statement:
  25         DN(K)=DF(1)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2942:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 35 K=M+2,NN
          ^^^^^^^^^^^^^^
./specfun.f:2943:1: DO loop currently ends at statement:
  35         DN(K)=RK(K)*DN(K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2945:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 40 J=1,(N+M+IP)/2
          ^^^^^^^^^^^^^^^^^^^^
./specfun.f:2946:1: DO loop currently ends at statement:
  40         R1=R1*(J+0.5D0*(N+M+IP))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2949:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 45 J=1,2*M+IP
          ^^^^^^^^^^^^^^^^
./specfun.f:2950:1: DO loop currently ends at statement:
  45         R=R*J
  ^^^^^^^^^^^^^^^^
./specfun.f:2953:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 50 K=2,NM
          ^^^^^^^^^^^^
./specfun.f:2957:1: DO loop currently ends at statement:
  50         SW=SU0
  ^^^^^^^^^^^^^^^^^
./specfun.f:2960:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 60 J=1,M
          ^^^^^^^^^^^
./specfun.f:2961:1: DO loop currently ends at statement:
  60         R2=2.0D0*C*R2*J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2963:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 65 J=1,(N-M-IP)/2
          ^^^^^^^^^^^^^^^^^^^^
./specfun.f:2964:1: DO loop currently ends at statement:
  65         R3=R3*J
  ^^^^^^^^^^^^^^^^^^
./specfun.f:2969:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 75 J=1,(N-M-IP)/2
          ^^^^^^^^^^^^^^^^^^^^
./specfun.f:2970:1: DO loop currently ends at statement:
  75         R4=4.0D0*R4*J
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:2972:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 80 J=1,M
          ^^^^^^^^^^^
./specfun.f:2973:1: DO loop currently ends at statement:
  80         R5=R5*(J+M)/C
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3009:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 I=1,NR-1
             ^^^^^^^^^^^^^^
./specfun.f:3010:1: DO loop currently ends at statement:
  15            P=P*(Z-X(I))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3013:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=2,N
             ^^^^^^^^^^^
./specfun.f:3017:1: DO loop currently ends at statement:
  20            F1=PF
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:3120:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:3124:1: DO loop currently ends at statement:
  10            CDY(K)=(1.0D+300,0.0D0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3160:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 30 K=1,K0
                ^^^^^^^^^^^^
./specfun.f:3163:1: DO loop currently ends at statement:
  30               CPZ=CPZ+CRP
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3166:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 35 K=1,K0
                ^^^^^^^^^^^^
./specfun.f:3169:1: DO loop currently ends at statement:
  35               CQZ=CQZ+CRQ
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3208:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 60 K=1,30
                ^^^^^^^^^^^^
./specfun.f:3211:1: DO loop currently ends at statement:
  60               CS0=CS0+CR0*W0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3216:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 65 K=1,30
                ^^^^^^^^^^^^
./specfun.f:3219:1: DO loop currently ends at statement:
  65               CS1=CS1+CR1*(2.0D0*W1+1.0D0/(K+1.0D0))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3243:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 70 K=2,N
             ^^^^^^^^^^^
./specfun.f:3247:1: DO loop currently ends at statement:
  70            CF1=CF
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:3257:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 75 K=M,0,-1
             ^^^^^^^^^^^^^^
./specfun.f:3261:1: DO loop currently ends at statement:
  75            CF1=CF
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:3264:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 80 K=0,N
             ^^^^^^^^^^^
./specfun.f:3265:1: DO loop currently ends at statement:
  80            CBJ(K)=CS*CBJ(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3268:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 85 K=1,N
          ^^^^^^^^^^^
./specfun.f:3269:1: DO loop currently ends at statement:
  85         CDJ(K)=-(K+V0)/Z*CBJ(K)+CBJ(K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3291:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 100 K=LB,1,-1
          ^^^^^^^^^^^^^^^^
./specfun.f:3294:1: DO loop currently ends at statement:
  100        CH1=CH0
  ^^^^^^^^^^^^^^^^^^
./specfun.f:3299:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 105 K=LB,1,-1
          ^^^^^^^^^^^^^^^^
./specfun.f:3302:1: DO loop currently ends at statement:
  105        CH1=CH0
  ^^^^^^^^^^^^^^^^^^
./specfun.f:3315:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 110 K=LB-1,0,-1
          ^^^^^^^^^^^^^^^^^^
./specfun.f:3319:1: DO loop currently ends at statement:
  110        CYL1=CYLK
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:3322:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 115 K=LB+1,N-1
          ^^^^^^^^^^^^^^^^^
./specfun.f:3326:1: DO loop currently ends at statement:
  115        CYL2=CYLK
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:3333:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 130 K=1,N
          ^^^^^^^^^^^^
./specfun.f:3334:1: DO loop currently ends at statement:
  130        CDY(K)=CBY(K-1)-(K+V0)/Z*CBY(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3377:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:3381:1: DO loop currently ends at statement:
  10            CDY(K)=(1.0D+300,0.0D0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3411:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:3414:1: DO loop currently ends at statement:
  25            CPZ=CPZ+CRP
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3417:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:3420:1: DO loop currently ends at statement:
  30            CQZ=CQZ+CRQ
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3448:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 45 K=1,30
                ^^^^^^^^^^^^
./specfun.f:3451:1: DO loop currently ends at statement:
  45               CS0=CS0+CR0*W0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3464:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 50 K=M,0,-1
          ^^^^^^^^^^^^^^
./specfun.f:3468:1: DO loop currently ends at statement:
  50         CF1=CF
  ^^^^^^^^^^^^^^^^^
./specfun.f:3470:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 55 K=0,N
          ^^^^^^^^^^^
./specfun.f:3471:1: DO loop currently ends at statement:
  55         CBJ(K)=CS*CBJ(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3494:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 70 K=1,N
          ^^^^^^^^^^^
./specfun.f:3495:1: DO loop currently ends at statement:
  70         CDJ(K)=-(K+V0)/Z*CBJ(K)+CBJ(K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3497:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 75 K=1,N
          ^^^^^^^^^^^
./specfun.f:3498:1: DO loop currently ends at statement:
  75         CDY(K)=CBY(K-1)-(K+V0)/Z*CBY(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3609:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 45 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:3611:1: DO loop currently ends at statement:
  45            Q0=Q0+B(K)*X**(-2*K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3618:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 50 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:3620:1: DO loop currently ends at statement:
  50            Q1=Q1+B1(K)*X**(-2*K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3790:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=1,10
             ^^^^^^^^^^^^
./specfun.f:3792:1: DO loop currently ends at statement:
  20            TI=TI+A(K)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3803:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=1,50
             ^^^^^^^^^^^^
./specfun.f:3810:1: DO loop currently ends at statement:
  25            TW=TK
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:3815:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=1,10
             ^^^^^^^^^^^^
./specfun.f:3817:1: DO loop currently ends at statement:
  35            TK=TK+A(K)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3854:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:3858:1: DO loop currently ends at statement:
  10            DY(K)=1.0D+300
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3896:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 30 K=1,K0
                ^^^^^^^^^^^^
./specfun.f:3899:1: DO loop currently ends at statement:
  30               PX=PX+RP
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3902:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 35 K=1,K0
                ^^^^^^^^^^^^
./specfun.f:3905:1: DO loop currently ends at statement:
  35               QX=QX+RQ
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3927:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 45 K=2,N
             ^^^^^^^^^^^
./specfun.f:3931:1: DO loop currently ends at statement:
  45            F1=F
  ^^^^^^^^^^^^^^^^^^
./specfun.f:3942:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 50 K=M,0,-1
             ^^^^^^^^^^^^^^
./specfun.f:3946:1: DO loop currently ends at statement:
  50            F1=F
  ^^^^^^^^^^^^^^^^^^
./specfun.f:3952:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 55 K=0,N
             ^^^^^^^^^^^
./specfun.f:3953:1: DO loop currently ends at statement:
  55            BJ(K)=CS*BJ(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3955:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 60 K=2,N
          ^^^^^^^^^^^
./specfun.f:3956:1: DO loop currently ends at statement:
  60         DJ(K)=-(K+V0)/X*BJ(K)+BJ(K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3985:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 80 K=1,30
                ^^^^^^^^^^^^
./specfun.f:3988:1: DO loop currently ends at statement:
  80               CS0=CS0+R0*W0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:3993:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 85 K=1,30
                ^^^^^^^^^^^^
./specfun.f:3996:1: DO loop currently ends at statement:
  85               CS1=CS1+R1*(2.0D0*W1+1.0D0/(K+1.0D0))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4002:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 90 K=2,N
          ^^^^^^^^^^^
./specfun.f:4006:1: DO loop currently ends at statement:
  90         BYV1=BYVK
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:4008:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 95 K=1,N
          ^^^^^^^^^^^
./specfun.f:4009:1: DO loop currently ends at statement:
  95         DY(K)=-(K+V0)/X*BY(K)+BY(K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4039:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:4041:2: DO loop currently ends at statement:
   10           DY(K) = 1.0D+300
   ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4045:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 40 K=1,NM
             ^^^^^^^^^^^^
./specfun.f:4046:2: DO loop currently ends at statement:
   40           DJ(K)=BJ(K-1)-K/X*BJ(K)
   ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4048:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 50 K=1,NM
             ^^^^^^^^^^^^
./specfun.f:4049:2: DO loop currently ends at statement:
   50           DY(K)=BY(K-1)-K*BY(K)/X
   ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4078:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=NMIN,N
             ^^^^^^^^^^^^^^
./specfun.f:4080:1: DO loop currently ends at statement:
  10            BY(K-NMIN)=-1.0D+300
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4099:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=M,0,-1
             ^^^^^^^^^^^^^^
./specfun.f:4109:1: DO loop currently ends at statement:
  15            F1=F
  ^^^^^^^^^^^^^^^^^^
./specfun.f:4111:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=NMIN,NM
             ^^^^^^^^^^^^^^^
./specfun.f:4112:1: DO loop currently ends at statement:
  20            BJ(K-NMIN)=BJ(K-NMIN)/S0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4135:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=1,4
             ^^^^^^^^^^^
./specfun.f:4137:1: DO loop currently ends at statement:
  25            Q0=Q0+B(K)*X**(-2*K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4146:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=1,4
             ^^^^^^^^^^^
./specfun.f:4148:1: DO loop currently ends at statement:
  30            Q1=Q1+B1(K)*X**(-2*K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4153:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=2,NM
             ^^^^^^^^^^^^
./specfun.f:4157:1: DO loop currently ends at statement:
  35            BJ1=BJK
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4161:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 45 K=KY,NM
          ^^^^^^^^^^^^^
./specfun.f:4165:1: DO loop currently ends at statement:
  45         BY1=BYK
  ^^^^^^^^^^^^^^^^^^
./specfun.f:4188:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 45 NR=1,N0
          ^^^^^^^^^^^^^
./specfun.f:4217:1: DO loop currently ends at statement:
  45         W(N+1-NR)=W(NR)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4192:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 I=1,NR-1
             ^^^^^^^^^^^^^^
./specfun.f:4193:1: DO loop currently ends at statement:
  15            P=P*(Z-X(I))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4197:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=2,N
             ^^^^^^^^^^^
./specfun.f:4201:1: DO loop currently ends at statement:
  20            F1=PF
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:4205:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 I=1,NR
             ^^^^^^^^^^^^
./specfun.f:4210:1: DO loop currently ends at statement:
  35            Q=Q+WP
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:4314:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:4318:1: DO loop currently ends at statement:
  10            DY(K)=1.0D+300
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4334:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=2,N
             ^^^^^^^^^^^
./specfun.f:4338:1: DO loop currently ends at statement:
  20            BJ1=BJK
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4349:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=M,0,-1
             ^^^^^^^^^^^^^^
./specfun.f:4353:1: DO loop currently ends at statement:
  30            F1=F
  ^^^^^^^^^^^^^^^^^^
./specfun.f:4359:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 40 K=0,NM
             ^^^^^^^^^^^^
./specfun.f:4360:1: DO loop currently ends at statement:
  40            BJ(K)=CS*BJ(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4362:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 50 K=2,NM
          ^^^^^^^^^^^^
./specfun.f:4363:1: DO loop currently ends at statement:
  50         DJ(K)=BJ(K-1)-K/X*BJ(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4366:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 60 K=2,NM
          ^^^^^^^^^^^^
./specfun.f:4370:1: DO loop currently ends at statement:
  60         F1=F
  ^^^^^^^^^^^^^^^
./specfun.f:4371:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 70 K=2,NM
          ^^^^^^^^^^^^
./specfun.f:4372:1: DO loop currently ends at statement:
  70         DY(K)=BY(K-1)-K*BY(K)/X
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4423:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=2,NA
             ^^^^^^^^^^^^
./specfun.f:4427:1: DO loop currently ends at statement:
  15            PD1=PDF
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4441:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 20 K=2,NA
                ^^^^^^^^^^^^
./specfun.f:4445:1: DO loop currently ends at statement:
  20               PD1=PD
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4455:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 25 K=NK-2,0,-1
                ^^^^^^^^^^^^^^^^^
./specfun.f:4459:1: DO loop currently ends at statement:
  25               F0=F
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4468:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 30 K=M,0,-1
                ^^^^^^^^^^^^^^
./specfun.f:4472:1: DO loop currently ends at statement:
  30               F0=F
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4474:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 35 K=0,NA
                ^^^^^^^^^^^^
./specfun.f:4475:1: DO loop currently ends at statement:
  35               DV(K)=S0*DV(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4531:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=1,20
             ^^^^^^^^^^^^
./specfun.f:4536:1: DO loop currently ends at statement:
  30            A1=AF
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:4539:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=1,10
             ^^^^^^^^^^^^
./specfun.f:4541:1: DO loop currently ends at statement:
  35            BF=BF+A(2*K)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4544:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 40 K=1,10
             ^^^^^^^^^^^^
./specfun.f:4546:1: DO loop currently ends at statement:
  40            BG=BG+A(2*K+1)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4571:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 35 NR=1,NT
          ^^^^^^^^^^^^^
./specfun.f:4597:1: DO loop currently ends at statement:
  35         ZO(NR)=Z
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:4581:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 I=1,NR-1
             ^^^^^^^^^^^^^^
./specfun.f:4582:1: DO loop currently ends at statement:
  20            ZP=ZP*(Z-ZO(I))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4585:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 I=1,NR-1
             ^^^^^^^^^^^^^^
./specfun.f:4591:1: DO loop currently ends at statement:
  30            ZQ=ZQ+ZW
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4621:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 10 K=2,M1
                ^^^^^^^^^^^^
./specfun.f:4622:1: DO loop currently ends at statement:
  10               GA=GA*K
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4631:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 15 K=1,M
                ^^^^^^^^^^^
./specfun.f:4632:1: DO loop currently ends at statement:
  15               R=R*(Z-K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4650:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=25,1,-1
             ^^^^^^^^^^^^^^^
./specfun.f:4651:1: DO loop currently ends at statement:
  20            GR=GR*Z+G(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4760:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:4762:1: DO loop currently ends at statement:
  10            DL(K)=0.0D0
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4769:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=0,N
             ^^^^^^^^^^^
./specfun.f:4778:1: DO loop currently ends at statement:
  25            IF (K.GE.1) DL(K-1)=-0.5D0*X/K*BK
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4800:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 40 K=M,0,-1
          ^^^^^^^^^^^^^^
./specfun.f:4805:1: DO loop currently ends at statement:
  40         F1=F
  ^^^^^^^^^^^^^^^
./specfun.f:4807:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 45 K=0,NM
          ^^^^^^^^^^^^
./specfun.f:4808:1: DO loop currently ends at statement:
  45         BL(K)=BL(K)/BG
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4810:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 50 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:4812:1: DO loop currently ends at statement:
  50         BL(K)=R0*BL(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4814:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 55 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:4815:1: DO loop currently ends at statement:
  55         DL(K)=2.0D0*K/X*(BL(K-1)-BL(K))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4872:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=1,20
             ^^^^^^^^^^^^
./specfun.f:4873:1: DO loop currently ends at statement:
  10            DK(2*K)=K*(B-K)*X/(A+2.0D0*K-1.0D0)/(A+2.0D0*K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4874:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=0,20
             ^^^^^^^^^^^^
./specfun.f:4875:1: DO loop currently ends at statement:
  15            DK(2*K+1)=-(A+K)*(A+B+K)*X/(A+2.D0*K)/(A+2.0*K+1.0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4877:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=20,1,-1
             ^^^^^^^^^^^^^^^
./specfun.f:4878:1: DO loop currently ends at statement:
  20            T1=DK(K)/(1.0D0+T1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4882:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=1,20
             ^^^^^^^^^^^^
./specfun.f:4883:1: DO loop currently ends at statement:
  25            FK(2*K)=K*(A-K)*(1.0D0-X)/(B+2.*K-1.0)/(B+2.0*K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4884:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=0,20
             ^^^^^^^^^^^^
./specfun.f:4885:1: DO loop currently ends at statement:
  30            FK(2*K+1)=-(B+K)*(A+B+K)*(1.D0-X)/
  ^
./specfun.f:4888:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=20,1,-1
             ^^^^^^^^^^^^^^^
./specfun.f:4889:1: DO loop currently ends at statement:
  35            T2=FK(K)/(1.0D0+T2)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4923:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 J=MJ,IC+1,-1
          ^^^^^^^^^^^^^^^^^^
./specfun.f:4924:1: DO loop currently ends at statement:
  10         T1=-Q*Q/((2.0D0*J+L)**2-B+T1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4938:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 J=J0,JF
             ^^^^^^^^^^^^^
./specfun.f:4939:1: DO loop currently ends at statement:
  15            T2=-Q*Q/((2.0D0*J-L-L0)**2-B+T2)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:4972:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 K=2,N
          ^^^^^^^^^^^
./specfun.f:4981:1: DO loop currently ends at statement:
  10         CP1=CPF
  ^^^^^^^^^^^^^^^^^^
./specfun.f:5001:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 K=0,N
          ^^^^^^^^^^^
./specfun.f:5003:1: DO loop currently ends at statement:
  10         QD(K)=0.0D0
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5005:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=0,N
             ^^^^^^^^^^^
./specfun.f:5007:1: DO loop currently ends at statement:
  15            QD(K)=1.0D+300
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5022:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 20 K=2,M
          ^^^^^^^^^^^
./specfun.f:5025:1: DO loop currently ends at statement:
  20         QF1=QM0
  ^^^^^^^^^^^^^^^^^^
./specfun.f:5033:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 25 K=2,N
                ^^^^^^^^^^^
./specfun.f:5037:1: DO loop currently ends at statement:
  25               QF1=QF2
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5041:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=2,M
             ^^^^^^^^^^^
./specfun.f:5044:1: DO loop currently ends at statement:
  30            QG1=QM1
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5051:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 35 K=2,N
                ^^^^^^^^^^^
./specfun.f:5055:1: DO loop currently ends at statement:
  35               QH1=QH2
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5062:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 45 L=2,N
                ^^^^^^^^^^^
./specfun.f:5076:1: DO loop currently ends at statement:
  45               QH1=Q1L
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5067:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 40 K=2,M
                   ^^^^^^^^^^^
./specfun.f:5071:1: DO loop currently ends at statement:
  40                  QF1=QMK
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5086:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 50 K=KM,0,-1
             ^^^^^^^^^^^^^^^
./specfun.f:5090:1: DO loop currently ends at statement:
  50            QF1=QF0
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5091:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 55 K=0,N
             ^^^^^^^^^^^
./specfun.f:5092:1: DO loop currently ends at statement:
  55            QM(K)=QM(K)*QM0/QF0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5095:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 60 K=0,N
             ^^^^^^^^^^^
./specfun.f:5096:1: DO loop currently ends at statement:
  60            QM(K)=(-1)**M*QM(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5099:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 65 K=1,N
          ^^^^^^^^^^^
./specfun.f:5100:1: DO loop currently ends at statement:
  65         QD(K)=(K*X*QM(K)-(K+M)*QM(K-1))/(X*X-1.0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5135:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=1,KM
             ^^^^^^^^^^^^
./specfun.f:5141:1: DO loop currently ends at statement:
  15            CF(K)=CF(K)*CT**K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5139:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 10 I=LF-1,L0,-1
                ^^^^^^^^^^^^^^^^^^
./specfun.f:5140:1: DO loop currently ends at statement:
  10               CF(K)=CF(K)*CT2+A(I)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5144:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=1,KM
             ^^^^^^^^^^^^
./specfun.f:5145:1: DO loop currently ends at statement:
  20            CSI=CSI+CF(K)*VR**K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5149:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=1,KM
             ^^^^^^^^^^^^
./specfun.f:5150:1: DO loop currently ends at statement:
  25            CSK=CSK+(-1)**K*CF(K)*VR**K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5254:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 I=1,10
          ^^^^^^^^^^^^
./specfun.f:5262:1: DO loop currently ends at statement:
  10         EL3=EL3+W(I)*(F1+F2)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5298:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=1,20
             ^^^^^^^^^^^^
./specfun.f:5300:1: DO loop currently ends at statement:
  25            EI=EI+R
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5436:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 20 I=1,8
                ^^^^^^^^^^^
./specfun.f:5440:1: DO loop currently ends at statement:
  20               SUM2=SUM2+R2
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5456:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 I=1,LA-1
             ^^^^^^^^^^^^^^
./specfun.f:5460:1: DO loop currently ends at statement:
  35            A=A+1.0D0
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5523:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=1,NM
             ^^^^^^^^^^^^
./specfun.f:5525:1: DO loop currently ends at statement:
  10            HF=HF+R
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5532:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=1,NM
             ^^^^^^^^^^^^
./specfun.f:5534:1: DO loop currently ends at statement:
  15            HF=HF+R
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5562:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 30 J=1,ABS(M)-1
                ^^^^^^^^^^^^^^^^^^
./specfun.f:5563:1: DO loop currently ends at statement:
  30               GM=GM*J
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5565:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 35 J=1,ABS(M)
                ^^^^^^^^^^^^^^^^
./specfun.f:5566:1: DO loop currently ends at statement:
  35               RM=RM*J
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5575:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 40 K=1,M-1
                   ^^^^^^^^^^^^^
./specfun.f:5577:1: DO loop currently ends at statement:
  40                  F0=F0+R0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5578:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 45 K=1,M
                   ^^^^^^^^^^^
./specfun.f:5579:1: DO loop currently ends at statement:
  45                  SP0=SP0+1.0D0/(A+K-1.0)+1.0/(B+K-1.0)-1.0/K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5581:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 55 K=1,250
                   ^^^^^^^^^^^^^
./specfun.f:5591:1: DO loop currently ends at statement:
  55                  HW=F1
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5584:21: portability: A DO loop should terminate with an END DO or CONTINUE
                      DO 50 J=1,M
                      ^^^^^^^^^^^
./specfun.f:5585:1: DO loop currently ends at statement:
  50                     SM=SM+(1.0D0-A)/((J+K)*(A+J+K-1.0))+1.0/
  ^
./specfun.f:5597:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 65 K=1,M-1
                   ^^^^^^^^^^^^^
./specfun.f:5599:1: DO loop currently ends at statement:
  65                  F0=F0+R0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5600:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 70 K=1,M
                   ^^^^^^^^^^^
./specfun.f:5601:1: DO loop currently ends at statement:
  70                  SP0=SP0+1.0D0/K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5603:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 80 K=1,250
                   ^^^^^^^^^^^^^
./specfun.f:5612:1: DO loop currently ends at statement:
  80                  HW=F1
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5606:21: portability: A DO loop should terminate with an END DO or CONTINUE
                      DO 75 J=1,M
                      ^^^^^^^^^^^
./specfun.f:5607:1: DO loop currently ends at statement:
  75                     SM=SM+1.0D0/(J+K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5628:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 90 K=1,250
                ^^^^^^^^^^^^^
./specfun.f:5634:1: DO loop currently ends at statement:
  90               HW=HF
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5647:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 100 K=1,250
             ^^^^^^^^^^^^^^
./specfun.f:5651:1: DO loop currently ends at statement:
  100           HW=HF
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:5707:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=1,M
             ^^^^^^^^^^^
./specfun.f:5709:1: DO loop currently ends at statement:
  10            CHG=CHG+CR
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5751:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 20 I=1,8
                   ^^^^^^^^^^^
./specfun.f:5755:1: DO loop currently ends at statement:
  20                  CS2=CS2+CR2
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5777:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 35 I=1,LA-1
                ^^^^^^^^^^^^^^
./specfun.f:5781:1: DO loop currently ends at statement:
  35               A=A+1.0D0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5856:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=1,NM
             ^^^^^^^^^^^^
./specfun.f:5858:1: DO loop currently ends at statement:
  10            ZHF=ZHF+ZR
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5864:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=1,NM
             ^^^^^^^^^^^^
./specfun.f:5866:1: DO loop currently ends at statement:
  15            ZHF=ZHF+ZR
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5878:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 20 K=1,500
                ^^^^^^^^^^^^^
./specfun.f:5882:1: DO loop currently ends at statement:
  20               ZW=ZHF
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5897:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 30 J=1,ABS(M)-1
                   ^^^^^^^^^^^^^^^^^^
./specfun.f:5898:1: DO loop currently ends at statement:
  30                  GM=GM*J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5900:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 35 J=1,ABS(M)
                   ^^^^^^^^^^^^^^^^
./specfun.f:5901:1: DO loop currently ends at statement:
  35                  RM=RM*J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5910:21: portability: A DO loop should terminate with an END DO or CONTINUE
                      DO 40 K=1,M-1
                      ^^^^^^^^^^^^^
./specfun.f:5912:1: DO loop currently ends at statement:
  40                     ZF0=ZF0+ZR0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5913:21: portability: A DO loop should terminate with an END DO or CONTINUE
                      DO 45 K=1,M
                      ^^^^^^^^^^^
./specfun.f:5914:1: DO loop currently ends at statement:
  45                     SP0=SP0+1.0D0/(A+K-1.0D0)+1.0/(B+K-1.0D0)-1.D0/K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5916:21: portability: A DO loop should terminate with an END DO or CONTINUE
                      DO 55 K=1,500
                      ^^^^^^^^^^^^^
./specfun.f:5929:1: DO loop currently ends at statement:
  55                     ZW=ZF1
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5935:21: portability: A DO loop should terminate with an END DO or CONTINUE
                      DO 65 K=1,M-1
                      ^^^^^^^^^^^^^
./specfun.f:5938:1: DO loop currently ends at statement:
  65                     ZF0=ZF0+ZR0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5939:21: portability: A DO loop should terminate with an END DO or CONTINUE
                      DO 70 K=1,M
                      ^^^^^^^^^^^
./specfun.f:5940:1: DO loop currently ends at statement:
  70                     SP0=SP0+1.0D0/K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5942:21: portability: A DO loop should terminate with an END DO or CONTINUE
                      DO 80 K=1,500
                      ^^^^^^^^^^^^^
./specfun.f:5952:1: DO loop currently ends at statement:
  80                     ZW=ZF1
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5946:24: portability: A DO loop should terminate with an END DO or CONTINUE
                         DO 75 J=1,M
                         ^^^^^^^^^^^
./specfun.f:5947:1: DO loop currently ends at statement:
  75                        SM=SM+1.0D0/(J+K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5968:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 90 K=1,500
                   ^^^^^^^^^^^^^
./specfun.f:5974:1: DO loop currently ends at statement:
  90                  ZW=ZHF
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:5986:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 100 K=1,1500
                ^^^^^^^^^^^^^^^
./specfun.f:5990:1: DO loop currently ends at statement:
  100              ZW=ZHF
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6009:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 110 K=1,500
                ^^^^^^^^^^^^^^
./specfun.f:6014:1: DO loop currently ends at statement:
  110              ZW=ZHF
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6037:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 120 K=1,MAB-1
                ^^^^^^^^^^^^^^^^
./specfun.f:6042:1: DO loop currently ends at statement:
  120              ZF0=ZF0+ZR*G0/GCBK
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6046:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 125 J=1,MAB
                ^^^^^^^^^^^^^^
./specfun.f:6047:1: DO loop currently ends at statement:
  125              SP=SP+1.0D0/J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6050:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 130 J=1,MAB
                ^^^^^^^^^^^^^^
./specfun.f:6051:1: DO loop currently ends at statement:
  130              SQ=SQ*(B+J-1.0D0)*(B-C+J)/J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6057:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 145 K=1,10000
                ^^^^^^^^^^^^^^^^
./specfun.f:6073:1: DO loop currently ends at statement:
  145              W0=WS
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6061:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 135 J=K+1,K+MAB
                   ^^^^^^^^^^^^^^^^^^
./specfun.f:6062:1: DO loop currently ends at statement:
  135                 RK2=RK2*(B+J-1.0D0)*(B-C+J)/J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6066:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 140 J=K+1,K+MAB
                   ^^^^^^^^^^^^^^^^^^
./specfun.f:6067:1: DO loop currently ends at statement:
  140                 SJ2=SJ2+1.0D0/J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6159:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 35 K=1,16
                ^^^^^^^^^^^^
./specfun.f:6161:1: DO loop currently ends at statement:
  35               SU1=SU1+A(K)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6164:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 40 K=1,16
                ^^^^^^^^^^^^
./specfun.f:6166:1: DO loop currently ends at statement:
  40               SU2=SU2+A(K)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6172:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 45 K=1,8
                ^^^^^^^^^^^
./specfun.f:6174:1: DO loop currently ends at statement:
  45               SU3=SU3+A(2*K)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6177:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 50 K=1,7
                ^^^^^^^^^^^
./specfun.f:6179:1: DO loop currently ends at statement:
  50               SU4=SU4+A(2*K+1)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6213:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:6217:1: DO loop currently ends at statement:
  10            DK(K)=-1.0D+300
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6235:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=2,N
             ^^^^^^^^^^^
./specfun.f:6239:1: DO loop currently ends at statement:
  15         H1=H
  ^^^^^^^^^^^^^^^
./specfun.f:6250:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=M,0,-1
             ^^^^^^^^^^^^^^
./specfun.f:6254:1: DO loop currently ends at statement:
  20            F1=F
  ^^^^^^^^^^^^^^^^^^
./specfun.f:6256:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=0,NM
             ^^^^^^^^^^^^
./specfun.f:6257:1: DO loop currently ends at statement:
  25            BI(K)=S0*BI(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6261:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 30 K=2,NM
          ^^^^^^^^^^^^
./specfun.f:6265:1: DO loop currently ends at statement:
  30         G1=G
  ^^^^^^^^^^^^^^^
./specfun.f:6266:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 40 K=2,NM
          ^^^^^^^^^^^^
./specfun.f:6268:1: DO loop currently ends at statement:
  40         DK(K)=-BK(K-1)-K/X*BK(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6305:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:6309:1: DO loop currently ends at statement:
  10            CDY(K)=(1.0D+300,0.0D0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6327:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=M,0,-1
             ^^^^^^^^^^^^^^
./specfun.f:6341:1: DO loop currently ends at statement:
  15            CF1=CF
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:6347:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=0,NM
             ^^^^^^^^^^^^
./specfun.f:6348:1: DO loop currently ends at statement:
  20            CBJ(K)=CBJ(K)/CS0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6363:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=1,4
             ^^^^^^^^^^^
./specfun.f:6364:1: DO loop currently ends at statement:
  25            CP0=CP0+A(K)*Z**(-2*K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6366:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=1,4
             ^^^^^^^^^^^
./specfun.f:6367:1: DO loop currently ends at statement:
  30            CQ0=CQ0+B(K)*Z**(-2*K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6375:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=1,4
             ^^^^^^^^^^^
./specfun.f:6376:1: DO loop currently ends at statement:
  35            CP1=CP1+A1(K)*Z**(-2*K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6378:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 40 K=1,4
             ^^^^^^^^^^^
./specfun.f:6379:1: DO loop currently ends at statement:
  40            CQ1=CQ1+B1(K)*Z**(-2*K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6384:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 45 K=2,NM
             ^^^^^^^^^^^^
./specfun.f:6388:1: DO loop currently ends at statement:
  45            CBJ1=CBJK
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6391:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 50 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:6392:1: DO loop currently ends at statement:
  50         CDJ(K)=CBJ(K-1)-K/Z*CBJ(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6405:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 60 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:6406:1: DO loop currently ends at statement:
  60         CDY(K)=CBY(K-1)-K/Z*CBY(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6437:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:6441:1: DO loop currently ends at statement:
  10            DK(K)=-1.0D+300
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6458:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 15 K=M,0,-1
          ^^^^^^^^^^^^^^
./specfun.f:6464:1: DO loop currently ends at statement:
  15         F1=F
  ^^^^^^^^^^^^^^^
./specfun.f:6466:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 20 K=0,NM
          ^^^^^^^^^^^^
./specfun.f:6467:1: DO loop currently ends at statement:
  20         BI(K)=S0*BI(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6481:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 25 K=1,K0
                ^^^^^^^^^^^^
./specfun.f:6483:1: DO loop currently ends at statement:
  25               BKL=BKL+R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6489:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 35 K=2,NM
          ^^^^^^^^^^^^
./specfun.f:6493:1: DO loop currently ends at statement:
  35         G1=G
  ^^^^^^^^^^^^^^^
./specfun.f:6496:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 40 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:6498:1: DO loop currently ends at statement:
  40         DK(K)=-BK(K-1)-K/X*BK(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6523:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 I=0,N
          ^^^^^^^^^^^
./specfun.f:6526:1: DO loop currently ends at statement:
  10         PD(J,I)=0.0D0
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6524:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 J=0,M
          ^^^^^^^^^^^
./specfun.f:6526:1: DO loop currently ends at statement:
  10         PD(J,I)=0.0D0
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6530:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 I=1,N
             ^^^^^^^^^^^
./specfun.f:6532:1: DO loop currently ends at statement:
  15            PD(0,I)=0.5D0*I*(I+1.0D0)*X**(I+1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6549:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 30 I=1,M
          ^^^^^^^^^^^
./specfun.f:6550:1: DO loop currently ends at statement:
  30         PM(I,I)=-LS*(2.0D0*I-1.0D0)*XQ*PM(I-1,I-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6551:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 35 I=0,MIN(M,N-1)
          ^^^^^^^^^^^^^^^^^^^^
./specfun.f:6552:1: DO loop currently ends at statement:
  35         PM(I,I+1)=(2.0D0*I+1.0D0)*X*PM(I,I)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6559:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 45 J=1,N
          ^^^^^^^^^^^
./specfun.f:6560:1: DO loop currently ends at statement:
  45         PD(0,J)=LS*J*(PM(0,J-1)-X*PM(0,J))/XS
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6747:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 50 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:6748:1: DO loop currently ends at statement:
  50            CP0=CP0+A(K)*Z1**(-2*K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6750:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 55 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:6751:1: DO loop currently ends at statement:
  55            CQ0=CQ0+B(K)*Z1**(-2*K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6757:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 60 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:6758:1: DO loop currently ends at statement:
  60            CP1=CP1+A1(K)*Z1**(-2*K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6760:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 65 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:6761:1: DO loop currently ends at statement:
  65            CQ1=CQ1+B1(K)*Z1**(-2*K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6852:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 25 K=M,0,-1
                ^^^^^^^^^^^^^^
./specfun.f:6861:1: DO loop currently ends at statement:
  25               XF0=XF
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6869:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 30 K=1,12
                ^^^^^^^^^^^^
./specfun.f:6871:1: DO loop currently ends at statement:
  30               XF=XF+XR
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:6874:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 35 K=1,12
                ^^^^^^^^^^^^
./specfun.f:6876:1: DO loop currently ends at statement:
  35               XG=XG+XR
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7053:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 50 K=1,KMAX
             ^^^^^^^^^^^^^^
./specfun.f:7057:1: DO loop currently ends at statement:
  50            DK(K)=-(6.0D0*K+1.0D0)/(6.0D0*K-1.0D0)*CK(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7062:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 55 K=1,KM
                ^^^^^^^^^^^^
./specfun.f:7065:1: DO loop currently ends at statement:
  55               SAD=SAD+DK(K)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7069:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 60 K=1,KM
                ^^^^^^^^^^^^
./specfun.f:7072:1: DO loop currently ends at statement:
  60               SBD=SBD+DK(K)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7085:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 65 K=1,KM
                ^^^^^^^^^^^^
./specfun.f:7088:1: DO loop currently ends at statement:
  65               SDA=SDA+DK(2*K)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7092:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 70 K=1,KM2
                ^^^^^^^^^^^^^
./specfun.f:7095:1: DO loop currently ends at statement:
  70               SDB=SDB+DK(2*K+1)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7144:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 5 K1=NM,K,-1
                   ^^^^^^^^^^^^^^^
./specfun.f:7145:1: DO loop currently ends at statement:
  5                   CK(K1)=CK(K1)*1.0D-100
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7163:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 10 J=3,KB+1
                   ^^^^^^^^^^^^^^
./specfun.f:7168:1: DO loop currently ends at statement:
  10                  F2=F
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7175:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 25 K=1,KB
          ^^^^^^^^^^^^
./specfun.f:7176:1: DO loop currently ends at statement:
  25         SU1=SU1+CK(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7178:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 30 K=KB+1,NM
          ^^^^^^^^^^^^^^^
./specfun.f:7179:1: DO loop currently ends at statement:
  30         SU2=SU2+CK(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7181:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 35 J=1,(N+M+IP)/2
          ^^^^^^^^^^^^^^^^^^^^
./specfun.f:7182:1: DO loop currently ends at statement:
  35         R1=R1*(J+0.5D0*(N+M+IP))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7184:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 40 J=1,(N-M-IP)/2
          ^^^^^^^^^^^^^^^^^^^^
./specfun.f:7185:1: DO loop currently ends at statement:
  40         R2=-R2*J
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:7191:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 45 K=1,KB
          ^^^^^^^^^^^^
./specfun.f:7192:1: DO loop currently ends at statement:
  45         CK(K)=FL/FS*S0*CK(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7193:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 50 K=KB+1,NM
          ^^^^^^^^^^^^^^^
./specfun.f:7194:1: DO loop currently ends at statement:
  50         CK(K)=S0*CK(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7226:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 35 K=0,NM-1
          ^^^^^^^^^^^^^^
./specfun.f:7247:1: DO loop currently ends at statement:
  35         CK(K+1)=FAC*SUM/R1
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7230:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 I=I1,I1+2*M-1
             ^^^^^^^^^^^^^^^^^^^
./specfun.f:7231:1: DO loop currently ends at statement:
  10            R=R*I
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:7233:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 I=I2,I2+K-1
             ^^^^^^^^^^^^^^^^^
./specfun.f:7234:1: DO loop currently ends at statement:
  15            R=R*(I+0.5D0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7236:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 I=K+1,NM
             ^^^^^^^^^^^^^^
./specfun.f:7243:1: DO loop currently ends at statement:
  20            SW=SUM
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:7245:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 I=2,M+K
             ^^^^^^^^^^^^^
./specfun.f:7246:1: DO loop currently ends at statement:
  30            R1=R1*I
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7303:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 35 NR=1,NT
          ^^^^^^^^^^^^^
./specfun.f:7335:1: DO loop currently ends at statement:
  35         ZO(NR)=Z
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:7319:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 I=1,NR-1
             ^^^^^^^^^^^^^^
./specfun.f:7320:1: DO loop currently ends at statement:
  20            ZP=ZP*(Z-ZO(I))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7323:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 I=1,NR-1
             ^^^^^^^^^^^^^^
./specfun.f:7329:1: DO loop currently ends at statement:
  30            ZQ=ZQ+ZW
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7396:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 J=1,M-1
             ^^^^^^^^^^^^^
./specfun.f:7397:1: DO loop currently ends at statement:
  10            RG=RG*(V*V-J*J)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7400:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 J=1,M
             ^^^^^^^^^^^
./specfun.f:7401:1: DO loop currently ends at statement:
  15            R0=.5D0*R0*XQ/J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7408:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=1,NV-M
             ^^^^^^^^^^^^^^
./specfun.f:7411:1: DO loop currently ends at statement:
  20            PMV=PMV+R
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7431:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 35 J=1,M
                   ^^^^^^^^^^^
./specfun.f:7432:1: DO loop currently ends at statement:
  35                  R2=R2*QR*J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7435:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 40 K=1,M-1
                   ^^^^^^^^^^^^^
./specfun.f:7437:1: DO loop currently ends at statement:
  40                  S0=S0+R1
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7443:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 45 J=1,M
                ^^^^^^^^^^^
./specfun.f:7444:1: DO loop currently ends at statement:
  45               S1=S1+(J*J+V*V)/(J*(J*J-V*V))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7450:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 50 J=1,M
                   ^^^^^^^^^^^
./specfun.f:7451:1: DO loop currently ends at statement:
  50                  S=S+((K+J)**2+V*V)/((K+J)*((K+J)**2-V*V))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7453:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 55 J=1,K
                   ^^^^^^^^^^^
./specfun.f:7454:1: DO loop currently ends at statement:
  55                  S2=S2+1.0D0/(J*(J*J-V*V))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7576:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 K=1,10
          ^^^^^^^^^^^^
./specfun.f:7579:1: DO loop currently ends at statement:
  10         GI=GI-A(K)*T*DSIN((2.0D0*K-1.0D0)*TH)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7583:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 J=0,NA-1
             ^^^^^^^^^^^^^^
./specfun.f:7585:1: DO loop currently ends at statement:
  15            GI1=GI1+DATAN(Y/(X+J))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7644:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:7648:1: DO loop currently ends at statement:
  10         SW=SU1
  ^^^^^^^^^^^^^^^^^
./specfun.f:7651:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 20 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:7655:1: DO loop currently ends at statement:
  20         SW=SU1
  ^^^^^^^^^^^^^^^^^
./specfun.f:7696:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 J=1,150
          ^^^^^^^^^^^^^
./specfun.f:7704:1: DO loop currently ends at statement:
  10         H0=HU
  ^^^^^^^^^^^^^^^^
./specfun.f:7786:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 K=9,1,-1
          ^^^^^^^^^^^^^^
./specfun.f:7787:1: DO loop currently ends at statement:
  10         GL0=GL0*X2+A(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7790:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=1,N
             ^^^^^^^^^^^
./specfun.f:7792:1: DO loop currently ends at statement:
  15            X0=X0-1.0D0
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7825:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=2,N
             ^^^^^^^^^^^
./specfun.f:7830:1: DO loop currently ends at statement:
  15            Q1=QF
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:7942:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:7943:1: DO loop currently ends at statement:
  35            BI0=BI0+A(K)*XR**K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7946:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 40 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:7947:1: DO loop currently ends at statement:
  40            BI1=BI1+B(K)*XR**K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:7956:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 65 K=1,50
             ^^^^^^^^^^^^
./specfun.f:7961:1: DO loop currently ends at statement:
  65            WW=BK0
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:7971:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 75 K=1,8
             ^^^^^^^^^^^
./specfun.f:7972:1: DO loop currently ends at statement:
  75            BK0=BK0+A1(K)*XR2**K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8013:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=2,N
             ^^^^^^^^^^^
./specfun.f:8017:1: DO loop currently ends at statement:
  10            CF1=CF
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:8031:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 15 K=2,N0
                ^^^^^^^^^^^^
./specfun.f:8035:1: DO loop currently ends at statement:
  15               CF1=CF
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8044:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 20 K=NM1,0,-1
                   ^^^^^^^^^^^^^^^^
./specfun.f:8048:1: DO loop currently ends at statement:
  20                  CFA=CF
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8053:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 25 K=M,0,-1
                   ^^^^^^^^^^^^^^
./specfun.f:8057:1: DO loop currently ends at statement:
  25                  CFB=CF
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8059:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 30 K=0,N0
                   ^^^^^^^^^^^^
./specfun.f:8060:1: DO loop currently ends at statement:
  30                  CPB(K)=CS0*CPB(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8066:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=1,N
             ^^^^^^^^^^^
./specfun.f:8067:1: DO loop currently ends at statement:
  35            CPD(K)=-0.5D0*Z*CPB(K)+K*CPB(K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8069:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 40 K=1,N0
             ^^^^^^^^^^^^
./specfun.f:8070:1: DO loop currently ends at statement:
  40            CPD(K)=0.5D0*Z*CPB(K)-CPB(K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8196:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 K=2,N
          ^^^^^^^^^^^
./specfun.f:8205:1: DO loop currently ends at statement:
  10         P1=PF
  ^^^^^^^^^^^^^^^^
./specfun.f:8237:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 5 I=1,251
          ^^^^^^^^^^^^
./specfun.f:8238:1: DO loop currently ends at statement:
  5          FC(I)=0.0D0
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8251:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 7666 I=1,251
                ^^^^^^^^^^^^^^^
./specfun.f:8252:2: DO loop currently ends at statement:
   7666            FC(I)=FNAN
   ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8295:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 7 I=1,251
             ^^^^^^^^^^^^
./specfun.f:8296:2: DO loop currently ends at statement:
   7            FC(I)=FNAN
   ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8347:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=1,KM
             ^^^^^^^^^^^^
./specfun.f:8348:1: DO loop currently ends at statement:
  30            FC(K)=S0*FC(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8366:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 40 K=1,KM
             ^^^^^^^^^^^^
./specfun.f:8367:1: DO loop currently ends at statement:
  40            FC(K)=S0*FC(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8409:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 65 K=1,KM
             ^^^^^^^^^^^^
./specfun.f:8410:1: DO loop currently ends at statement:
  65            FC(K)=S0*FC(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8437:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 90 J=1,KM
             ^^^^^^^^^^^^
./specfun.f:8438:1: DO loop currently ends at statement:
  90            FC(J)=-FC(J)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8466:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:8468:1: DO loop currently ends at statement:
  10            DI(K)=0.0D0
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8486:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=M,0,-1
             ^^^^^^^^^^^^^^
./specfun.f:8490:1: DO loop currently ends at statement:
  15            F1=F
  ^^^^^^^^^^^^^^^^^^
./specfun.f:8492:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=0,NM
             ^^^^^^^^^^^^
./specfun.f:8493:1: DO loop currently ends at statement:
  20            SI(K)=CS*SI(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8496:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 25 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:8497:1: DO loop currently ends at statement:
  25         DI(K)=SI(K-1)-(K+1.0D0)/X*SI(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8542:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 L1=4,200,2
          ^^^^^^^^^^^^^^^^
./specfun.f:8547:1: DO loop currently ends at statement:
  10         H1=HL
  ^^^^^^^^^^^^^^^^
./specfun.f:8569:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 40 L2=5,160,2
          ^^^^^^^^^^^^^^^^
./specfun.f:8574:1: DO loop currently ends at statement:
  40         D2=DL
  ^^^^^^^^^^^^^^^^
./specfun.f:8625:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 J=1,2*M+IP
          ^^^^^^^^^^^^^^^^
./specfun.f:8626:1: DO loop currently ends at statement:
  10         R0=R0*J
  ^^^^^^^^^^^^^^^^^^
./specfun.f:8630:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 15 K=2,NM
          ^^^^^^^^^^^^
./specfun.f:8634:1: DO loop currently ends at statement:
  15         SW=SUC
  ^^^^^^^^^^^^^^^^^
./specfun.f:8640:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 J=1,NM
             ^^^^^^^^^^^^
./specfun.f:8643:1: DO loop currently ends at statement:
  25            SW1=SUM
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8645:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 J=1,(N+M+IP)/2
             ^^^^^^^^^^^^^^^^^^^^
./specfun.f:8646:1: DO loop currently ends at statement:
  35            R1=R1*(J+0.5D0*(N+M+IP))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8648:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 40 J=1,M
             ^^^^^^^^^^^
./specfun.f:8649:1: DO loop currently ends at statement:
  40            R2=2.0D0*C*R2*J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8651:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 45 J=1,(N-M-IP)/2
             ^^^^^^^^^^^^^^^^^^^^
./specfun.f:8652:1: DO loop currently ends at statement:
  45            R3=R3*J
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8670:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 50 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:8682:1: DO loop currently ends at statement:
  50         SW=R1F
  ^^^^^^^^^^^^^^^^^
./specfun.f:8687:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 60 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:8699:1: DO loop currently ends at statement:
  60         SW=SUD
  ^^^^^^^^^^^^^^^^^
./specfun.f:8917:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=1,25
             ^^^^^^^^^^^^
./specfun.f:8922:1: DO loop currently ends at statement:
  15            HU=HU+R
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:8949:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:8952:1: DO loop currently ends at statement:
  10         GW=GF0
  ^^^^^^^^^^^^^^^^^
./specfun.f:8956:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 20 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:8963:1: DO loop currently ends at statement:
  20         GW=GD0
  ^^^^^^^^^^^^^^^^^
./specfun.f:9015:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=1,16
             ^^^^^^^^^^^^
./specfun.f:9020:1: DO loop currently ends at statement:
  30            A1=AF
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:9023:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=1,8
             ^^^^^^^^^^^
./specfun.f:9025:1: DO loop currently ends at statement:
  35            BF=BF+A(2*K)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9028:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 40 K=1,8
             ^^^^^^^^^^^
./specfun.f:9030:1: DO loop currently ends at statement:
  40            BG=BG+A(2*K+1)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9057:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:9059:1: DO loop currently ends at statement:
  10            DY(K)=1.0D+300
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9068:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 15 K=2,N
          ^^^^^^^^^^^
./specfun.f:9073:1: DO loop currently ends at statement:
  15         RF1=RF2
  ^^^^^^^^^^^^^^^^^^
./specfun.f:9076:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 25 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:9077:1: DO loop currently ends at statement:
  25         DY(K)=-K*RY(K)/X+RY(K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9119:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 J=1,N1
             ^^^^^^^^^^^^
./specfun.f:9120:1: DO loop currently ends at statement:
  10            R=(0.5D0/J-1.0D0)*R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9235:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 85 K=1,KM
             ^^^^^^^^^^^^
./specfun.f:9246:1: DO loop currently ends at statement:
  85            QN0=QN0+FAC*RS
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9336:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 15 K=1,150
          ^^^^^^^^^^^^^
./specfun.f:9343:1: DO loop currently ends at statement:
  15         H0=HM1
  ^^^^^^^^^^^^^^^^^
./specfun.f:9350:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 25 M=1,N
          ^^^^^^^^^^^
./specfun.f:9352:1: DO loop currently ends at statement:
  25         IF (B.LT.0.0) S0=S0+(1.0D0-A)/(M*(A+M-1.0D0))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9357:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 50 K=1,150
          ^^^^^^^^^^^^^
./specfun.f:9378:1: DO loop currently ends at statement:
  50         H0=HM2
  ^^^^^^^^^^^^^^^^^
./specfun.f:9361:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 30 M=1,K
                ^^^^^^^^^^^
./specfun.f:9362:1: DO loop currently ends at statement:
  30               S1=S1-(M+2.0D0*A-2.0D0)/(M*(M+A-1.0D0))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9363:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 35 M=1,N
                ^^^^^^^^^^^
./specfun.f:9364:1: DO loop currently ends at statement:
  35               S2=S2+1.0D0/(K+M)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9366:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 40 M=1,K+N
                ^^^^^^^^^^^^^
./specfun.f:9367:1: DO loop currently ends at statement:
  40               S1=S1+(1.0D0-A)/(M*(M+A-1.0D0))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9368:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 45 M=1,K
                ^^^^^^^^^^^
./specfun.f:9369:1: DO loop currently ends at statement:
  45               S2=S2+1.0D0/M
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9387:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 60 K=1,N-1
          ^^^^^^^^^^^^^
./specfun.f:9389:1: DO loop currently ends at statement:
  60         HM3=HM3+R
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:9450:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 I=1,NR-1
             ^^^^^^^^^^^^^^
./specfun.f:9451:1: DO loop currently ends at statement:
  20            ZP=ZP*(Z-ZO(I))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9634:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 J=1,NM
             ^^^^^^^^^^^^
./specfun.f:9637:1: DO loop currently ends at statement:
  10            SW=SUM
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:9685:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 20 K=M,0,-1
          ^^^^^^^^^^^^^^
./specfun.f:9690:1: DO loop currently ends at statement:
  20         F1=F
  ^^^^^^^^^^^^^^^
./specfun.f:9691:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 25 K=0,N
          ^^^^^^^^^^^
./specfun.f:9692:1: DO loop currently ends at statement:
  25         BJ(K+1)=BJ(K+1)/(BS-F)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9695:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 30 K=1,N
          ^^^^^^^^^^^
./specfun.f:9697:1: DO loop currently ends at statement:
  30         FJ(K+1)=(K*K/(X*X)-1.0D0)*BJ(K+1)-DJ(K+1)/X
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9724:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:9726:1: DO loop currently ends at statement:
  10            DJ(K)=0.0D0
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9751:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=M,0,-1
             ^^^^^^^^^^^^^^
./specfun.f:9755:1: DO loop currently ends at statement:
  15            F1=F
  ^^^^^^^^^^^^^^^^^^
./specfun.f:9759:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=0,NM
             ^^^^^^^^^^^^
./specfun.f:9760:1: DO loop currently ends at statement:
  20            SJ(K)=CS*SJ(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9762:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 25 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:9763:1: DO loop currently ends at statement:
  25         DJ(K)=SJ(K-1)-(K+1.0D0)*SJ(K)/X
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9811:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 K=2,N
          ^^^^^^^^^^^
./specfun.f:9826:1: DO loop currently ends at statement:
  10         DY1=DYN
  ^^^^^^^^^^^^^^^^^^
./specfun.f:9859:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 15 M=1,NT
          ^^^^^^^^^^^^
./specfun.f:9887:1: DO loop currently ends at statement:
  15         RT=RT+4.44D0
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9975:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,NM
             ^^^^^^^^^^^^
./specfun.f:9977:1: DO loop currently ends at statement:
  10            CHD1(K)=CDJ(K)+CI*CDY(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9981:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=0,NM
             ^^^^^^^^^^^^
./specfun.f:9984:1: DO loop currently ends at statement:
  15            CFAC=CFAC*CI
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9990:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=0,NM
             ^^^^^^^^^^^^
./specfun.f:9993:1: DO loop currently ends at statement:
  20            CFAC=CFAC*CF1
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:9995:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=0,NM
             ^^^^^^^^^^^^
./specfun.f:9997:1: DO loop currently ends at statement:
  25            CHD2(K)=CDJ(K)-CI*CDY(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10000:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=0,NM
             ^^^^^^^^^^^^
./specfun.f:10004:1: DO loop currently ends at statement:
  30            CHD2(K)=CDJ(K)-CI*CDY(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10172:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:10176:1: DO loop currently ends at statement:
  10            DK(K)=1.0D+300
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10209:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:10211:1: DO loop currently ends at statement:
  25            SUM=SUM+R
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10224:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 30 K=M,0,-1
          ^^^^^^^^^^^^^^
./specfun.f:10228:1: DO loop currently ends at statement:
  30         F1=F
  ^^^^^^^^^^^^^^^
./specfun.f:10230:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 35 K=0,N
          ^^^^^^^^^^^
./specfun.f:10231:1: DO loop currently ends at statement:
  35         BI(K)=CS*BI(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10233:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 40 K=1,N
          ^^^^^^^^^^^
./specfun.f:10234:1: DO loop currently ends at statement:
  40         DI(K)=-(K+V0)/X*BI(K)+BI(K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10241:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 45 K=1,50
                ^^^^^^^^^^^^
./specfun.f:10247:1: DO loop currently ends at statement:
  45               WW=WA
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10257:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 55 K=1,120
                ^^^^^^^^^^^^^
./specfun.f:10263:1: DO loop currently ends at statement:
  55               WW=WA
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10270:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 65 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:10272:1: DO loop currently ends at statement:
  65            SUM=SUM+R
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10278:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 70 K=2,N
          ^^^^^^^^^^^
./specfun.f:10282:1: DO loop currently ends at statement:
  70         BK1=BK2
  ^^^^^^^^^^^^^^^^^^
./specfun.f:10284:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 80 K=1,N
          ^^^^^^^^^^^
./specfun.f:10285:1: DO loop currently ends at statement:
  80         DK(K)=-(K+V0)/X*BK(K)-BK(K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10315:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 5 I=1,NM
             ^^^^^^^^^^^
./specfun.f:10316:1: DO loop currently ends at statement:
  5             DF(I)=0D0
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10349:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 12 K1=K,NM
                   ^^^^^^^^^^^^^
./specfun.f:10350:1: DO loop currently ends at statement:
  12                  DF(K1)=DF(K1)*1.0D-100
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10367:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 20 J=3,KB+1
                   ^^^^^^^^^^^^^^
./specfun.f:10377:1: DO loop currently ends at statement:
  20                  F2=F
  ^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10371:24: portability: A DO loop should terminate with an END DO or CONTINUE
                         DO 15 K1=1,J
                         ^^^^^^^^^^^^
./specfun.f:10372:1: DO loop currently ends at statement:
  15                        DF(K1)=DF(K1)*1.0D-100
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10385:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 40 J=M+IP+1,2*(M+IP)
          ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10386:1: DO loop currently ends at statement:
  40         R1=R1*J
  ^^^^^^^^^^^^^^^^^^
./specfun.f:10388:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 45 K=2,KB
          ^^^^^^^^^^^^
./specfun.f:10390:1: DO loop currently ends at statement:
  45           SU1=SU1+R1*DF(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10393:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 50 K=KB+1,NM
          ^^^^^^^^^^^^^^^
./specfun.f:10397:1: DO loop currently ends at statement:
  50         SW=SU2
  ^^^^^^^^^^^^^^^^^
./specfun.f:10399:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 60 J=1,(M+N+IP)/2
          ^^^^^^^^^^^^^^^^^^^^
./specfun.f:10400:1: DO loop currently ends at statement:
  60         R3=R3*(J+0.5D0*(N+M+IP))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10402:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 65 J=1,(N-M-IP)/2
          ^^^^^^^^^^^^^^^^^^^^
./specfun.f:10403:1: DO loop currently ends at statement:
  65         R4=-4.0D0*R4*J
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10405:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 70 K=1,KB
          ^^^^^^^^^^^^
./specfun.f:10406:1: DO loop currently ends at statement:
  70         DF(K)=FL/FS*S0*DF(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10407:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 75 K=KB+1,NM
          ^^^^^^^^^^^^^^^
./specfun.f:10408:1: DO loop currently ends at statement:
  75         DF(K)=S0*DF(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10478:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 30 K=1,K0
                ^^^^^^^^^^^^
./specfun.f:10481:1: DO loop currently ends at statement:
  30               PX=PX+RP
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10484:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 35 K=1,K0
                ^^^^^^^^^^^^
./specfun.f:10487:1: DO loop currently ends at statement:
  35               QX=QX+RQ
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10543:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 75 K=1,K0
                ^^^^^^^^^^^^
./specfun.f:10545:1: DO loop currently ends at statement:
  75               VSL=VSL+R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10573:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 100 K=1,K0
                ^^^^^^^^^^^^^
./specfun.f:10575:1: DO loop currently ends at statement:
  100              SUM=SUM+R
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10622:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:10626:1: DO loop currently ends at statement:
  10            CDK(K)=1.0D+300
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10658:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:10660:1: DO loop currently ends at statement:
  25            CS=CS+CR
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10671:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 30 K=M,0,-1
          ^^^^^^^^^^^^^^
./specfun.f:10675:1: DO loop currently ends at statement:
  30         CF1=CF
  ^^^^^^^^^^^^^^^^^
./specfun.f:10677:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 35 K=0,N
          ^^^^^^^^^^^
./specfun.f:10678:1: DO loop currently ends at statement:
  35         CBI(K)=CS*CBI(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10714:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 60 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:10716:1: DO loop currently ends at statement:
  60            CS=CS+CR
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10738:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 80 K=1,N
          ^^^^^^^^^^^
./specfun.f:10740:1: DO loop currently ends at statement:
  80         CDK(K)=-(K+V0)/Z*CBK(K)-CBK(K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10784:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:10788:1: DO loop currently ends at statement:
  10            CDK(K)=1.0D+300
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10820:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:10822:1: DO loop currently ends at statement:
  25            CS=CS+CR
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10833:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 30 K=M,0,-1
          ^^^^^^^^^^^^^^
./specfun.f:10837:1: DO loop currently ends at statement:
  30         CF1=CF
  ^^^^^^^^^^^^^^^^^
./specfun.f:10839:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 35 K=0,N
          ^^^^^^^^^^^
./specfun.f:10840:1: DO loop currently ends at statement:
  35         CBI(K)=CS*CBI(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10878:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 60 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:10880:1: DO loop currently ends at statement:
  60            CS=CS+CR
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10888:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 65 K=2,N
          ^^^^^^^^^^^
./specfun.f:10892:1: DO loop currently ends at statement:
  65         CG1=CGK
  ^^^^^^^^^^^^^^^^^^
./specfun.f:10907:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 75 K=1,N
          ^^^^^^^^^^^
./specfun.f:10909:1: DO loop currently ends at statement:
  75         CDK(K)=-(K+V0)/Z*CBK(K)-CBK(K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10941:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=1,80
             ^^^^^^^^^^^^
./specfun.f:10947:1: DO loop currently ends at statement:
  10            WA0=WA
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:10953:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=M,0,-1
             ^^^^^^^^^^^^^^
./specfun.f:10957:1: DO loop currently ends at statement:
  15            CF0=CF
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:10976:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=1,20
             ^^^^^^^^^^^^
./specfun.f:10978:1: DO loop currently ends at statement:
  20            CF=CF+CR
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:10981:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=1,12
             ^^^^^^^^^^^^
./specfun.f:10983:1: DO loop currently ends at statement:
  25            CG=CG+CR
  ^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11038:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=M,0,-1
             ^^^^^^^^^^^^^^
./specfun.f:11047:1: DO loop currently ends at statement:
  25            F0=F
  ^^^^^^^^^^^^^^^^^^
./specfun.f:11054:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=1,20
             ^^^^^^^^^^^^
./specfun.f:11056:1: DO loop currently ends at statement:
  30            F=F+R
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:11059:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=1,12
             ^^^^^^^^^^^^
./specfun.f:11061:1: DO loop currently ends at statement:
  35            G=G+R
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:11094:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:11096:1: DO loop currently ends at statement:
  10            DJ(K)=0.0D0
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11115:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=M,0,-1
             ^^^^^^^^^^^^^^
./specfun.f:11119:1: DO loop currently ends at statement:
  15            F1=F
  ^^^^^^^^^^^^^^^^^^
./specfun.f:11122:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 K=0,NM
             ^^^^^^^^^^^^
./specfun.f:11123:1: DO loop currently ends at statement:
  20            RJ(K)=CS*RJ(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11126:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 25 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:11127:1: DO loop currently ends at statement:
  25         DJ(K)=-K*RJ(K)/X+RJ(K-1)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11154:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 40 NR=1,N/2
          ^^^^^^^^^^^^^^
./specfun.f:11188:1: DO loop currently ends at statement:
  40         W(N+1-NR)=W(NR)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11162:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=2,N
             ^^^^^^^^^^^
./specfun.f:11166:1: DO loop currently ends at statement:
  15            F1=HF
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:11168:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 20 I=1,NR-1
             ^^^^^^^^^^^^^^
./specfun.f:11169:1: DO loop currently ends at statement:
  20            P=P*(Z-X(I))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11172:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 I=1,NR-1
             ^^^^^^^^^^^^^^
./specfun.f:11178:1: DO loop currently ends at statement:
  30            Q=Q+WP
  ^^^^^^^^^^^^^^^^^^^^
./specfun.f:11185:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=1,N
             ^^^^^^^^^^^
./specfun.f:11186:1: DO loop currently ends at statement:
  35            R=2.0D0*R*K
  ^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11297:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=2,N
             ^^^^^^^^^^^
./specfun.f:11298:1: DO loop currently ends at statement:
  10            EN(K)=1.0D0/(K-1.0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11305:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 15 J=1,L-1
                ^^^^^^^^^^^^^
./specfun.f:11306:1: DO loop currently ends at statement:
  15               RP=-RP*X/J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11308:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 20 M=1,L-1
                ^^^^^^^^^^^^^
./specfun.f:11309:1: DO loop currently ends at statement:
  20               PS=PS+1.0D0/M
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11315:18: portability: A DO loop should terminate with an END DO or CONTINUE
                   DO 25 J=1,M
                   ^^^^^^^^^^^
./specfun.f:11316:1: DO loop currently ends at statement:
  25                  R=-R*X/J
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11326:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 50 L=1,N
             ^^^^^^^^^^^
./specfun.f:11331:1: DO loop currently ends at statement:
  50            EN(L)=EXP(-X)*T
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11328:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 45 K=M,1,-1
                ^^^^^^^^^^^^^^
./specfun.f:11329:1: DO loop currently ends at statement:
  45               T0=(L+K-1.0D0)/(1.0D0+K/(X+T0))
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11354:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:11356:1: DO loop currently ends at statement:
  10            DK(K)=-1.0D+300
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11363:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 15 K=2,N
          ^^^^^^^^^^^
./specfun.f:11368:1: DO loop currently ends at statement:
  15         F1=F
  ^^^^^^^^^^^^^^^
./specfun.f:11371:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 25 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:11372:1: DO loop currently ends at statement:
  25         DK(K)=-SK(K-1)-(K+1.0D0)/X*SK(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11393:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 K=2,N
          ^^^^^^^^^^^
./specfun.f:11396:1: DO loop currently ends at statement:
  10         E1=EK
  ^^^^^^^^^^^^^^^^
./specfun.f:11417:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=2,M1
             ^^^^^^^^^^^^
./specfun.f:11418:1: DO loop currently ends at statement:
  10            GA=GA*K
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11422:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=1,M
             ^^^^^^^^^^^
./specfun.f:11423:1: DO loop currently ends at statement:
  15            GA=0.5D0*GA*(2.0D0*K-1.0D0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11483:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=KV,NA
             ^^^^^^^^^^^^^
./specfun.f:11487:1: DO loop currently ends at statement:
  15            F1=F
  ^^^^^^^^^^^^^^^^^^
./specfun.f:11498:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 20 K=KV-2,0,-1
                ^^^^^^^^^^^^^^^^^
./specfun.f:11502:1: DO loop currently ends at statement:
  20               F0=F
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11509:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 25 K=M,0,-1
                ^^^^^^^^^^^^^^
./specfun.f:11513:1: DO loop currently ends at statement:
  25               F0=F
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11515:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 30 K=0,NA
                ^^^^^^^^^^^^
./specfun.f:11516:1: DO loop currently ends at statement:
  30               VV(K)=S0*VV(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11529:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 35 K=2,NA
                ^^^^^^^^^^^^
./specfun.f:11533:1: DO loop currently ends at statement:
  35               F1=F
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11612:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 K=KM,0,-1
             ^^^^^^^^^^^^^^^
./specfun.f:11616:1: DO loop currently ends at statement:
  25            CQF1=CQF0
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11617:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=0,N
             ^^^^^^^^^^^
./specfun.f:11618:1: DO loop currently ends at statement:
  30            CQM(0,K)=CQ0*CQM(0,K)/CQF0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11621:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=KM,0,-1
             ^^^^^^^^^^^^^^^
./specfun.f:11625:1: DO loop currently ends at statement:
  35            CQF1=CQF0
  ^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11627:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 40 K=0,N
             ^^^^^^^^^^^
./specfun.f:11628:1: DO loop currently ends at statement:
  40            CQM(1,K)=CQ10*CQM(1,K)/CQF0
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11640:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 50 J=1,N
          ^^^^^^^^^^^
./specfun.f:11641:1: DO loop currently ends at statement:
  50         CQD(0,J)=LS*J*(CQM(0,J-1)-Z*CQM(0,J))/ZS
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11672:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 5 I=1,N-M+1
             ^^^^^^^^^^^^^^
./specfun.f:11673:1: DO loop currently ends at statement:
  5             EG(I)=(I+M)*(I+M-1.0D0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11681:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 I=1,NM
             ^^^^^^^^^^^^
./specfun.f:11691:1: DO loop currently ends at statement:
  10            G(I)=K*(K-1.0)/((DK2-3.0)*(DK2-1.0))*CS
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11692:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 15 K=2,NM
             ^^^^^^^^^^^^
./specfun.f:11694:1: DO loop currently ends at statement:
  15            F(K)=E(K)*E(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11707:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 25 I=1,ICM
             ^^^^^^^^^^^^^
./specfun.f:11709:1: DO loop currently ends at statement:
  25            H(I)=XB
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11781:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:11785:1: DO loop currently ends at statement:
  10            CDK(K)=-(1.0D+300,0.0D0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11804:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 15 K=M,0,-1
          ^^^^^^^^^^^^^^
./specfun.f:11810:1: DO loop currently ends at statement:
  15         CF1=CF
  ^^^^^^^^^^^^^^^^^
./specfun.f:11812:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 20 K=0,NM
          ^^^^^^^^^^^^
./specfun.f:11813:1: DO loop currently ends at statement:
  20         CBI(K)=CS0*CBI(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11827:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 25 K=1,K0
                ^^^^^^^^^^^^
./specfun.f:11829:1: DO loop currently ends at statement:
  25               CBKL=CBKL+CR
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11835:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 35 K=2,NM
          ^^^^^^^^^^^^
./specfun.f:11839:1: DO loop currently ends at statement:
  35         CG1=CG
  ^^^^^^^^^^^^^^^^^
./specfun.f:11854:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 50 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:11856:1: DO loop currently ends at statement:
  50         CDK(K)=-CBK(K-1)-K/Z*CBK(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11888:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:11892:1: DO loop currently ends at statement:
  10            CDK(K)=(1.0D+300,0.0D0)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11915:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 45 K=M,0,-1
          ^^^^^^^^^^^^^^
./specfun.f:11919:1: DO loop currently ends at statement:
  45         CF1=CF
  ^^^^^^^^^^^^^^^^^
./specfun.f:11921:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 50 K=0,NM
          ^^^^^^^^^^^^
./specfun.f:11922:1: DO loop currently ends at statement:
  50         CBI(K)=CS*CBI(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11923:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 60 K=2,NM
          ^^^^^^^^^^^^
./specfun.f:11929:1: DO loop currently ends at statement:
  60         CBK(K)=CKK
  ^^^^^^^^^^^^^^^^^^^^^
./specfun.f:11930:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 70 K=2,NM
          ^^^^^^^^^^^^
./specfun.f:11932:1: DO loop currently ends at statement:
  70         CDK(K)=-CBK(K-1)-K/Z*CBK(K)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:12004:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 30 K=1,KM
          ^^^^^^^^^^^^
./specfun.f:12015:1: DO loop currently ends at statement:
  30         W1=F1R
  ^^^^^^^^^^^^^^^^^
./specfun.f:12018:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 40 K=1,KM
          ^^^^^^^^^^^^
./specfun.f:12032:1: DO loop currently ends at statement:
  40         W2=D1R
  ^^^^^^^^^^^^^^^^^
./specfun.f:12036:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 55 K=1,KM
          ^^^^^^^^^^^^
./specfun.f:12047:1: DO loop currently ends at statement:
  55         W1=F2R
  ^^^^^^^^^^^^^^^^^
./specfun.f:12050:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 65 K=1,KM
          ^^^^^^^^^^^^
./specfun.f:12064:1: DO loop currently ends at statement:
  65         W2=D2R
  ^^^^^^^^^^^^^^^^^
./specfun.f:12146:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 30 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:12147:1: DO loop currently ends at statement:
  30            CBI0=CBI0+A(K)*ZR**K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:12150:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 35 K=1,K0
             ^^^^^^^^^^^^
./specfun.f:12151:1: DO loop currently ends at statement:
  35            CBI1=CBI1+B(K)*ZR**K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:12159:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 40 K=1,50
             ^^^^^^^^^^^^
./specfun.f:12164:1: DO loop currently ends at statement:
  40            CW=CS
  ^^^^^^^^^^^^^^^^^^^
./specfun.f:12175:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 50 K=1,10
             ^^^^^^^^^^^^
./specfun.f:12176:1: DO loop currently ends at statement:
  50            CBK0=CBK0+A1(K)*ZR2**K
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:12237:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=1,8
             ^^^^^^^^^^^
./specfun.f:12239:1: DO loop currently ends at statement:
  10            PSI=PSI-A(K)*Z2**(-K)*DSIN(2.0D0*K*TH)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:12243:15: portability: A DO loop should terminate with an END DO or CONTINUE
                DO 20 K=1,N
                ^^^^^^^^^^^
./specfun.f:12245:1: DO loop currently ends at statement:
  20               RI=RI+Y/((X0-K)**2.0D0+Y*Y)
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:12280:12: portability: A DO loop should terminate with an END DO or CONTINUE
             DO 10 K=0,N
             ^^^^^^^^^^^
./specfun.f:12282:1: DO loop currently ends at statement:
  10            DY(K)=1.0D+300
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:12293:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 15 K=2,N
          ^^^^^^^^^^^
./specfun.f:12298:1: DO loop currently ends at statement:
  15         F1=F
  ^^^^^^^^^^^^^^^
./specfun.f:12300:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 25 K=1,NM
          ^^^^^^^^^^^^
./specfun.f:12301:1: DO loop currently ends at statement:
  25         DY(K)=SY(K-1)-(K+1.0D0)*SY(K)/X
  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
./specfun.f:12325:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 10 N=1,40
          ^^^^^^^^^^^^
./specfun.f:12332:1: DO loop currently ends at statement:
  10         B0=B
  ^^^^^^^^^^^^^^^
./specfun.f:12335:9: portability: A DO loop should terminate with an END DO or CONTINUE
          DO 20 J=N,1,-1
          ^^^^^^^^^^^^^^
./specfun.f:12339:1: DO loop currently ends at statement:
  20         DN=D
  ^^^^^^^^^^^^^^^
/usr/local/clang18/bin/clang -shared -L/usr/local/clang/lib64 -L/usr/local/clang18/lib -L/usr/local/clang18/lib/x86_64-unknown-linux-gnu -L/usr/local/gcc13/lib64 -L/usr/local/lib64 -o GeoModels.so CompositeLikelihood2.o CompositeLikelihood2_ani.o CompositeLikelihoodCond2.o CompositeLikelihoodCond2_ani.o CorrelationFunction.o Distributions.o GeoModels_init.o Godambe.o Gradient.o HypergeometricFunctions.o TB.o Utility.o mvn.o mvndstpack.o specfun.o weightedleastsquare.o -L/data/gannet/ripley/R/R-clang/lib -lRlapack -L/data/gannet/ripley/R/R-clang/lib -lRblas -lFortran_main -lFortranRuntime -lFortranDecimal -lm -lFortran_main -lFortranRuntime -lFortranDecimal -lm
make[1]: Leaving directory '/data/gannet/ripley/R/packages/tests-clang/GeoModels/src'
make[1]: Entering directory '/data/gannet/ripley/R/packages/tests-clang/GeoModels/src'
make[1]: Leaving directory '/data/gannet/ripley/R/packages/tests-clang/GeoModels/src'
installing to /data/gannet/ripley/R/packages/tests-clang/GeoModels.Rcheck/00LOCK-GeoModels/00new/GeoModels/libs
** R
** data
** byte-compile and prepare package for lazy loading
** help
*** installing help indices
** building package indices
** testing if installed package can be loaded from temporary location
** checking absolute paths in shared objects and dynamic libraries
** testing if installed package can be loaded from final location
** testing if installed package keeps a record of temporary installation path
* DONE (GeoModels)