Ncl--垂直速度叠加风场的纬向剖面图
作者:第八星系-木零
邮箱:1692089531@qq.com
数据:nc数据
f=addfile("E:\data\fnl_2022022700_00.nc","r")
;读取数据 (time,lev,lat,lon)
lat =f->lat_0
lon =f->lon_0
lev =f->lv_ISBL0
vv =f->VVEL_P0_L100_GLL0 ({85000:20000}, {36.53}, {96:108})
U =f->UGRD_P0_L100_GLL0({80000:21000}, {36.53}, {97:107})
V =f->VGRD_P0_L100_GLL0({80000:21000}, {36.53}, {97:107})
wks=gsn_open_wks("png","vv")
;*******************************
res=True
res@gsnMaximize = True ;最大化
res@gsnDraw = False ;暂不输出
res@gsnFrame = False
res@cnFillOn = True ;颜色填充
res@cnFillPalette = "NCV_blu_red";色板
res@cnLinesOn = False ;关闭等值线
res@cnSmoothingOn = True ;平滑
res@cnSmoothingTensionF = -2.5
res@gsnPresHgtHeightLabelOn = False
res@tiYAxisString = ""
res@gsnRightString = ""
res@gsnLeftString = ""
;*******************************
;刻度
res@tmXTOn = False
res@tmXBMode = "Explicit"
res@tmXBValues = ispan(98,106,4)
res@tmXBLabels =(/"98~S~o~N~E","102~S~o~N~E","106~S~o~N~E"/)
res@tmXTMinorOn = False
res@tmYLMajorLengthF = 0.007 ;Y轴左刻度线长度
res@tmYRValues = False ;关掉右侧值
res@tmXBMinorOn = False ;去掉小刻度
;*******************************
;色标
res@pmLabelBarHeightF = 0.1 ;宽度
res@pmLabelBarWidthF = 0.69 ;长度
res@pmLabelBarParallelPosF = 0.6 ;水平位置
res@pmLabelBarOrthogonalPosF = 0.001 ;垂直位置
res@lbLabelFontHeightF = 0.0179 ;字体大小
res@lbTitleOn = True
res@lbTitleString = "Pa/s" ;色标单位
res@lbTitleFontHeightF = 0.0176 ;字体大小
res@lbTitleDirection = "Across" ;横向显示
res@lbTitlePosition = "Right" ;位置
res@lbTitleJust = "CenterLeft"
;*******************************
;等值线
res@cnLevelSelectionMode = "ManualLevels" ;等值线设置
res@cnMinLevelValF = -0.8 ;最小值
res@cnMaxLevelValF = 0.8 ;最大值
res@cnLevelSpacingF = 0.1 ;间隔
res@cnInfoLabelOn = False ;关闭等值线信息显示
;*******************************
;风羽
res@vcRefMagnitudeF = 10. ; 风羽大小
res@vcRefLengthF = 0.029 ; 风羽长度
res@vcGlyphStyle = "WindBarb" ; 类型
res@vcMinDistanceF = 0.03 ; 风羽上下间的距离
res@vcWindBarbScaleFactorF = 2.5
res@vcRefAnnoOn = False ; 关闭参考风矢
res@vcWindBarbTickLengthF = 0.4 ; 风羽的长短(默认0.3)
res@vcWindBarbLineThicknessF = 2. ; 风羽粗细
vv = smth9_Wrap(vv,0.5,-0.25,False)
plot= gsn_csm_pres_hgt_vector(wks,vv,U,V,res)
draw(plot)
frame(wks)
微信搜索“第八星系人造大气理论爱好者”公众号,关注获取文章数据
页:
[1]