grads处理micaps第一类数据

筱筱筱筱筱雅痴

格点、温度、露点温度的ctl都是正确的，因为可以画出温度场剖面图，但根据下面的脚本画相对湿度剖面图时，一直报错，改来改去都不对，希望大神们帮忙看一下
'open J:\lw-data\grads\jiaoben\grid.ctl'
'open J:\lw-data\grads\jiaoben\temp.ctl'
'open J:\lw-data\grads\jiaoben\dtemp.ctl'
i=1
while(i<=560)
'set grads off'
'set lon 110 120'
'set lat 35 45'
'set t 'i''
'set mpdset hires'
'set map 1 1 9'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'
'define Es=6.11*pow(10.0,7.5*dtemp.3/(237.7+dtemp.3))'
'define E=6.11*pow(10.0,7.5*temp.2/(237.7+temp.2))'
'define rh=(E/Es)*100'
'set gxout shaded'
'cnbasemap oacres(g,rh)'
'cbarn.gs'
'q time'
x=subwrd(result,3)
'draw title 'x''
'printim J:\lw-data\grads\gs\micaps_rh_'x'.png white'
'c'
i=i+1
endwhile
'reinit'
;

river

我觉得插值函数的前提是一定范围内一些点上的数据是已经存在的，而你的rh是在计算过程中的，一个点一点的计算。所以插值函数才会报错吧。你可以试试把rh输出成二进制文件，配上描述文件，再来画应该可以。
楼主也可以去看下这个帖子http://bbs.06climate.com/forum.p ... p;extra=&page=1

筱筱筱筱筱雅痴

river 发表于 2018-4-9 21:23
我觉得插值函数的前提是一定范围内一些点上的数据是已经存在的，而你的rh是在计算过程中的，一个点一点的计 ...

好的，万分感谢，我再来试一下

robot814

计算方法还是可以进一步探讨