diff -r 000000000000 -r 0c6405ab2ff4 taylor/01.ncl --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/taylor/01.ncl Mon Jan 26 22:08:20 2009 -0500 @@ -0,0 +1,57 @@ +;********************************** +; taylor_3.ncl +;********************************** +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" +load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" +load "./taylor_diagram.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 = (/ "Case A", "Case B" /) + nCase = dimsizes(case ) ; # of Cases [Cases] + +; variables compared + var = (/ "SLP","Tsfc","Prc","Prc 30S-30N","LW","SW", "U300", "Guess" /) + nVar = dimsizes(var) ; # of Variables + +; "Case A" + CA_ratio = (/1.230, 0.988, 1.092, 1.172, 1.064, 0.966, 1.079, 0.781/) + CA_cc = (/0.958, 0.973, 0.740, 0.743, 0.922, 0.982, 0.952, 0.433/) + +; "Case B" + CB_ratio = (/1.129, 0.996, 1.016, 1.134, 1.023, 0.962, 1.048, 0.852 /) + CB_cc = (/0.963, 0.975, 0.801, 0.814, 0.946, 0.984, 0.968, 0.647 /) + +; arrays to be passed to taylor plot + ratio = new ((/nCase, nVar/),typeof(CA_cc) ) + cc = new ((/nCase, nVar/),typeof(CA_cc) ) + + ratio(0,:) = CA_ratio + ratio(1,:) = CB_ratio + + cc(0,:) = CA_cc + cc(1,:) = CB_cc + +;********************************** +; create plot +;********************************** + + res = True ; default taylor diagram + + res@Markers = (/16, 16/) ; make all solid fill + res@Colors = (/"red", "blue" /) + res@varLabels = var + res@caseLabels = case + + wks = gsn_open_wks("ps","taylor") + plot = taylor_diagram(wks,ratio,cc,res) + + frame(wks) + +end