;***********************************************
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   = -500              ; min height to be plotted
  ymax   = 5000              ; max height to be plotted
  ncol   = 42                ; number of columns of data (Each ascii file has the same number of columns)
  wbcrit = 0                 ; minimum distance in meters between wind barbs 

  fils = systemfunc("ls *.csv")
  nfiles = dimsizes(fils)
  
  wks = gsn_open_wks("pdf","barb05_dec")
  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")

     do j=0,ncol-3,3
   	z    = TestData (:,j)
   	wdir = TestData (:,j+1)
   	wspd = TestData (:,j+2)
     end do

     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.7
        res@vpHeightF = 0.45
        res@vpXF = 0.15
        res@vpYF = 0.9
        res@trYMinF = ymin
        res@trYMaxF = ymax
	res@tmXBMode = "Explicit"
        res@tmXBValues = ispan(0,nfiles-1,1)
        res@tmXBLabels = (/"11Z Dec","17Z Dec","23Z Dec","05Z Dec"/)
        res@tiYAxisString = "height (m)"
        res@tiXAxisString = "time (h)"
        res@tiMainString = "Wind Vertical Profile at Cotonou in December, 1999"
        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)
     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 frame
end