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

皮特

lqouc 发表于 2015-11-17 14:38
问题肯定是解决了，不过时间太久远了，我也记不住了。
可以详细描述下你的问题，带上ctl，gs和图。
从 ...

你好，首先感谢你的回复。
我的问题和你类似，如下：
1、插值后用来画图或转二进制的GS
'reinit'
'open F:\Mstation\grid.ctl'
'open F:\Mstation\prec.ctl'
'set gxout fwrite'
'set fwrite F:\Mstation\dailyprecip2015.grd'（转二进制)
*'enable print F:\Mstation\rain6.gmf'（画图）
*n=1
*while(n<=184)

'set x 1 27'
'set y 1 17'
'set z 1'
*'set t 'n''
'set t 1'
'define a=oacres(g,p.2)'
*'define a=oacres(g,p.2,15,7,4,2,1)'
'define a1=maskout(a,g-0.5)'
'define aa=smth9(a1)'
'd a'
*n=n+1
*endwhile
'disable fwrite'
'print'
'reinit'

1、对某日的降水量插值后，直接画图

2、插值后转成二进制grd，再写ctl画图（出现9e+1的结果）
DSET F:\Mstation\dailyprecip2015.grd
title sample grid data
undef -999.0
xdef  27 linear 97.00 1.000
ydef  17 linear 18.00 1.000
zdef 1 linear 1000 1
tdef 184 linear 1jan1951 1yr
vars 1
aa 0 99 grid data
endvars

请问楼主当时怎么解决的，谢谢了！

lqouc

皮特 发表于 2015-11-17 15:51
你好，首先感谢你的回复。
我的问题和你类似，如下：
1、插值后用来画图或转二进制的GS

额，这不就是我上面回复里告诉你的么。
所有的脚本和计算都是对的，只是你的ctl里面undef要改成-9.99e8。
原因是经过grads计算后得到的所有结果的缺测值都会被修改为-9.99e8，如果原始数据没有缺测就看不出来，有的话就会出现你图中那样了。所以以后给输出的数据配ctl的时候注意下undef就行了。

皮特

lqouc 发表于 2015-11-17 16:12
额，这不就是我上面回复里告诉你的么。
所有的脚本和计算都是对的，只是你的ctl里面undef要改成-9.99e8 ...

哥们，你简直大神啊，ok了，感谢！这样把站点数据写成一个格点数据，处理什么都方便！

lqouc

皮特 发表于 2015-11-17 17:39
哥们，你简直大神啊，ok了，感谢！这样把站点数据写成一个格点数据，处理什么都方便！

不用客气~

lujianjun

{:5_213:}{:5_213:}{:5_213:}

zhaohaoran999

学习学习！！！

陀思妥耶夫斯基

膜拜大神！受教了