求教GrADS站点格式数据cressman插值后输出的问题

lqouc

在处理多变量多时次的站点数据时想要利用GrADS将站点数据插值（按照清风大大的教程），绘制了填色图，都没问题，但是想要把差值出来的格点数据输出，结果得到的数据全是0.请高手指点，代码如下：
'set gxout fwrite'
'set fwrite d:\station\14\14avevis.dat'
'open d:\station\14\grid.ctl'
'open d:\station\14\14avevis.ctl'
'set lon 70 140'
'set lat 15 55'
t=1
while(t<=131)
'set t 't''
'set mpdset hires'
'define a1=oacres(g,vis.2)'
'define a2=oacres(g,tem.2)'
'define a3=oacres(g,td.2)'
'define a4=oacres(g,wd.2)'
'define a5=oacres(g,ws.2)'
'd a1'
'd a2'
'd a3'
'd a4'
'd a5'
t=t+1
endwhile
'disable fwrite'
'reinit'

river

是不是和数据输出的精度设置有关呢

lqouc

river 发表于 2013-1-9 18:12

                               
登录/注册后可看大图

是不是和数据输出的精度设置有关呢

数据的输出时1*1的，跟插值时那个网格场的是一样的，没有地方设置，我觉得应该也不用设置吧。

river

lqouc 发表于 2013-1-9 18:16

                               
登录/注册后可看大图

数据的输出时1*1的，跟插值时那个网格场的是一样的，没有地方设置，我觉得应该也不用设置吧。

我说的不是网格的精度，是说数据的精度，就是保留几位小数点那个意思。有个命令可以设置，但是我一时想不起来了

lqouc

river 发表于 2013-1-9 18:40

                               
登录/注册后可看大图

我说的不是网格的精度，是说数据的精度，就是保留几位小数点那个意思。有个命令可以设置，但是我一时想不 ...

额，好吧，我也不知道为什么，以前转nc的数据都很顺利，那我去查查。

lqouc

虚心求教啊……问题仍然没解决

空谷幽竹

是不是Cressman插值半径的问题呢，楼主可以试试oacres(g,vis.2,1.0)

shangww

想请问下这个R是什么意义

皮特

你好，你解决了吗？我也遇到这个问题了，插值后转的dat，画出来的图带e的！

lqouc

皮特 发表于 2015-11-17 14:13
你好，你解决了吗？我也遇到这个问题了，插值后转的dat，画出来的图带e的！

问题肯定是解决了，不过时间太久远了，我也记不住了。
可以详细描述下你的问题，带上ctl，gs和图。
从当前的描述来看ctl里面的undef应该是-9.99e8，不知道这个你设置的对不对。如果改了还不行，看考我第二行的话。