关于线性倾向估计的FORTRAN程序问题

呆萌萱

我用魏凤英老师的fortran程序计算线性倾向估计
spring600.txt
  272.9870   
   272.9634   
   272.9492   
   272.9411   
   272.9330   
   272.9222   
   272.9073   
   272.8877   
   272.8614   
   272.8297   
   272.7959   
   272.7628   
   272.7351   
   272.7155   
   272.7054   
   272.7007   
   272.6993   
   272.6986   
   272.7027   
   272.7101   
   272.7195   
   272.7297   
   272.7371   
   272.7419   
   272.7473   
   272.7594   
   272.7776   
   272.8019   
   272.8229   
   272.8384   
   272.8472   
   272.8560   
   272.8708   
   272.8951   
   272.9269   
   272.9607   

year.txt     1980
   1981
   1982
   1983
   1984
   1985
   1986
   1987
   1988
   1989
   1990
   1991
   1992
   1993
   1994
   1995
   1996
   1997
   1998
   1999
   2000
   2001
   2002
   2003
   2004
   2005
   2006
   2007
   2008
   2009
   2010
   2011
   2012
   2013
   2014
   2015

程序如下
C        THIS IS A PROGRAM OF LINEAR REGRESSION BETWEEN
C        X AND T(TIME OR YEAR),i.e. X=A+BT
C        PROGRAM LR
        PARAMETER (N=34)!常数声明
        DIMENSION X(N),IT(N),RX(N)!定义数组
        REAL A,B,R
C        ***************************************************
C        * N:     SAMPLE SIZE                              *
C        * X:     ARRAY OF X                               *
C        * IT:    ARRAY OF YEAR                            *
C        ***************************************************
        OPEN(2,FILE='spring600.TXT')
        READ(2,*)(X(I),I=1,N)
        OPEN(3,FILE='year.TXT')
        READ(3,*)(IT(I),I=1,N)
        CALL SLR(N,X,IT,RX,A,B,R)
        OPEN(6,FILE='LR600.TXT',STATUS='NEW')
        WRITE(6,10)
  10        FORMAT(20X,'LINEAR TREND ANALYSIS'/)
          WRITE(6,20)A,B,R
  20        FORMAT(5X,'CONSTANT OF REGRESSION A=',F10.4/
     &        5X,'COEFFICIENT OF REGRESSION B=',F10.4/
     &        5X,'CORRELATION COEFFICIENT R=',F10.4/5X,'--------------------')
             DO 30 I=1,N
        WRITE(6,40)IT(I),X(I),RX(I)
  40        FORMAT(I6,F7.2,F7.2)
  30        CONTINUE
          STOP
        END
C********************************************************************
        SUBROUTINE SLR(N,X,IT,RX,A,B,R)
        DIMENSION X(N),IT(N),RX(N)
        REAL A,B,R
        CX=0.0
        CT=0.0
        CXC=0.0
        CTC=0.0
        CXT=0.0
        DO 10 I=1,N
        CX=CX+X(I)
        CT=CT+IT(I)
        CXC=CXC+X(I)*X(I)
        CTC=CTC+IT(I)*IT(I)
  10        CXT=CXT+X(I)*IT(I)
          DX=CXC-CX*CX/FLOAT(N)!FLOAT强制转换成浮点
        DT=CTC-CT*CT/FLOAT(N)
        DXT=CXT-CX*CT/FLOAT(N)
        B=DXT/DT
        A=CX/FLOAT(N)-B*CT/FLOAT(N)
        R=DXT/SQRT(DX*DT)
        DO 20 I=1,N
  20        RX(I)=A+B*IT(I)
        RETURN
        END
C        ---------------------END-------------------------------------

但是运行出来的错误如图所示，这是什么问题啊？

Masterpiece

【系统】标号3的文件，本宝宝读着读着就到头了，然而循环这个老家伙还让本宝宝继续读

呆萌萱

Masterpiece 发表于 2016-4-27 21:30
【系统】标号3的文件，本宝宝读着读着就到头了，然而循环这个老家伙还让本宝宝继续读

那这个问题该怎么改啊  我试过它读到2014年就不读了

Masterpiece

呆萌萱 发表于 2016-4-27 21:45
那这个问题该怎么改啊  我试过它读到2014年就不读了

1980年到2015年一共有36年，当N大于36的时候就可能报错了
要么是你的year.txt里边少了2015

呆萌萱

Masterpiece 发表于 2016-4-27 21:48
1980年到2015年一共有36年，当N大于36的时候就可能报错了
要么是你的year.txt里边少了2015

我的year.txt里就是1980-20115，没少数据，但是就是读year.txt有错啊，怎么办啊  好苦恼