ncl多年平均的风场剖面图组
作者：第八星系-欣妹儿
邮箱：3035245582@qq.com
利用uwnd.mon.mean.nc数据，绘制沿45°N纬线和沿95°E经线的多年平均U风的垂直剖面等值线填色图（既画等值线又画填色图）并放置在一个图组中，即一张图分别包含左图（沿45°N的垂直剖面）和右图（沿95°E的垂直剖面）两张子图。设置等值线等级范围为-10m/s~30m/s，等值线间隔为5m/s，标注子图序号，设置公共色标等使图形更美观。
完整代码
lat1=45
lon1=95
begin
f=addfile("/cygdrive/d/ncl/experiment-ncl/5shangji/data/uwnd.mon.mean.nc","r")
YYYYMM=cd_calendar(f->time,-1)
;print(YYYYMM)
uwnd1=f->uwnd(:,:,{lat1},:)
uwnd2=f->uwnd(:,:,:,{lon1})
z0=dim_avg_n_Wrap(uwnd1,0)
z1=dim_avg_n_Wrap(uwnd2,0);从0维开始 时间是第0维
;vardp0=stdatmus_z2tdp(z0)
;vardp1=stdatmus_z2tdp(z1)
;copy_VarCoords(z0,vardp0(0,:))
;copy_VarCoords(z1,vardp1(0,:))
;printVarSummary(z0)
wks=gsn_open_wks("png","/cygdrive/d/ncl/experiment-ncl/5shangji/5-1")
res1=True
res1@gsnDraw=False
res1@gsnFrame=False
res1@gsnLeftString="U_wind"
res1@cnLevelSelectionMode="ManualLevels"
res1@cnLevelSpacingF = 5.0
res1@cnMaxLevelValF=30.0
res1@cnMinLevelValF=-10.0
res1@cnLineLabelsOn=True
res1@cnFillOn=True
res1@cnFillPalette="MPL_rainbow"
;创建图组
plot=new(2,graphic)
plot(0)=gsn_csm_pres_hgt(wks,z0,res1)
plot(1)=gsn_csm_pres_hgt(wks,z1,res1)
pres=True
pres@txString="Mean wind speed vertical profile"
gsn_panel(wks,plot,(/1,2/),pres)
end
微信搜索“第八星系人造大气理论爱好者”公众号，关注获取文章数据
页:
[1]