新手关于用fortran进行mk检验的问题

张俊俊

之前要用fortran进行mk检验，然后就下了个mk检验的代码的安装包。安装包打开是这样的 打开第二个程序代码是
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=?,NYEAR=?')
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=' ')
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=' ',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
这样的，不知道有没有错，然后是要做65年从1951年开始的改了N=65,NYEAR=1951，然后OPEN(2,FILE=' ')这里路径是有要求的吗，是要建专门的txt文本还是怎么回事？反正建了两个文本，将数据输入到了第一个文本里，就是1.txt里，然后输入的是OPEN(2,FILE='1.txt ')，然后OPEN(3,FILE=' ')这里输入的是OPEN(2,FILE='mk.txt ')，点了f5运行，然后出现了 输入65 1951 出现这个，意思是不是我的路径出了问题？该怎么解决呢？我不是专门搞这个的，也从来没接触过这个，就是临时抱佛脚，所以有很多不理解的希望大家不要喷我，主要比较急，也没太多时间研究，希望各位大神能帮忙解答

daqidongq

通道2的1.text是什么？难道不是1.txt?提示这个错误还有可能是你读取该文件的行列数有问题

张俊俊

daqidongq 发表于 2017-7-21 08:54
通道2的1.text是什么？难道不是1.txt?提示这个错误还有可能是你读取该文件的行列数有问题

txt那个是打错了，行列数是提取文件里的行列数吗？就是说每一行的数据有数量上的要求还是什么，比如不能超过多少个，还是横纵有什么要求，刚用，不太懂

daqidongq

张俊俊 发表于 2017-7-21 20:47
txt那个是打错了，行列数是提取文件里的行列数吗？就是说每一行的数据有数量上的要求还是什么，比如不能 ...

就是你要提取的文件里的行列数