求教画连续时次图，用一个ctl修改后出错原因

tszhaoweiwen

我下载了gdas2013052000.grib2.ctl
gdas2013052006.grib2.ctl
gdas2013052012.grib2.ctl
3个时次的资料 然后想用同一个GS画图  
ctl文件：
dset ^gdas2013052000f%f2.grib2
index ^gdas2013052000.grib2.idx
undef 9.999E+20
title GDAS 00Z20may2013 downloaded May 20 07:04 UTC
*  produced by g2ctl v0.0.4i
dtype grib2
ydef 181 linear -90.000000 1
xdef 360 linear 0.000000 1.000000
tdef 2 linear 00Z20may2013 6hr  *这里要改成tdef 3 linear 00Z20may2013 6hr？？

GS脚本
'reinit'
'open C:\Users\sta288\Desktop\516\gdas2013052000.grib2.ctl'
'set mpdset cnworld'
'set rbcols auto'
'set rgb 41 255 255 255'
'set rgb 42 150 210 250'
'set rgb 43 80 165 245'
'set rgb 44 255 250 170'
'set rgb 45 255 192 60'
'set rgb 46 255 96 0'
'set rgb 55 255 226 171'
'set rgb 56 145 76 57'
'set rgb 57 176 135 56'
'set parea 1 10.0 0.5 7.90'


'set t 1'

'set lon 89.0 130.0'
'set lat 9.0 40.0'
'set lev 925'

'define t=TMPprs'
'define rh=RHprs'
'a=17.2693882'
'b=35.86'
'define es=6.1078*exp((a*(T-273.16))/(T-b))'
'define e=rh*(es/100.0)'
'define q=(0.622*e)/(lev-(0.378*e))'
*'define qdiv=mag(UGRDprs,VGRDprs)*2.5*1000*q/9.8'
*'define qdiv=1000*q/9.8'
'define qdiv=mag(UGRDprs,VGRDprs)*2.5*1000*q/9.8'
'set lon 100 120'
'set lat 15 30'
'set gxout shaded'
'set xlopts 1 4 0.2'
'set ylopts 1 4 0.2'
'set xlint 5'
'set ylint 5'
'set lev 850'
'set clevs 10 20 30 40 50'
'set rbcols 41 42 43 44 45 46'
'set grads off'
'd qdiv'
'run cbarn.gs 1 1 10.05 4.25'
*'set gxout stream'
*'set gxout barb'
'set strmden 2'
'set grads off'
'd UGRDprs*2.5;VGRDprs*2.5'
*'draw title 'qtime()'_'qlev()'hPa vaporflux div'
'printim E:/qdiv7.png  white'

把set t 2 就出错了 只能画 t1的图 请问为什么？

传说中的谁

改ctl之后gribmap了没有？

liutaoz21

帮顶一下~水平不够，没看大明白~

river

楼主啊，你是批量描述三个时次的资料是吧？为什么不加option template呢，没有这个项的话，GrADS可不会认为是批量描述的。我最近发了两个关于批量描述的帖子，楼主都没有仔细看啊

carpenter3

加一句 options yrev template 试试

loftydongshi

river 发表于 2013-5-28 19:24

                               
登录/注册后可看大图

楼主啊，你是批量描述三个时次的资料是吧？为什么不加option template呢，没有这个项的话，GrADS可不会认为 ...

没有这个选项是不能批量的，所言极是啊，river

tszhaoweiwen

river 发表于 2013-5-28 19:24

                               
登录/注册后可看大图

楼主啊，你是批量描述三个时次的资料是吧？为什么不加option template呢，没有这个项的话，GrADS可不会认为 ...

option template ctl文件就有啊 你说的贴我看了 是教用几年资料然后逐日打开  下面附ctl
dset D:\OpenGrADS\data\gdas2013052000f%f2.grib2
index D:\OpenGrADS\data\gdas2013052000.grib2.idx
undef 9.999E+20
title GDAS 00Z20may2013 downloaded May 20 07:04 UTC
*  produced by g2ctl v0.0.4i
dtype grib2
ydef 181 linear -90.000000 1
xdef 360 linear 0.000000 1.000000
tdef 2 linear 00Z20may2013 6hr
* PROFILE hPa
zdef 26 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000 7000 5000 3000 2000 1000
options pascals template
vars 139
no4LFTXsfc   0,1,0   0,7,193 ** surface Best (4 layer) Lifted Index [K]

tszhaoweiwen

loftydongshi 发表于 2013-5-29 10:15

                               
登录/注册后可看大图

没有这个选项是不能批量的，所言极是啊，river

您好，tdef 2 linear 00Z20may2013 6hr options pascals template ctl里是这样的 他们是配套使用的啊 还是不能批处理  

tszhaoweiwen

传说中的谁 发表于 2013-5-28 19:08

                               
登录/注册后可看大图

改ctl之后gribmap了没有？

版主，论坛的贴都教改tdef 2 linear 00Z20may2013 6hr options pascals template  您提到的gribmap如何改 不是很懂

传说中的谁

tszhaoweiwen 发表于 2013-5-29 22:03

                               
登录/注册后可看大图

版主，论坛的贴都教改tdef 2 linear 00Z20may2013 6hr options pascals template  您提到的gribmap如何改 ...

看起来楼主对于grads不是很了解。。。。ctl是用来描述数据文件的，但并不直接指向数据文件，而是需要通过索引文件来指向数据文件，index ^gdas2013052000.grib2.idx，这个就是所谓的索引文件，你只修改ctl，虽然看起来它可以描述多个文件了，但实际上索引文件还是旧的，所以实际效果跟没修改的时候是一样的。至于怎么使用gribmap，还请楼主自己捉摸一下，可以参考http://bbs.06climate.com/forum.php?mod=viewthread&tid=2465