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