expansion half circle of Taylor Diagram,

From: antonio avio <antonioavio_at_nyahnyahspammersnyahnyah>
Date: Wed Jul 04 2012 - 16:25:13 MDT

HI all,
I just wonder, anybody has a modified Taylor diagram with half circle
covering negative
correlation for assessing the models skill?

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "/home/NCL/Taylor_Diagram/taylor_diagram_cam.ncl"
;**********************************
begin
;**********************************
; Assume the following have already been computed:
; _ratio are the ratio: Case_Variance/Reference_Variance
; _cc are the cross correlation coef of Case to Reference
;**********************************
; Cases [Model]
  case = (/ "Data Assimilation"/)
  nCase = dimsizes(case ) ; # of Cases [Cases]

; variables compared
  var = (/ "F=1","F=2","F=3","F=4","F=5","F=6","F=7","F=8","F=9"/)
  nVar = dimsizes(var) ; # of Variables

; "Case A"
  CA_ratio =
(/1.4267,0.5652,1.3196,0.5984,0.5370,0.3953,0.5224,1.9355,4.2984/)
  CA_cc =
(/0.2175,-0.1005,0.09408,-0.0853,0.007545,0.004018,-0.03651,-0.09119,-0.02975/)
  BA = (/2.7 ,-1.5 , 17.31,-20.11, 12.5 , 8.341, 5, 2, 3 /)
; arrays to be passed to taylor plot
  ratio = new ((/nCase, nVar/),typeof(CA_cc) )
  cc = new ((/nCase, nVar/),typeof(CA_cc) )
  bias = new ((/nCase, nVar/),typeof(CA_cc) )

  ratio(0,:) = CA_ratio
  cc(0,:) = CA_cc
  bias(0,:) = BA

;**********************************
; Put the ratios and pattern correlations into
; arrays for plotting
;**********************************

  nDataSets = 1 ; number of datasets
  npts = dimsizes( CA_ratio)
  ratio = new ((/nDataSets, npts/),typeof(CA_ratio) )
  cc = new ((/nDataSets, npts/),typeof(CA_cc ) )

  ratio(0,:) = CA_ratio
  cc(0,:) = CA_cc

 ; ratio(1,:) = t_rat
 ; cc(1,:) = t_cc

;**********************************
; create plot
;**********************************

  res = True ; default taylor diagram

  res@Markers = (/10, 10/) ; make all solid fill
  res@Colors = (/"red", "blue" /)
  res@varLabels = var
  res@caseLabels = case
  res@gsMarkerSizeF = 0.0060

res@varLabelsYloc=1.45
res@caseLabelsFontHeightF = 0.09 ; make slight larger [default=0.12
]
res@varLabelsFontHeightF = 0.009 ; make slight smaller
 [default=0.013]
res@centerDiffRMS = True
;res@ccRays_color = "LightGray" ; default is black
 res@stnRad = (/ 0.5, 1.5 /) ; additional standard radii
 res@ccRays = (/ 0.6, 0.9 /) ; correlation rays

  wks = gsn_open_wks("pdf","taylor")
  plot = taylor_diagram_cam(wks,ratio,cc,bias,res)
  frame(wks)

end

Thanks a lot for your help,
Antonio.

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jul 4 16:25:26 2012

This archive was generated by hypermail 2.1.8 : Mon Jul 09 2012 - 10:45:32 MDT