为什么ncl 画气候平均风场 空白

冒失鬼

脚本如下：

begin

  yrStrt = 1980

  yrLast = 2011

  season = "JFM"    ; choose Dec-Jan-Feb seasonal mean

  f      = addfile ("/disk_raid5/data/NCEP2-2.5/monthly/pressure/uwnd.mon.mean.nc", "r")

  TIME   = f->time

  YYYY   = cd_calendar(TIME,-1)/100                 ; entire file

  iYYYY  = ind(YYYY.ge.yrStrt .and. YYYY.le.yrLast)

  uwnd   = short2flt(f->uwnd(11:394,7,0:36,:))           ;

  UWND         = month_to_season (uwnd, season)

  nyrs        = dimsizes(UWND&time)

  time  = TIME(11:394)

  lat=f->lat(0:36)

  lon   = f->lon

  level = f->level

; calculate climatology  

  UWND_cli = dim_avg_n_Wrap(UWND, 0)

ff      = addfile ("/disk_raid5/data/NCEP2-2.5/monthly/pressure/vwnd.mon.mean.nc", "r")

TIME2   = ff->time

YYYY2   = cd_calendar(TIME2,-1)/100                 ; entire file

iYYYY2  = ind(YYYY2.ge.yrStrt .and. YYYY2.le.yrLast)

vwnd   = short2flt(ff->vwnd(11:394,7,0:36,:))           ;VWND         = month_to_season (vwnd, season)

VWND         = month_to_season (vwnd, season)

nyrs        = dimsizes(VWND&time)

time  = TIME(11:394)

lat=ff->lat(0:36)

lon   = ff->lon

level = ff->level

; calculate climatology

VWND_cli = dim_avg_n_Wrap(VWND, 0)

wks = gsn_open_wks("ps","UV300")

gsn_define_colormap(wks,"temp_diff_18lev")       ; choose colormap

  

  

res2 = True

res2@gsnDraw=False

res2@gsnFrame=False

res2@vcRefMagnitudeF         = 2.5             ; define vector ref mag

res2@vcRefLengthF            = 0.05            ; define length of vec ref

res2@vcGlyphStyle            = "LineArrow"    ; turn on curly vectors

res2@vcMinDistanceF          = 0.025

res2@vcRefAnnoFontHeightF    = 0.015

res2@vcRefAnnoSide    = "Top"

res2@vcRefAnnoPerimOn    = False

res2@vcRefAnnoString1  = "2.5m/s"

res2@vcRefAnnoString2On  = False

res2@vcRefAnnoOrthogonalPosF   =  -0.16

plot  =gsn_csm_vector(wks,UWND_cli,VWND_cli,res2); panel plot only resources

  

  

end

运行没有报错，请各位大神看看哪里出问题啦

songwei

最后在end前加两句
draw(plot)
frame(wks)

冒失鬼

songwei 发表于 2017-9-27 09:11
最后在end前加两句
draw(plot)
frame(wks)

我加了之后，就报错：
fatal:VectorPlotDraw: VVECTR - VECTOR NDC LENGTH TOO GREAT
fatal:VectorPlotDraw: error drawing vectors
fatal:VectorPlotDraw: draw error
fatal:PlotManagerDraw: error in plot draw
fatal:_NhlPlotManagerDraw: Draw error

冒失鬼

songwei 发表于 2017-9-27 09:11
最后在end前加两句
draw(plot)
frame(wks)

而且还是空白..

songwei

冒失鬼 发表于 2017-9-27 10:56
而且还是空白..

看看第几行出错

冒失鬼

songwei 发表于 2017-9-27 18:36
看看第几行出错

画出来了，我把res2@gsnDraw=False
res2@gsnFrame=False 注释掉；又把参考值改了一下就好了。谢谢～