循环or判断的问题

yuqianhao0315

volcanic11000.f90 (833 Bytes, 下载次数: 1)

2014-10-5 16:35 上传

点击文件名下载附件

volcanic11000.txt (15.16 KB, 下载次数: 1)

2014-10-5 16:35 上传

点击文件名下载附件

各位假期工作在前线的劳模朋友们好！小妹遇到了困扰了两天的程序问题，求解答~

数据说明：
数据一共有940行。第一列是年份(yearBP),每个值要用1950减一下，才表示公元的年份；第二列和第三列是两种变量值。想求10年平均的两种变量值，比如1980年-1990年的平均值，算作1985年这一时刻的变量值。

程序说明：
要求10年平均，根据1950-yr(i)来判断这行数据是属于哪个时间段的。
比如第一行，是1985.26年，判断它是1980到1990年这个时间段的，再往下走，第二三四行都是这个时间段的，后两列表示两种变量，求他们前4行的平均，分别记做1985年的totalsulsum和sulsum。
看看时间范围，是从1985.26到3.9年，最后能求出变量的时间点是1985、1975、1965、……、15、5年，也就是说输出结果共199行。
1985表示第一个时间点，year=1990时，200-year/10的值就是1，表示了输出结果的第一行；
1975表示第二个时间点，year=1980时，200-year/10的值是2，表示输出结果的第二行；
……
5表示第199个时间点，year=10时，200-year/10的值是199，表示输出结果的最后一行。

问题：
输出结果没有预想的做了10年平均，而是都是同一个值，if没有起到约束作用。下面贴出数据和fortran程序，求大神帮找出问题所在~

lqouc

没几行的程序就请楼主直接把代码贴出来吧

yuqianhao0315

lqouc 发表于 2014-10-5 16:54
没几行的程序就请楼主直接把代码贴出来吧

好滴！
program l
implicit none

real,dimension(940)::yr,totalsul,sul
real,dimension(199)::n,totalsulsum,sulsum
integer::i,m,year

open(1,file='volcanic11000.txt')
open(2,file='volcanic11000new.txt')

read(1,*)(yr(i),totalsul(i),sul(i),i=1,940)
write(*,*)(yr(i),totalsul(i),sul(i),i=1,940)

n(1:199)=0.0
totalsulsum(1:199)=0.0
sulsum(1:199)=0.0

do year=1990,10,-10
do i=1,940
  if((year-10)<=(1950-yr(i))<year)then
    n(200-year/10)=n(200-year/10)+1
    totalsulsum(200-year/10)=totalsulsum(200-year/10)+totalsul(i)
    sulsum(200-year/10)=sulsum(200-year/10)+sul(i)
   endif
  enddo
totalsulsum(200-year/10)=totalsulsum(200-year/10)/n(200-year/10)
sulsum(200-year/10)=sulsum(200-year/10)/n(200-year/10)
enddo

m=5
do i=199,1,-1
write(2,*)m,totalsulsum(i),sulsum(i)
m=m+10
enddo

end program l

排山倒海

if((year-10)<=(1950-yr(i))<year)then
这个明显是c语言的规则啊，不是fortran的规则吧
要改成 （year-10)<=(1950-yr(i)).and.(1950-yr(i))<year

yuqianhao0315

排山倒海 发表于 2014-10-5 18:48
if((year-10)

嗯嗯，改过了，果然是这样！好厉害啊你，赞一个赞一个~

topmad

不会用C语言的表示一下子还看懵了

yuqianhao0315

topmad 发表于 2014-10-6 14:40
不会用C语言的表示一下子还看懵了

我也不会，只是按照逻辑就写成这样了……哈哈，楼上是高手的耶！