fnl资料用grads作图，怎么会像这样啊

林宝宝

CTL里面描述的是全球的数据啊，我的gs范围给的是90-140和20-45，怎么出来的图就少了很多啊。我试过将经纬度再缩小，但是图形也跟着往里面缩小，我去，这是什么问题造成的啊。我的gs是这样的：
'reinit'
'open F:\fnl\fnl_20130525_06_00.ctl'
'set mpdset cnworld'
'set grid off'
'set display color white'
'c'
'set grads off'

'set csmooth on'
'define t=tmpprs'
'define rh=rhprs'
'define prs=lev'
'define es=(6.112*exp(17.67*(t-273.15)/(t-29.65)))'
'define q=rh*(0.62197*es/(prs-es))/100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(t)-log(e)-4.805)'
'define theta=t*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'

'set lat 20 45'
'set lon 90 140'
'set lev 1000 '
'set clevs 330 335 340 345 350 355 360 365 370 375 380'
'set ccols 14 9 4 11 5 10 7 12 8 2 6'

'd eqt'
'cbarn'
'draw title eqt of 2013052506'
'printim F:\fnl\23-18-117E-eqt.jpg'
'print'

林宝宝

谁来救救我啊

mofangbao

gs哪去了？

林宝宝

清风救我啊，用的单位的电脑：XP系统，grads 1.9版本。。。

mofangbao

这个问题有人问过  我也没做过这个  但是答案貌似是把计算区域扩大一些就好了  具体的忘了  你搜索一下吧

三江魂

'set lat 20 45'
'set lon 90 140'
'set lev 1000 '
以上放到
'set csmooth on'
之前试试。

njzqxt

'set display color white'
'c'
'set lon 1 359'
'set lat 1 89'
'set grads off'
试试。也可能那时的资料不全吧。

xiayu1990426

有可能是海拔高度的问题，你的数据是1000hpa，对应没数值的那些地方是高原，肯定没值啊

白兔糖0217

新手新手。。。。望不吝赐教。这些都是在类似DOS的窗口下直接敲入就可以出来图形的吗？还是还要借助于其他什么工具？？？

river

估计是资料的问题，gs没有问题的。这个计算不涉及边界的问题，经纬度什么时候设置好都行；但是如果涉及到插值，想求涡度 散度什么的，就会出现少一圈的情况，这个时候就需要在计算的时候设置一个比较大的经纬度范围，画图的时候再设置成自己要的那个范围