求教！！任意直线画剖面分辨率问题

diandian

参考论坛的帖子画任意一条直线的垂直剖面，但是
while(lon<=lon2)
  lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
  'collect 1 gr2stn(dbz,'lon','lat')'
  lon =lon+0.05
endwhile
精度分辨率取不同值时，有的可以画出图，有的则不能画出图，如0.05、0.2可以出图，但0.06、0.5、。。。。很多都画不出图，且取0.2和0.05市画出的图也不一样，如下：第一张0.2，第二张0.05

其实，自己也模模糊糊懂一点，就是分辨率不同，第一张0.2导致下一个值为101.5+0.2=101.7，所以第一张没有画出第二张101.7之前的图，但是第二张是0.05，为什么没有画出第一张101.7之后的图呢？
还有就是经度分辨率如何取值，上面的值都是试出来的，总感觉这样画出来不是很精确，且0.05可以，0.06就不行，这都是为什么啊？？
请教大神解答一二，不胜感激！！

diandian

完整gs如下：
'reinit'
'open d1.ctl'
'set grads off'
'set xlopts 1 1 0.15'
'set ylopts 1 1 0.15'
'set zlog on'
*'set csmooth on'

'set t 17'
'set x 1'
'set y 1'
*'set lon 111 118'
*'set lat 30 33'
'set z 1 19'
lon1 = 101.5
lon2 = 102.2
lat1 = 36.65
lat2 = 35.25
lon = lon1
'collect 1 free'
while(lon<=lon2)
  lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
  'collect 1 gr2stn(dbz,'lon','lat')'
  lon =lon+0.05
endwhile

'set cmin 0.01'
'set t 17'
'set lon 101.5 102.2'
'set xlint 0.1'
'set z 1 19'
*'set xaxis 'lon1' 'lon2
'set clab on'
'set gxout shaded'
*'set ccols 0 0 7 0'
'd coll2gr(1,-u)'
'cbarn_interp.gs'

'set cmin 0.01'
'set gxout contour'
'set cthick 1'
'd coll2gr(1,-u)'

'printim qxiepo.png white'

'reinit'

river

这个最好根据资料的分辨率来设置吧，我觉得

diandian

river 发表于 2015-5-18 21:10
这个最好根据资料的分辨率来设置吧，我觉得

资料是wrf 经ARWpost后处理得到的dat格式数据，ctl中
xdef  967 linear   68.18367   0.06756757
ydef  552 linear   15.68879   0.06756757
但是，这个分辨率也试过，画不出图，还是不知为什么？？

river

diandian 发表于 2015-5-18 22:25
资料是wrf 经ARWpost后处理得到的dat格式数据，ctl中
xdef  967 linear   68.18367   0.06756757
ydef  ...

画不出图有没有什么具体的提示？

river

diandian 发表于 2015-5-18 22:25
资料是wrf 经ARWpost后处理得到的dat格式数据，ctl中
xdef  967 linear   68.18367   0.06756757
ydef  ...

也可以去看一下官方文档http://grads.iges.org/grads/gadoc/gadocindex.html

diandian

river 发表于 2015-5-19 08:44
画不出图有没有什么具体的提示？

没有报错，今天又试了一下， 0.06756757分辨率又可以了，我同学也遇到这样，就像看运气一样~~

diandian

river 发表于 2015-5-19 08:44
画不出图有没有什么具体的提示？

好吧，刚说完可以，现在0.06756757又不可以了

river

diandian 发表于 2015-5-19 10:18
好吧，刚说完可以，现在0.06756757又不可以了

你把你的grads升级到2.1再试试

diandian

river 发表于 2015-5-19 10:36
你把你的grads升级到2.1再试试

我发现再重新打开grads之后就可以了~~在服务器上画的，感觉也可能是版本问题