对micaps站点数据进行插值，并且想在图上标注每个站点的风向风速和计算出来的相对湿度

筱筱筱筱筱雅痴

我的gs脚本：
'open J:\lw-data\grads\jiaoben\grid.ctl'
'open J:\lw-data\grads\jiaoben\rh.ctl'
'open J:\lw-data\grads\jiaoben\wind.ctl'
i=1
while(i<=560)
'set grads off'
'set grid off'
'set lon 112 120'
'set lat 36 43'
'set t 'i''
'set mpdset hires'
'set map 1 1 9'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'
'set gxout shaded'
'cnbasemap oacres(g,rh.2)'
'cbarn.gs'
'set gxout stnmark'
'set digsiz 0.2'
*'set stid on'
'd rh.3'
'set gxout barb'
'set ccolor 1'
'set cthick 8'
'set digsiz 0.1'
'd -sin(winddir.3*3.14159/180)*wind.3*2.5;-cos(winddir.3*3.14159/180)*wind.3*2.5'
'q time'
m=subwrd(result,3)
'draw title 'm''
'set parea 1.5 10 1.5 8'
'q w2xy 117.12 39.04'
x=subwrd(result,3)
y=subwrd(result,6)
'set string 2 tc 8 0'
'set strsiz 0.2 0.2'
'draw string 'x' 'y' Tianjin'
'q w2xy 117.12 39.04'
x=subwrd(result,3)
y=subwrd(result,6)
'draw mark 3 'x' 'y' 0.1'
'printim J:\lw-data\grads\gs\deal\micaps_rh_'m'.png white'
'c'
i=i+1
endwhile
'reinit'
;
脚本中的wind.ctl：
dset J:/lw-data/xlw/Micaps/diamond1_batch/wind.dat
dtype station
stnmap J:/lw-data/xlw/Micaps/diamond1_batch/wind.map
undef 9999
title relative humidity
tdef 560 linear 02z01Nov2016 3hr
vars 3
rh 0 99 xiangduishidu
wind 0 99 fengsu
winddir 0 99 fengxiang
endvars

筱筱筱筱筱雅痴

有人知道嘛，教教我

micolian

楼主， 有个软件叫 “迈特因否”，我音译的，具体的你自己去找，可以直接对MICAPS数据进行画图的

四叶草

看楼主的脚本，湿度设置成填色显示，这样具体的数值肯定标不上去嘛