魏凤英 mk检验程序 不能运行

xiaohundanpipi

从论坛里下了魏凤英老师的mk检验程序，但是不管是自己的数据还是坛子里的测试数据，运行之后只能到这一步，就不能往下运行了，有没有人知道是为什么？是fortran 版本的问题吗?这个程序好像是fortran77的，我的软件是fortran6.0.
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)
        WRITE(*,10)
  10        FORMAT(2X,'N=91,NYEAR=1860')
        READ(*,*)N,NYEAR
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='1.txt')
        READ(2,*)(Y(I),I=1,N)
        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='mk.txt',STATUS='new')
        A=1.96
        B=-1.96
        DO 40 I=1,N
        WRITE(3,50)NYEAR+I-1,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


xiaohundanpipi

http://bbs.06climate.com/forum.php?mod=viewthread&tid=12989 后来用了晴天的这个程序就可以运行了，但是我还是不知道为什么。

ab3720587

我也在用这个方法，到时候我试试这程序看看是什么问题

言深深

ab3720587 发表于 2013-5-15 11:03

                               
登录/注册后可看大图

我也在用这个方法，到时候我试试这程序看看是什么问题

无语···这是等待输入的部分，你输入“91，1860”按回车就可以了···

善人/jw

LZ肯定没好好看程序的具体，按回车

xiaohundanpipi

言深深 发表于 2013-5-15 11:49

                               
登录/注册后可看大图

无语···这是等待输入的部分，你输入“91，1860”按回车就可以了···

原来是这样... 没仔细看程序就往里面带，好惭愧，谢谢言深深

xiaohundanpipi

善人/jw 发表于 2013-5-15 11:58

                               
登录/注册后可看大图

LZ肯定没好好看程序的具体，按回车

嗯嗯...果然是这样，好惭愧，下次注意

for_apply

我也在用这个方法，但是也搞不懂

xiaohundanpipi

for_apply 发表于 2013-5-22 15:01

                               
登录/注册后可看大图

我也在用这个方法，但是也搞不懂

怎么搞不懂，先把魏凤英老师的书下下来看看，上面MK检验部分说的挺清楚的，论坛里就有下载的。

for_apply

xiaohundanpipi 发表于 2013-5-25 10:50

                               
登录/注册后可看大图

怎么搞不懂，先把魏凤英老师的书下下来看看，上面MK检验部分说的挺清楚的，论坛里就有下载的。

谢谢，我已经弄出来了，谢谢指点