Re: expansion half circle of Taylor Diagram,

To my knowledge, no one has done a "half circle."
However, the original code has been modified to include
indications of bias and the negative correlations
are added at the bottom.

http://www.ncl.ucar.edu/Applications/taylor.shtml
Example 8

Download: taylor_diagram_cam.ncl

On 7/4/12 4:25 PM, antonio avio wrote:
> 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
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jul 4 17:07:45 2012