【求助】企图用fnl资料画925散度场，并叠加micaps第3类资料的6小时降水量，出...

clare

说明：p.s.主题分类选择不了  脚本编辑
我借鉴论坛中好几个画micaps第3类资料的帖子，主要为@luoziwuhui的http://bbs.06climate.com/forum.php?mod=viewthread&tid=13648
尝试grads画图。
然后想模仿@传说中的谁 http://bbs.06climate.com/forum.php?mod=viewthread&tid=893帖子最后的出图效果。
本人使用了36个时次的micaps3的6小时降水资料；生成的格点场为0.25x0.25;  没有使用超赞的meteoinfo，因为对底图不需要做修饰。
但是一旦叠加散度场就会成酱紫：
如图：


单独画降水：


对应的原始数据：
      

对应的ctl:
         

画图使用的gs:


'reinit'
'open ..\fnl.ctl'
'open rain6h.ctl'
'open grid.ctl'
'enable print r6.gmf'
********************************************
i=1
while(i<=28)
'set lat 20 44'
'set lon 100 124'
'set t 'i''
'set mpdset cnworld cnriver'
'set parea 1 10 1 8'
'set xlopts  1 4 0.20'
'set ylopts  1 4 0.20'
'set clopts -1 -1 0.13'
'set xlab off'
'set ylab off'
'set grads off'
'set grid off'
'set csmooth on'
'set gxout shaded'
'set cmin 10'
'set clevs  10 15 20 25 30 35 40 50 60 70 80 90 100 120'
'set rgb 20 255 255 255'
'set rgb 21 205 236 251'
'set rgb 22 150 212 243'
'set rgb 23 106 173 220'
'set rgb 24 72 148 181'
'set rgb 25 113 195 179'
'set rgb 26 73 168 113'
'set rgb 27 116 193 75'
'set rgb 28 206 219 86'
'set rgb 29 249 195 77'
'set rgb 30 246 130 51'
'set rgb 31 232 78 41'
'set rgb 32 215 39 40'
'set rgb 33 181 26 33'
'set rgb 34 146 21 25'

'set ccols 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 '
'set xlpos -20'                                                                           
'set ylpos -20'
'cnbasemap  oacres(g.3,r.2)'
*'cnbasemap  oacres(g.2,r.1,2,50,20,10,5,1)'


'set gxout contour'
'set cthick 5'
'set ccolor rainbow'
'set csmooth on'
'set cint 5'
'set lat 20 44'
'set lon 100 124'
'set lev 925'
'define div= hdivg(UGRDprs,VGRDprs)*1e6'
'd div'



'run axis.gs -type b  -position o  -lfont 5 -lsize 0.15 -lthick 0.6 -interval 4 -sinterval 4 -suffix `3.'
'run axis.gs -type L  -position o  -lfont 5 -lsize 0.15 -lthick 0.6 -interval 4 -sinterval 4 -suffix `3.'  
*'cbarm 0.8 1 9.8  4.8'

'q dims'
rec=sublin(result,5)
recs=subwrd(rec,6)
' draw title Micaps'recs'-6h rain'

'print'
'c'
i=i+1
endwhile
'disable print'
;
错误信息提示：

river

可能是格点场分辨率的问题，你的fnl资料是1*1的，格点场也生成1*1的试一下

造字仓颉

非常感谢楼主，希望好好学习一下，提高一下技术。

dych

我是来学习的

abaike

我是来学习的