Re: expansion half circle of Taylor Diagram,

From: Wilhelm Haygen <w.haygen_at_nyahnyahspammersnyahnyah>
Date: Thu Jul 05 2012 - 00:01:16 MDT

Dear Dennis,
Thanks for the explaination. Yeah i used that one. I look forward to the
new taylor diagram with half circle expansion. Thanks a lot.

On Thursday, July 5, 2012, Dennis Shea wrote:

> 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
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Jul 5 00:01:31 2012

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