| 查看: 700 | 回复: 2 | ||
[求助]
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),X1(1000),X2(1000),X3(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='D:\fortran\xxbb\MK1.txt') READ(2,*) X1(I),X2(I),X3(I) DO 60 I=1,N 60 Y(I)=(X1(I)+X2(I)+X3(I))/3 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='D:\fortran\xxbb\mk2.txt',STATUS='NEW') A=1.96 B=-1.96 DO 40 I=1,N WRITE(3,50)NYEAR+I-1,Y(I),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 !GE是>=，GT是>，LE是<=,LT是< IF(Y(I).LE.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  发自小木虫Android客户端 |
回复此楼» 猜你喜欢
面上的水很深
已经有17人回复
-
国社科通讯评审
已经有5人回复
-
5个%2F，非固定段22个字母
已经有9人回复
-
生命口C13面上会评时间
已经有11人回复
-
内心匮乏
已经有16人回复
-
生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？
已经有9人回复
-
投英文期刊如何查重
已经有6人回复
-
不要再数国自然申请书的 filecode 的分隔符个数了
已经有21人回复
-
伯胺上甲基求助
已经有7人回复
-
跨出版社商投稿
已经有5人回复
seanbsd
木虫 (小有名气)
- 应助: 8 (幼儿园)
- 金币: 3652.1
- 红花: 1
- 帖子: 207
- 在线: 45.6小时
- 虫号: 3744514
- 注册: 2015-03-17
- 专业: 电气科学与工程
2楼2017-04-07 07:12:03
秋犯
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 578.8
- 沙发: 7
- 帖子: 331
- 在线: 37.5小时
- 虫号: 3427738
- 注册: 2014-09-19
- 性别: GG
- 专业: 环境污染化学
3楼2017-04-07 21:08:16
10