Modifications to scoring and graphics production for the final version of code for the C-LAMP paper in GCB.
1 ;**********************************
3 ;**********************************
4 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
5 load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
6 load "./taylor_diagram.ncl"
7 ;**********************************
9 ;**********************************
10 ; Assume the following have already been computed:
11 ; _ratio are the ratio: Case_Variance/Reference_Variance
12 ; _cc are the cross correlation coef of Case to Reference
13 ;**********************************
16 case = (/ "Case A", "Case B" /)
17 nCase = dimsizes(case ) ; # of Cases [Cases]
20 var = (/ "SLP","Tsfc","Prc","Prc 30S-30N","LW","SW", "U300", "Guess" /)
21 nVar = dimsizes(var) ; # of Variables
24 CA_ratio = (/1.230, 0.988, 1.092, 1.172, 1.064, 0.966, 1.079, 0.781/)
25 CA_cc = (/0.958, 0.973, 0.740, 0.743, 0.922, 0.982, 0.952, 0.433/)
28 CB_ratio = (/1.129, 0.996, 1.016, 1.134, 1.023, 0.962, 1.048, 0.852 /)
29 CB_cc = (/0.963, 0.975, 0.801, 0.814, 0.946, 0.984, 0.968, 0.647 /)
31 ; arrays to be passed to taylor plot
32 ratio = new ((/nCase, nVar/),typeof(CA_cc) )
33 cc = new ((/nCase, nVar/),typeof(CA_cc) )
41 ;**********************************
43 ;**********************************
45 res = True ; default taylor diagram
47 res@Markers = (/16, 16/) ; make all solid fill
48 res@Colors = (/"red", "blue" /)
52 wks = gsn_open_wks("ps","taylor")
53 plot = taylor_diagram(wks,ratio,cc,res)