grads画图出现彩虹条

黄鸿惠

我在利用u、v风求相对涡度并转dat，
想知道为什么会出现转了dat后出来的图会变成这样，不转dat直接出图就不会
感觉是我转dat的时候出问题了，脚本写的是不是不对，劳烦大神们帮我看看。

黄鸿惠

第一次提问没经验，我把三个gs放这里了
这是计算相对涡度并转dat
'reinit'
'sdfopen F:\taifeng\chazhi\xwodu\uwnd.mon.mean.nc'
'sdfopen F:\taifeng\chazhi\xwodu\vwnd.mon.mean.nc'
'set lon 100 270'
'set lat -10 30'
'set lev 500'
'set parea 0.5 10.8 0 8.5'

'rvo=hcurl(uwnd.1,vwnd.2)'
'd rvo'

'set gxout fwrite'
'set fwrite  F:\taifeng\chazhi\xwodu\rvo1.dat'
tt=1
while (tt<=858)
'd rvo'
tt=tt+1
endwhile
'disable fwrite '

这是对应的ctl
dset  F:\taifeng\chazhi\xwodu\rvo1.dat
title monthly mean 500hpa relatively vorticity
undef 1e+20
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 levels 500
tdef 858 linear 00z01JAN1948 1mo
vars 1
rvo 1 99 the monthly mean 500hpa relatively vorticity
endvars

这是调用生成的ctl的gs
'reinit'
'open F:\taifeng\chazhi\xwodu\rvo.ctl'
'set lon 100 270'
'set lat -20 30'

'd rvo'

'draw title CP-EP rvo Difference Map and T-Test'
'printim F:\taifeng\chazhi\xwodu\Cp_Ep_rvo_3.png x1000 y800 white'

四叶草

明显不是一个数据，数据有问题，你先提取一个时次的看看，还有看看原文件的描述文件，或者搜索家园里的相关帖子学习下。nc数据本人也没这样提取过，给不出太多建议。