麻烦大家帮忙看个程序

jeffggg

parameter (n=3650,m=73,nl=160)
real A(n),S(n)
integer i,j,k,o
character*5 needname(nl)
open(10,file='H:\filename.txt')
do i=1,nl
print*,i
read(10,'(a5)') needname(i)
enddo
close(10)
write(*,*) 'read OK!'
do i=1,nl
open(30,file='H:\AAA\'//needname(i)//'.grd',position='Append',form='binary')
open(20,file='H:\'//needname(i)//'.grd',form='binary')
do j=1,n
  read(20) A(j)
   do k=1,73
     do o=k,3577,73
     if(j<=3650)then
      A(k)=A(k)+A(o+73)
  endif
     enddo
      S(k)=A(k)/73
   enddo
enddo
      write(30) S(k)
      close(30)
   close(20)
enddo   
有160个站的GRD文件 ，每个文件里是这个站点50年的 每侯的数据 即50*73=3650个 现在我通过批处理的方式 把每个GRD里的数据求个气候平均，即 50个第一侯求平均 50个第二侯求平均依次类推，最后再输出为GRD文件 但是我现在这个程序能运行 但是 运行后生成的每个GRD只有4字节的数据 肯定不对 而且额画图也不对，正常应该改输出73*4=292字节的数据的
麻烦大家啦   

mofangbao

你是想输出s吧，写s(k)当然只输出了一个元素，提问前先仔细检查啊，但愿你是粗心造成，要不然实在是基础有点欠缺

jeffggg

mofangbao 发表于 2013-6-9 17:46

                               
登录/注册后可看大图

你是想输出s吧，写s(k)当然只输出了一个元素，提问前先仔细检查啊，但愿你是粗心造成，要不然实在是基础有点 ...

清风版主 我编程基础确实很差 但是我不理解 我前面定义的S是数组啊 为什么最后输出值输出S呢？

loftydongshi

S(k)这不是表示输出的是一个数么，你把write那行提到  S(k)=A(k)/73下面就对了！太马虎了