;***********************************************
; 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 = 16.            ; min height to be plotted
  ymax = 10767.         ; max height to be plotted
  ncol = 4                ; 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
     gg=1
     t = asciiread("WS.txt", -1 , "float")
                                          ; print(t)
                                          ; print(dimsizes(t))
     nlvl = dimsizes(t)/ncol              ; figure out number of levels in the ascii file   
            print(nlvl)
     delete(t)
     TestData = asciiread("WS.txt" ,(/nlvl,ncol/), "float")
;    printVarSummary(TestData)
;     printVarSummary(t)
     z    = TestData (:,1)
                                         ;print(z)
     wdir = TestData (:,2)
     wspd = TestData (:,3)

     wdir@_FillValue = 999
     wspd@_FillValue = 99
;                                         print(z+"  "+wspd+"  "+wdir)
     delete(TestData)

     rad = 4.0*atan(1.0)/180.
     u = -wspd*sin(rad*wdir)
     v = -wspd*cos(rad*wdir)
                                          print(u+ "    "+v)
        
;    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"


nfiles=1
        xc = (/-20.5,nfiles-0.5,nfiles-0.5,-0.5,-20.5/)
        yc = (/ymin,ymin,ymax,ymax,ymin/)
             print(xc)

     
        plot = gsn_csm_xy(wks,xc,yc,res)                          ; draw blank xy plot by outlining min/max of X/Y axes
        contours = ispan(6,30,2)*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
        if (ismissing(wspd(hh)).eq.True) then 
            print("stop")
         else
        wmsetp("col", GetFillColorIndex(contours,colors,wspd(hh)))   ; assign barb color based on wind speed value
         end if
        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 frame
end