ncl 风矢量图组
作者：第八星系-欣妹儿邮箱：3035245582@qq.com
利用uwnd.mon.mean.nc、 vwnd.mon.mean.nc数据，绘制2001-2010年200hPa (60°S-60°N、0°E -180°E)范围内4月和10月平均风速分布（即二维风矢量图）并放置在一个图组中，即一张图分别包含上图（4月二维风矢量图）和下图（10月二维风矢量图）两张子图。并根据所学内容将图画的更美观。
读取数据
begin
f1=addfile("/cygdrive/d/ncl/experiment-ncl/5shangji/data/uwnd.mon.mean.nc","r")
f2=addfile("/cygdrive/d/ncl/experiment-ncl/5shangji/data/vwnd.mon.mean.nc","r")
YYYYMM1=cd_calendar(f1->time,-1) ;将时间转为年月日
YYYYMM2=cd_calendar(f2->time,-1)
计算UV月平均
uwnd1=clmMonTLL(f1->uwnd(:,{200},{-60:60},{0:180}))
vwnd1=clmMonTLL(f2->vwnd(:,{200},{-60:60},{0:180}))
; printVarSummary(uwnd1)
画图
wks=gsn_open_wks("png","/cygdrive/d/ncl/experiment-ncl/5shangji/5-2")
res=True
res@gsnDraw=False
res@gsnFrame=False
;res@vcRefAnnoAngleF=0.05
res@vcRefMagnitudeF=20
res@vcRefLengthF=0.045
res@vcMinDistanceF=0.017
res@vcGlyphStyle="CurlyVector"
res@gsnAddCyclic=False
res@mpMinLatF=-60
res@mpMaxLatF=60
res@mpMinLonF=0
res@mpMaxLonF=180
;创建组图
plot=new(2,graphic)
res@gsnLeftString="Apr.wind"
plot(0)=gsn_csm_vector_map(wks,uwnd1(3,:,:),vwnd1(3,:,:),res)
res@gsnLeftString="Sep.wind"
plot(1)=gsn_csm_vector_map(wks,uwnd1(9,:,:),vwnd1(9,:,:),res)
pres=True
pres@txString=" wind velocity vector"
gsn_panel(wks,plot,(/2,1/),pres)
end
微信搜索“第八星系人造大气理论爱好者”公众号，关注获取文章数据
页:
[1]