【求助】多时次的垂直剖面图

quentin011099

各位大侠，想画这样的多时次的垂直剖面图。我是先在fortran里面定义了一个四维数组，把所有数据都存进去了。然后用下面的GS文件作图，可是做不出来，求各位大侠的帮助，先谢谢啦~~~
FORTRAN读取和存储数据：  
data time/'12071920.000','12072008.000','12072020.000','12072108.000','12072120.000','12072208.000','12072220.000','12072220.000'/
  data level/'1000','925','850','700','500','400','300','250','200','150','100'/
!*************************************READ SE*********************************************!
do iz=1,nz
do it=1,nt
  open(1,file='f:\program\micaps\pou\thetse\'//trim(level(iz))//'\'//time(it),err=101)  
   do i=1,4          !!!前四行进行空读
    read(1,*)           
   enddo
   do iy=ny,1,-1
    read(1,*) (var(ix,iy,iz,it),ix=1,nx)
   enddo
  close(1)
enddo
enddo

!********************************WRITE****************************************************!
open(1,file='f:\program\micaps\pou\dat\721.dat',form="binary",err=102)
do it=1,nt
  write(1) (((var(ix,iy,iz,it),ix=1,nx),iy=1,ny),iz=1,nz)
enddo
close(1)

CTL文件：
dset f:\program\micaps\pou\dat\721.dat
undef 9999
title 12年7月19日-22日假相当位温剖面图
xdef 53 linear 30 2.5
ydef 29 linear 10 2.5
zdef 11 levels 1000 925 850 700 500 400 300 250 200 150 100
tdef 8 linear 08z19jul2012 12hr
vars 1
se 1 99 diamond4tograds
endvars

GS文件
'open f:\program\micaps\pou\ctl\721.ctl'
'set grads off'
'set grid off'
'set lon 105'
'set lat 18 55'
'set lev 1000 100'
'set t 1 8'
'set gxout contour'
'd se'
'printim f:\program\micaps\pou\map\721'it'.png x1000 y800 white'
'c'
'reinit'

mofangbao

这种方式用来分享资料和代码还不错，但是用来提问就不合适了
你需要自己先提炼 定位你的出错点，然后介绍你的东西
要不然别人得花多少时间帮你从头到尾看？

quentin011099

mofangbao 发表于 2012-8-10 16:42

                               
登录/注册后可看大图

这种方式用来分享资料和代码还不错，但是用来提问就不合适了
你需要自己先提炼 定位你的出错点，然后介绍你 ...

啊，那我回去再好好研究一下，看看问题出在哪里，先谢谢了···

kanbaoyun

学习了，学习了

啊珏珏珏

学习学习 正在找这类型的资料~~！