关于MK检验 fortran程序的问题

我是潘潘潘

在气象家园下载的好心人提供的程序：
C        THIS IS A PROGRAM FOR DETECTING ABRUPT CLIMATIC CHANGE
C        BY USING MANN-KENDALL TECHNIQUE  
        PROGRAM MK
        DIMENSION Y(1000),YY(1000),U(1000),UF(1000),UB(1000),
     &        M(1000),MD(1000)

        N=39
          NYEAR=1969
C        ***************************************************
C        * N:     SAMPLE SIZE                              *
C        * NYEAR: FIRST YEAR OF THE TIME SERIES            *
C        * Y(N):  ORIGINAL TIME SERIES                     *
C        * UF(N): ORIGINAL SERIES OF U(LN)                 *
C        * UB(N): COUNTER SERIES OF U(LN)                  *
C        * A,B:   CRITICAL VALUE 1.96 AND -1.96            *
C        ***************************************************            
        OPEN(2,FILE='c:\mk\shuju year.txt')
        READ(2,*)(Y(I),I=1,N)
         do i=1,N
          Y(i)=Y(i)/10.0
         enddo
        CALL SMK(Y,M,MD,UF,N)
        DO 20 I=1,N
  20        YY(I)=Y(N+1-I)
        CALL SMK(YY,M,MD,U,N)
        DO 30 I=1,N
  30        UB(I)=-U(N+1-I)
        OPEN(3,FILE='c:\mk\r-mk.grd',form='binary')
        A=1.96
        B=-1.96
        DO 40 I=1,N
        WRITE(3)UF(I),UB(I),A,B
  50        FORMAT(1X,I4,4F8.2)
  40        CONTINUE
        CLOSE(3)
        STOP
        END
C***********************************************************
        SUBROUTINE SMK(Y,M,MD,U,N)
        DIMENSION Y(N),M(N),MD(N),U(N)
        M(1)=0
        DO 10 I=2,N
        M(I)=0
        MD(I)=0
        DO 20 J=1,I-1
        IF(Y(I).LT.Y(J))GOTO 20
        M(I)=M(I)+1
  20        CONTINUE
        MD(I)=MD(I-1)+M(I)
  10        CONTINUE
        U(1)=0.0
        DO 30 I=2,N
        E=I*(I-1)/4.00
        VAR=I*(I-1)*(2*I+5)/72.00
        U(I)=(MD(I)-E)/SQRT(VAR)
  30        CONTINUE
        RETURN
        END


在我改了里面的文件名称 运行以后是这样的：

forrtl: severe (59): list-directed I/O syntax error, unit 2, file c:\mk\shuju ye
ar.txt
Image              PC        Routine            Line        Source
m-k.exe            0040BA39  Unknown               Unknown  Unknown

Incrementally linked image--PC correlation disabled.
Press any key to continue


想问下 这个是哪里出错导致的呢。。。。求

小心小安静

我也是这样的错误，想问楼主找到原因了吗？

我是潘潘潘

小心小安静 发表于 2018-4-21 17:41
我也是这样的错误，想问楼主找到原因了吗？

找到了！是自己的txt数据文件格式有问题 论坛里面有人发过检验用的txt 你按照那个的格式改自己的数据吧 程序本身没有问题的

笨拙的土豆

格式就是要竖着排列吧，但是我的还是不行啊，奇怪

地表植被

笨拙的土豆 发表于 2018-11-14 15:58
格式就是要竖着排列吧，但是我的还是不行啊，奇怪

是不是初始值没改

地表植被

请问楼主你是用什么画的图

kevin520kaka

这个程序怎么做的呀？求大神具体指导一下

Datian

kevin520kaka 发表于 2019-8-20 17:14
这个程序怎么做的呀？求大神具体指导一下

直接运行就完了