oacres循环插值问题

zytcerulean

站点观测的降水资料储存格式为
站点（stn） 纬度（lat） 经度（lon） 降水量（r）
共有2001年6月30天的资料 储存为grd格式 分别为20010601 20010602......20010630
已生成.map文件 10092002.map
现在需要将30天的降水插值图画出并输出插值后各个格点上的降水量
*******************grid.ctl***********为格点信息
dset D:\extreme\climate01-10\6grd\gridt500.dat
title data file to draw china map
undef 0
xdef 116 linear 70.3125  0.5625
ydef 73 linear 14.625  0.5625
zdef 1 linear 1 1
tdef 1 linear 17Z13Jul2009 12hr
vars  1
prc 0 99 afd
ENDVARS
**********************20110724.tcl********************
dset D:\extreme\climate01-10\6grd\200106%d2.grd
options template
dtype station
stnmap D:\extreme\climate01-10\6grd\10092002.map
undef -999
title rain
tdef 1 linear 17Z13Jul2009 12hr
vars 1
hight  0 99  station hight
endvars
*********************.gs文件*******************
'reinit'
'open D:\extreme\climate01-10\6grd\grid.ctl'
'open D:\extreme\climate01-10\6grd\20110724.ctl'
day=1
t=1
while(day<31)
'enable print D:\extreme\climate01-10\6grd\r200106'%day%'.gmf'
'set rgb 20 153 255 153'
'set rgb 30  51 204  51'
'set rgb 40 102 204 255'
'set rgb 50   0   0 255'
'set rgb 60 255   0 255'
'set rgb 70 150   0   0'
'set rgb 22 255 232 120'
'set rgb 26 255  50   0'
'set rgb 48  30 110 235'
'set t ' t

'set mpdset cnworld'
'set vpage 0 11 0 8.5'
'set strsiz 0.25'
'set string 1 c'
'set ccolor rainbow'
'set grads off'
'set grid off'
'set parea 1.5 9.5 0.8 7.8'

'set xlabs  73|81|90|99|108|117|126|135`3.`1E'
'set ylabs  18|24|30|36|42|48|54`3.`1N'
'set xlopts 1 1 0.225'
'set ylopts 1 1 0.225'
'set clab on'
'define A=MASKOUT(OACRES(prc,hight.2,10,8,6,3,1,1,1),prc(t=1)-0.1)'
'set gxout shaded'
'set clevs 5 10 25 30 40 50 75 100 '
'set ccols 0 20 30 22 40 48 50 60 70'
'd A'
'cbarn 1 0 '
'set gxout contour'
'set clevs 5 10 15 25 30 40 50 75 106 '
'd A'
'print'
'disable print'

day=day+1
t=t+1
endwhile
;
现在，运行完.gs文件出现错误提示 Warning from oacres:less than two stations
grid will be all missing values
can not contour grid——all undefined values
can not plot color bar：no shading information
can not contour grid——all undefined values
之前用一天的grd文件画一幅图是正确的 但是循环起来就出现问题……请大家帮忙看一下是什么原因 谢谢啦

mofangbao

有人出现过类似的，自己论坛搜索下，给过查错脚本的

river

楼主看下这个帖子http://bbs.06climate.com/forum.php?mod=viewthread&tid=10006&extra=&page=1 ，参考一下

zytcerulean

谢谢 2l的帖子很有帮助~我的思路有问题

meifhaisxiac

我也出现了“Warning from oacres:less than two stations”的错误，原来的写法是
'set z 6'
'DEFINE a=OACRES(g(z=6),u.2(z=6)) '
之后改成
'DEFINE a=OACRES(g(z=6),u.2(z=6)) '
就没有出错了，貌似有多层变量的时候，高度层的信息需要在变量名后面直接指定grads才认得。

samvelly

{:eb502:}{:eb502:}

quilowlin

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