;**********************************
; taylor_1.ncl
;**********************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"   
load "./taylor.ncl"
;**********************************
begin
;**********************************
; create some dummy data. all values must be between 0 and 1.65
; note on the resulting plot, the x-axis has no tickmark labels.
;**********************************
; "p" dataset                      
  p_xpts    = (/0.60, 0.50, 0.45, 0.75, 1.15/)
  p_ypts    = (/1.24, 0.75, 1.00, 0.93, 0.37/)

; "t" dataset
  t_xpts    = (/0.75, 0.64, 0.40, 0.85, 1.15 /)
  t_ypts    = (/0.24, 0.75, 0.47, 0.88, 0.73/)

  nDataSets  = 2                               ; number of datasets
  npts       = dimsizes(p_xpts)
  x          = new ((/nDataSets, npts/),typeof(p_xpts) )
  y          = new ((/nDataSets, npts/),typeof(p_ypts) )

  x(0,:) = p_xpts
  y(0,:) = p_ypts

  x(1,:) = t_xpts
  y(1,:) = t_ypts
;**********************************
; create plot
;**********************************
  res   = False                             ; default taylor diagram

  wks   = gsn_open_wks("ps","taylor") 
  plot  = taylor_diagram(wks,x,y,res)

end