implicit none
integer,parameter::n=365
real jg(365,1),js(365,1),qjs,a
integer::i,j,k,g
js=0.0
jg=0.0
a=0.75
k=1
open(11,file='801.txt')
read(11,*)(jg(i,1),i=1,n)
close(11)
do g=1,7
js(g,1)=jg(g,1)
end do
  do i=8,365
    if(jg(i,1)==0.0)then
      js(i,1)=jg(i,1)
      else if(abs(jg(i,1)).gt.0.0)then
        qjs=0
        do j=1,7
            qjs=jg(i-j,1)*a**(i-j)+qjs
            js(i,1)=jg(i,1)+qjs
        end do 
     end if
   end do
open(15,file='yxjs1980.txt')
write(15,*)(js(k,1),k=1,n)
close(15)
end
我用这个程序算个有效值，但是不晓得为什么计算出来的结果前面部分还有改变，后面的就和源数据一模一样，我是用前7天的值乘以个衰减的系数加上当天的数值，求达人帮忙看看，数据和程序都在附件里