我把8个时次的grd空间风场用一个ctl描述，为什么画出的图各个时次都是一样的呢？谢谢

大果子

v345.ctl (216 Bytes, 下载次数: 4)

2013-7-27 18:05 上传

点击文件名下载附件

lqouc

请直接贴出来你的ctl，还有如果数据是自己写的，把写数据的程序也直接贴出来。

大果子

生成数据的程序
program main
implicit none
real,parameter::nn=8,nx=144,ny=45
real::a(nx,ny),b(nx,ny),c(nx,ny),sum1(nx,ny),ave(nx,ny)
integer::n,i,j
    character*2::nam1(nn)=(/'1','2','3','4','5','6','7','8'/)
   
do n=1,nn
      open(1,file='d:\MJO\v200'//trim(nam1(n))//'-3.grd',form='binary')
      open(2,file='d:\MJO\v200'//trim(nam1(n))//'-4.grd',form='binary')
      open(3,file='d:\MJO\v200'//trim(nam1(n))//'-5.grd',form='binary')
      read(1) ((a(i,j),i=1,nx),j=1,ny)
   close(1)
      read(2) ((b(i,j),i=1,nx),j=1,ny)
   close(2)
   read(3) ((c(i,j),i=1,nx),j=1,ny)
   close(3)
     
   open(10,file='d:\MJO\v200'//trim(nam1(n))//'-345.grd',form='binary')
    open(11,file='d:\MJO\v200'//trim(nam1(n))//'-345.txt')

do j=1,ny
do i=1,nx
sum1(i,j)=a(i,j)+b(i,j)+c(i,j)
ave(i,j)=sum1(i,j)/3
write(10) ave(i,j)
    write(11,*) ave(i,j)
enddo
enddo
    close(10)
enddo
  end

大果子

ctl
dset d:\MJO\v200%d1-345.grd
undef -9.99+E33
title season ave
options template
xdef 144 linear 0.00 2.5
ydef 45 linear -40.00 2.5
zdef 1 levels 200
tdef 8 linear 01jan1981 1dy
vars 1
v 1 99 v-wind
endvars
;

river

大果子 发表于 2013-7-27 19:44

                               
登录/注册后可看大图

生成数据的程序
program main
implicit none

你不是分别写了二进制和十进制文件出来吗，先打开生成十进制的资料看看数据有没有问题，如果有问题那就是fortran程序的问题了。如果没有那就可能是ctl和画图的gs有问题，不过看你的ctl没什么大问题，就是ctl结尾不需要 ; 来结尾，去掉看看。如果也不是这个问题，那就把画图的gs也贴上来吧

大果子

river 发表于 2013-7-27 21:59

                               
登录/注册后可看大图

你不是分别写了二进制和十进制文件出来吗，先打开生成十进制的资料看看数据有没有问题，如果有问题那就是 ...

不好意思，是我的原始资料出了问题，现在已经解决了，谢谢！