关于M-K检验遇到问题

泷岛彗

本人第一次用Fortran，不太明白，按论坛里前辈们的说法在读入数据中选择了自己的文件，为什么运行后
没有结果呢，而是出现
[img]file:///C:/Documents%20and%20Settings/Administrator/Application%20Data/Tencent/Users/234391930/QQ/WinTemp/RichOle/HBS)]$ZM$~SN%7B~[BOT)@JGT.jpg[/img]
请高手们指导下菜鸟 急。。。
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=49,NYEAR=1956')
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='F:\MK\test.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='F:\MK\SHUCHU ',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


泷岛彗

这是运行后出现的 ，不太明白啊

G-mimzy

结果已经输出到 F:\MK\SHUCHU 文件里面了

泷岛彗

G-mimzy 发表于 2012-12-30 14:16

                               
登录/注册后可看大图

结果已经输出到 F:\MK\SHUCHU 文件里面了

你好 问题是没有这个SHUCHU文件啊？？？

G-mimzy

改成 F:\MK\SHUCHU.txt 试试

泷岛彗

G-mimzy 发表于 2012-12-30 14:34

                               
登录/注册后可看大图

改成 F:\MK\SHUCHU.txt 试试

还是没有 到底怎么回事啊 弄了快一小天了。。。

爱的侍者

泷岛彗 发表于 2012-12-30 14:39

                               
登录/注册后可看大图

还是没有 到底怎么回事啊 弄了快一小天了。。。

你的程序没有运行成功啊~ 成功了会输出 succeed termination 之类的语句

泷岛彗

爱的侍者 发表于 2012-12-30 14:57

                               
登录/注册后可看大图

你的程序没有运行成功啊~ 成功了会输出 succeed termination 之类的语句

可是提示没有错误啊
不是在build菜单下，编译(compile)、再点击叹号运行就是我上面图片出现的结果。。。。

G-mimzy

泷岛彗 发表于 2012-12-30 14:39

                               
登录/注册后可看大图

还是没有 到底怎么回事啊 弄了快一小天了。。。

不好意思，刚才没仔细看程序，你应该在运行窗口光标处输入 49,1956  就可以了

G-mimzy

G-mimzy 发表于 2012-12-30 16:12

                               
登录/注册后可看大图

不好意思，刚才没仔细看程序，你应该在运行窗口光标处输入 49,1956  就可以了

输完后回车，另外建议把  --,STATUS='NEW'---删掉，否则每次运行你都要先删除SHUCHU文件