调试该源代码求累计距平，有一个错误未能找出来

水哥哥8663

1. C  THIS IS A PROGRAM FOR CALCULATING CUMULATIVE
   
2. C  ANOMALY OF A TIME SERIES
   
3.   PROGRAM CA
   
4.   DIMENSION X(1000),X1(1000)
   
5.   WRITE(*,10)
   
6.   10  FORMAT(5X,'N=?,NYEAR=?')
   
7.   READ(*,*)N,NYEAR
   
8. C  ********************************************
   
9. C  * N:     SAMPLE SIZE                       *
   
10. C  * NYEAR: FIRST YEAR OF THE TIME SERIES     *
    
11. C    * X(N):  ORIGINAL TIME SERIES              *
    
12. C  * X1(N): CUMULATED ANOMALY SERIES          *
    
13. C  ********************************************
    
14.   OPEN(2,FILE=' ')
    
15.   READ(2,*)(X(I),I=1,N)
    
16.   C=0.0
    
17.   DO 20 I=1,N
    
18.   20  C=C+X(I)
    
19.   XM=C/N
    
20.   DO 30 I=1,N
    
21.   C=0.0
    
22.   DO 40 J=1,I
    
23.   40  C=C+(X(J)-XM)
    
24.   X1(I)=C
    
25.   30  CONTINUE
    
26.   OPEN(3,FILE=' ',STATUS='NEW')
    
27.   WRITE(3,50)
    
28.   50  FORMAT(30X,'CUMULATIVE ANOMALY'/)
    
29.   DO 60 I=1,N
    
30.   WRITE(3,70)NYEAR-1+I,X1(I)
    
31.   70  FORMAT(1X,I4,F8.2)
    
32.   60  CONTINUE
    
33.   STOP
    
34.   END
    

C  THIS IS A PROGRAM FOR CALCULATING CUMULATIVE<br /> C  ANOMALY OF A TIME SERIES<br />   PROGRAM CA<br />   DIMENSION X(1000),X1(1000)<br />   WRITE(*,10)<br />   10  FORMAT(5X,'N=?,NYEAR=?')<br />   READ(*,*)N,NYEAR<br /> C  ********************************************<br /> C  * N:     SAMPLE SIZE                       *<br /> C  * NYEAR: FIRST YEAR OF THE TIME SERIES     *<br /> C    * X(N):  ORIGINAL TIME SERIES              *<br /> C  * X1(N): CUMULATED ANOMALY SERIES          *<br /> C  ********************************************<br />   OPEN(2,FILE=' ')<br />   READ(2,*)(X(I),I=1,N)<br />   C=0.0<br />   DO 20 I=1,N<br />   20  C=C+X(I)<br />   XM=C/N<br />   DO 30 I=1,N<br />   C=0.0<br />   DO 40 J=1,I<br />   40  C=C+(X(J)-XM)<br />   X1(I)=C<br />   30  CONTINUE<br />   OPEN(3,FILE=' ',STATUS='NEW')<br />   WRITE(3,50)<br />   50  FORMAT(30X,'CUMULATIVE ANOMALY'/)<br />   DO 60 I=1,N<br />   WRITE(3,70)NYEAR-1+I,X1(I)<br />   70  FORMAT(1X,I4,F8.2)<br />   60  CONTINUE<br />   STOP<br />   END<br />

mofangbao

不会fortran就学吧  这种提问方式没有会给你答案的！因为根本就没法帮你看，你还不如把标题改为分享累计距平的源代码呢，程序都是要根据自己的需求修改的

水哥哥8663

mofangbao 发表于 2012-4-26 11:43

                               
登录/注册后可看大图

不会fortran就学吧  这种提问方式没有会给你答案的！因为根本就没法帮你看，你还不如把标题改为分享累计距平 ...

哦  晓得了