【求助】多层次多时次的剖面图

quentin011099

各位大侠，前几天我在论坛上发了一个帖子，我回去后也仔细的研究了一下，还是没搞明白如何画出多层次多时次的垂直剖面图（如附图），主要的问题是坐标图画得出来，但是不出图，数据不循环，肯定是GS文件写错了。我是先用fortran读取了11层次7时次的数据，然后GS文件如下，请各位大侠帮忙，在此谢谢了！！！

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

做个霸气的木头

这图不是对的吗？有什么问题？

abd

那楼主画过其他二维图吗？比如x,y平面的？ 其实时间垂直剖面的原理也完全一样，吧x当做时间，把y当做高度不就行了吗！?

quentin011099

做个霸气的木头 发表于 2012-8-19 12:39

                               
登录/注册后可看大图

这图不是对的吗？有什么问题？

这个不是我画的，2我就想画这样的，但是不会啊

quentin011099

abd 发表于 2012-8-19 14:25

                               
登录/注册后可看大图

那楼主画过其他二维图吗？比如x,y平面的？ 其实时间垂直剖面的原理也完全一样，吧x当做时间，把y当做高度不 ...

把Y当高度那时次和层次呢？请指点···

做个霸气的木头

quentin011099 发表于 2012-8-20 09:03

                               
登录/注册后可看大图

这个不是我画的，2我就想画这样的，但是不会啊

这种剖面图，其实就是让变量只在时间和层次上面变化，其他的维数都固定。你的gs应该没有问题，你把你提取数据的gs、提取出来数据的ctl，以及错的图都贴上来，大家才好帮你分析原因呐~

quentin011099

做个霸气的木头 发表于 2012-8-20 12:29

                               
登录/注册后可看大图

这种剖面图，其实就是让变量只在时间和层次上面变化，其他的维数都固定。你的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=1，ny
     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=ny,1,-1),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


画图的时候就是什么也没有，只有横纵坐标，，请您指点啊···

做个霸气的木头

quentin011099 发表于 2012-8-20 16:30

                               
登录/注册后可看大图

FORTRAN读取和存储数据：  
data time/'12071920.000','12072008.000','12072020.000','12072108.000', ...

提取的数据是11层，那么ctl倒数第二行se后面应该是11才对，这是基础性错误，楼主要细心啊！

abd

quentin011099 发表于 2012-8-20 09:04

                               
登录/注册后可看大图

把Y当高度那时次和层次呢？请指点···

Y就是层次，X就是时次啊！层次就是高度，时次不就是时间吗？？

quentin011099

做个霸气的木头 发表于 2012-8-20 18:38

                               
登录/注册后可看大图

提取的数据是11层，那么ctl倒数第二行se后面应该是11才对，这是基础性错误，楼主要细心啊！

啊，好抱歉，是我太粗心太过于着急没好好检查。谢谢啦，出图了，实在是太感谢了~~~