;***************************************
; rose_3.ncl
;***************************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" 
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/wind_rose.ncl"
;***************************************

begin
;***************************************
; generate test data and add attributes
;***************************************
;  wrData = wr_GenBogusData (200)
;  wspd   = wrData(0,:)
;  wdir   = wrData(1,:)
                                  
;Read data in
;**************************************************
  data_wspd = asciiread("wspdjn10.txt",-1,"integer")
  data_wdir = asciiread("wdirjn10.txt",-1,"integer")

  npts_s = dimsizes(data_wspd)
  npts_d = dimsizes(data_wdir)
  print(npts_s)
  print(npts_d)

  wspd = data_wspd(:)
  wdir = data_wdir(:)
                                  
  wspd@long_name = "Wind Speed"
  wspd@units     = "m/s"
  wdir@long_name = "Wind Direction"
;******************************

;******************************
; specify plot arguments
;******************************
  numPetals =  8                    ; N, NE, E, SE, S, SW, W, NW
  circFr    = 10. 
;  spdBounds         = (/ 10., 20., 30., 100. /)
  spdBounds          = (/ 1., 2., 5., 10. /)
;******************************
; generate color plot
;******************************
;  wks = gsn_open_wks("x11","rose")
  wks = gsn_open_wks("pdf","rose")
;  wks = gsn_open_wks("ps","rose")
  res               = True
  res@tiMainString  = "Wind Rose: Kigali Mean, January_1000LT"
  colorBounds       = (/ "blue", "green", "yellow", "red" /)
  wrColor           = WindRoseColor (wks,wspd,wdir,numPetals, \
                                     circFr,spdBounds,colorBounds,res)

end