处理好的12个月的高度平均场，为啥GRADS只画出一张图！！！！！

努力的叮宝

用FORTRAN编了一个求了35年月平均高度场，写成grd的格式，然后想用grads画图，所以又写了一个ctl，但是问题来了。。。画出的12个月的图都是一样的，我自己琢磨很久，不知道什么原因
~~~本人菜鸟，已经弄了很久。。。特此求助各位大神！！！

努力的叮宝

补上FORTRAN程序
open(12,file='d:\data\climate.grd',form='binary')
  do n=1,12
        do i=1,73
      do j=1,37
                 sum=0.0
            do t=n,420,12
                   sum=sum+var(i,j,t)
        end do
        ave(i,j,n)=sum/35
          end do
     end do
   end do
write(12) (((ave(i,j,t),i=1,ii),j=1,jj),t=1,12)
close(12)
end

river

看着没什么大问题啊，你的ctl里面时间项里的dims是什么意思

努力的叮宝

river 发表于 2014-10-18 16:42
看着没什么大问题啊，你的ctl里面时间项里的dims是什么意思

写错了，后来删了了。。。╮(╯▽╰)╭

river

努力的叮宝 发表于 2014-10-18 16:48
写错了，后来删了了。。。╮(╯▽╰)╭

那你把完全版的fortran贴上来看看

努力的叮宝

river 发表于 2014-10-18 17:00
那你把完全版的fortran贴上来看看

program ex_01
implicit none
integer i,j,t,n,ii,jj,nt,irec
parameter(ii=73,jj=37,nt=420)
real var(ii,jj,nt),sum,ave(ii,jj,12)

!读取数据
open(9,file='d:\data\z79_13.grd',form='binary')
    do t=1,420
        irec=1
    read(9,rec=irec) ((var(i,j,t),i=1,ii),j=1,jj)
    irec=irec+1
    end do
close(9)
write(*,*) var(1,1,3)

!计算1979-2013年月平均高度场
open(12,file='d:\data\climate.grd',form='binary')
  do n=1,12
        do i=1,73
      do j=1,37
                 sum=0.0
            do t=n,420,12
                   sum=sum+var(i,j,t)
        end do
        ave(i,j,n)=sum/35
          end do
     end do
   end do
write(12) (((ave(i,j,t),i=1,ii),j=1,jj),t=1,12)
close(12)
end

mofangbao

怎么画的？

努力的叮宝

river 发表于 2014-10-18 17:00
那你把完全版的fortran贴上来看看

解决了，数据读的有问题

努力的叮宝

mofangbao 发表于 2014-10-18 17:19
怎么画的？

就直接画啊，数据处理好了呀

萝卜都

学习{:5_213:}