| 查看: 707 | 回复: 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客户端 |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有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
20