NCL风廓线作图提问

aycc90

求助按照官网的风廓线脚本作图，底部出现红色的不知道是什么符号，很奇怪，附上官网脚本，求大神帮忙解答
;***********************************************; barb_5.ncl;***********************************************;; These files are loaded by default in NCL V6.2.0 and newer; 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  ymin = 2300.            ; min height to be plotted  ymax = 5500.            ; max height to be plotted  ncol = 7                ; number of columns of data (Each ascii file has the same number of columns)  fils = systemfunc("ls barb*.txt")  nfiles = dimsizes(fils)    wks = gsn_open_wks("png","barb")     ; send graphics to PNG file  gsn_define_colormap(wks,"WhBlGrYeRe")  do gg = 0,nfiles-1     t = asciiread(fils(gg) , -1 , "float")     nlvl = dimsizes(t)/ncol           ; figure out number of levels in the ascii file        delete(t)     TestData = asciiread(fils(gg) ,(/nlvl,ncol/), "float")     z    = TestData (:,2)     wspd = TestData (:,6)     wdir = TestData (:,5)     delete(TestData)     rad = 4.0*atan(1.0)/180.     u = -wspd*sin(rad*wdir)     v = -wspd*cos(rad*wdir)             if (gg.eq.0) then        res = True        res@gsnDraw = False        res@gsnFrame = False        res@vpWidthF = 0.8        res@vpHeightF = 0.5        res@vpXF = 0.1        res@vpYF = 0.9        res@trYMinF = ymin        res@trYMaxF = ymax        res@tiMainString = "Wind Profile"        xc = (/-0.5,nfiles-0.5,nfiles-0.5,-0.5,-0.5/)        yc = (/ymin,ymin,ymax,ymax,ymin/)             plot = gsn_csm_xy(wks,xc,yc,res)                          ; draw blank xy plot by outlining min/max of X/Y axes        contours = ispan(5,65,5)*1.                               ; set contours used to define barb color         colors = (/20,30,40,45,50,55,60,65,70,75,80,85,90,95/)    ; set colors     end if        wmsetp ("wdf", 1)                 ; meteorological dir     wmsetp("wbs",0.025)     wmsetp ("blw",2.0)    ; increase line thickness of wind barbs     do hh = 0,dimsizes(z)-1        wmsetp("col", GetFillColorIndex(contours,colors,wspd(hh)))   ; assign barb color based on wind speed value         if (z(hh).lt.ymax.and.z(hh).ge.ymin) then           wmbarb(wks, gg*1., z(hh), u(hh), v(hh))                   ; draw wind barb        end if     end do     delete([/z,wspd,wdir,u,v/])  end do;--------------------------------------------------------------------------;    Draw Labelbar;--------------------------------------------------------------------------   lbres                      = True  lbres@lbPerimOn            = False               ; no label bar box  lbres@lbOrientation        = "Horizontal"        ; orientation  lbres@vpWidthF             = 0.5                 ; size  lbres@vpHeightF            = 0.075  lbres@lbLabelFontHeightF   = 0.015               ; label font height  lbres@lbLabelAlignment     = "InteriorEdges"     ; where to label  lbres@lbMonoFillPattern    = True                ; fill sold  lbres@lbFillColors         = colors  gsn_labelbar_ndc (wks,dimsizes(contours)+1,sprintf("%3.0f",contours),0.25,0.32,lbres)    draw(plot)                        ; draw the plot  frame(wks)                        ; advance the frameend

Eudora99

查看缺省值