| 查看: 673 | 回复: 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客户端 |
» 猜你喜欢
国自然上会要求
已经有9人回复
博士申请
已经有5人回复
面上本子正文33页,违规吗?会被低分嘛?
已经有9人回复
评审有感
已经有6人回复
今年审到国自然15份,谈谈感受
已经有16人回复
上海大学实验技术岗位非升即走
已经有8人回复
考博自荐
已经有6人回复
青C资助名额大幅增加!
已经有16人回复
重磅!青年科学基金项目(C类)资助增幅预计超过50%
已经有10人回复
我在等一个没有答案的答案
已经有3人回复
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