grads将tbb与风场叠加

IAMQWW1

刚接触grads，想请教一下，如何用grads将tbb与风场进行叠加；
ctl如下

dset C:\Users\IAMQWW\Desktop\miss\tbb\tbb\chengxu\FY2D_TBB_IR1_OTG_20120620_0015.AWX
title 2012 twpice fy2d_tbb
options yrev
* yrev表示y轴反向
undef 0.0
xdef 1201 linear  45 0.1
ydef 1201 linear -60 0.1
zdef 1 linear 1000 1
tdef 1 linear 00z20JUN2012 1hr
vars  1
tbb   1    -1,40,1    Geopotential Heights (m)
endvars

gs：
'reinit'
file='C:\Users\IAMQWW\Desktop\miss\tbb\tbb\chengxu\tbb.gmf'
'enable print 'file
'sdfopen C:\Users\IAMQWW\Desktop\miss\tbb\tbb\chengxu\uwnd.2012.nc'
'sdfopen C:\Users\IAMQWW\Desktop\miss\tbb\tbb\chengxu\vwnd.2012.nc'

'set mpdset hires'
'set grads off'
'set grid off'
'set lev 500'
'set t 29'
'set lat 15 65'
'set lon 70 140'
'set gxout vector'
'd uwnd;vwnd.2'
'close 1'
'close 2'

'open C:\Users\IAMQWW\Desktop\miss\tbb\tbb\chengxu\a.ctl'
'set lat 15 65'
'set lon 70 140'
'd tbb+100'
'cbarn'
'print'
'clear'
'disable print'
'printim C:\Users\IAMQWW\Desktop\miss\tbb\tbb\chengxu\tbb.png white'
'reinit'
;

balfulosa

大概原因可能是tbb数据和nc的风场数据时间、层次对应不上。换个顺序，先d tbb，然后在set lev 、set t风场试一试。

IAMQWW1

balfulosa 发表于 2016-1-20 14:58
大概原因可能是tbb数据和nc的风场数据时间、层次对应不上。换个顺序，先d tbb，然后在set lev 、set t风场 ...

谢谢你的回答，我将nc中set t设置成了跟tbb是一个时次的了，但是出现下面的错误
C:\Users\IAMQWW\Desktop

听海落雪

不可以画。因为nc资料和tbb资料的分辨率不一样。

雪颜

听海落雪 发表于 2016-9-22 13:26
不可以画。因为nc资料和tbb资料的分辨率不一样。

可以画的.....

雪颜

dfile 1
dfile 2

听海落雪

雪颜 发表于 2016-9-22 14:11
dfile 1
dfile 2
并且nc的资料下载时请选择于TBB的范围一样

好的，我试试，谢谢。

听海落雪

雪颜 发表于 2016-9-22 14:11
dfile 1
dfile 2
并且nc的资料下载时请选择于TBB的范围一样

'reinit'

'open F:\Data\GrADS\Typhoon_test\TBB\FY2G_TBB_IR1.ctl'
*'set dfile 1'
'sdfopen F:\Data\GrADS\2016\20160123\hgt.2016.nc'
*'set dfile 2'

'set grads off'
'set grid off'

*'set vpage 1 10  1 8'
*'set grads off'

****** 蓝色系 ******
'set rgb 22  200 210 245'
'set rgb 23  110 150 240'
'set rgb 24   80 100 240'
'set rgb 25   20  60 230'
'set rgb 26    0  20 205'
'set rgb 27    0   0 160'

'set lon 80 110'
'set lat 10 30'

'set xlopts 1 5 0.28'
'set ylopts 1 5 0.28'
'set xlab on'
'set ylab on'
'set xlabs 80`3.|85`3.|90`3.|95`3.|100`3.|105`3.|110`3.`1E'
'set ylabs 10`3.|15`3.|20`3.|25`3.|30`3.`1N'

'set t 2'
'set gxout shaded'
'set csmooth on'
'set clevs -82 -72 -62 -52 -42 -32'
'set ccols 27 26 25 24 23 22 0'
'set cint 20'
'd tbb.1-173.15'

*'set t 558'
'set lev 500'
'set gxout contour'
'set csmooth on'
'set cint 2'
'set cstyle 1'
'set ccolor 1'
'set clopts 1  0.5  0.15'
'set clab forced'
'set clab masked'
'd hgt.2/10'

'set mpdset cnworld'
'set map 15 1 2'
'draw map'
'set mpdset yunn'
'set map 8 1 5'
'draw map'

'query w2xy 82.0 15.4'
x1=subwrd(result,3)
y1=subwrd(result,6)
'set string 1 c 5 0'
'set strsiz 0.16'
'draw wxsym 41 'x1' 'y1' 0.5'

'set string 1 c 6 0'
'set strsiz 0.58'
'draw string 1.8  7.3   (a)'

'set string 1 c 5 0'
'set strsiz 0.35'
'draw string 8.0 1.05   14BJZ19MAY'


*** run plotcbar rainbow.gs的路径 色标的 刻度值间隔 大小 垂直/水平 水平位置 垂直位置***
'run plotcbar F:\Data\GrADS\Typhoon_test\TBB\rainbow.gs 1 1.2 1 10.3 4.2'

'printim F:\Data\GrADS\Typhoon_test\test_03.png white'


当我把 TBB资料的时间间隔设置成与nc资料一样（即6小时），先画TBB，再hgt，就能出图了，而且二者的时间一致。我并没有用到“set dfile..."。

听海落雪

雪颜 发表于 2016-9-22 14:11
dfile 1
dfile 2
并且nc的资料下载时请选择于TBB的范围一样

这是我画的图。

雪颜

听海落雪 发表于 2016-9-22 22:49
这是我画的图。

dfile可以隔离画图，直接d就可以了，比较不容易出错，变量名后面不需要加.1 .2，特别是循环出图时好用

例如

'set dfile 1'
'set gxout shaded'
'd tbb'

'set dfile 2'
'set gxout contour'
'd hgt'