24小时热门版块排行榜

返回列表

Gimo

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 3
帖子: 1
在线: 20分钟
虫号: 6252786
注册: 2017-04-06
专业: 经济史

[求助] 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客户端

回复此楼