NCL画近10年的北半球平均纬向风图

rolemodel

脚本如下，但是会报错，我实在改不出来，大神帮忙看看啊load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
;load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
itimes = 200701
itimee = 201701
fi = addfile("/cygdrive/d/yanjiusheng/uwnd.mon.mean.nc","r")
var = fi->uwnd
time = fi->time
YYYYMM := cd_calendar(time,-1)
ist = ind(itimes.eq.YYYYMM)
ied = ind(itimee.eq.YYYYMM)
uwnd = var(ist:ied,:,:,:)
uwnd_zm = dim_avg_Wrap(uwnd)
wks = gsn_open_wks("png","/cygdrive/c/Data/zuoye")
res = True
res@gsnDraw = False
res@gsnFrame = False
res@gsnLeftString = ""
res@gsnRightString =  ""

res@tmXBMode = "Explicit"
res@tmXBValues = (/-90,-60,-30,0,30,60,90/)
res@tmXBLabels = (/"90~S~o~N~S","60~S~o~N~S","30~S~o~N~S","30~S~o~N~N","60~S~o~N~N","90~S~o~N~N"/)
res@tmXBLabelDeltaF = -0.6
res@tmXBLabelFontHeightF = 0.015
res@tmYLLabelFontHeightF = 0.015
res@tiYAxisString = "Level(hPa)"
res@cnFillOn = False
res@cnLinesOn = True
res@cnLineLabelsOn = False
res@cnLevelSpacingF = 0.05
res@gsnContourZeroLineThicknessF = 0
res@cnInfoLabelOn = False
res@lbLabelBarOn = False
res@cnFillOn = False
res@cnLinesOn = True
res@cnLevelSpacingF = 0.15
res@cnLineDashPattern = 2
res@cnLineThicknessF = 3.
res@cnInfoLabelOn = True

plot = gsn_csm_pres_hgt(wks,var(level|:,lat|:),res)
plot = ColorNegDashZeroPosContour(plot,"black","white","black")
draw(plot)
frame(wks)
end

南气斑织逍遥

uwnd = var(ist:ied,:,:,:)
uwnd_zm = dim_avg_Wrap(uwnd)

没定义层数，函数用法不对，看官网

rolemodel

南气斑织逍遥 发表于 2017-9-27 12:10
uwnd = var(ist:ied,:,:,:)
uwnd_zm = dim_avg_Wrap(uwnd)

大神，我想画那种左边坐标是level，X坐标是纬度的平均纬向风的图。需要定义层数吗？

南气斑织逍遥

rolemodel 发表于 2017-9-27 20:39
大神，我想画那种左边坐标是level，X坐标是纬度的平均纬向风的图。需要定义层数吗？

uu=dim_avg_n_wrap(u,(/0/))
uuu=dim_avg_n_wrap(uu,(/2/))
对时间维和经向做平均
画图的var(level|:,lat|:)时错误写法，你没有定义这些名字时什么，直接写uuu加你想要的时间维就可以

rolemodel

南气斑织逍遥 发表于 2017-9-28 05:01
uu=dim_avg_n_wrap(u,(/0/))
uuu=dim_avg_n_wrap(uu,(/2/))
对时间维和经向做平均

谢谢大神，已解决。学习了

小悦悦

您好，我想请问一下你用ncl画近10年的北半球平均纬向风图那个程序最后是怎么解决的？能把脚本发来看一下吗，谢谢了