questions regarding Taylor diagram scripts

From: Erik N <nobleeu_at_nyahnyahspammersnyahnyah>
Date: Wed Jun 06 2012 - 13:00:06 MDT

Hi.

  How can I have RMS values appear on the RMS cycles?

  How can I turn on/off the numbers that appear with each point plotted on
the diagram?

  Lastly, how can I make some of the point labels into strings ( e.g...one
point says NNRP2 whereas another point nearby says MERRA?)

Can all of this be solved with makeing changes in the required
taylor_diagram.ncl script or can thye be solved using the interface script
(below)?

Below is the script I used to make the following picture.
Basically, I want to keep the red numbers as they are (or make them
disappear), but put labels on the blue green and orange points instead of
numbers.

Any tips? Thank you ahead of time.
Sincerely,
Erik

Erik Noble | Research Assistant | NASA Goddard Institute for Space Studies
Direct: 1 (212) 678-5512 | www.giss.nasa.gov | 2880 Broadway, New York,
NY-10025

;**********************************
; taylor_1_Erik.ncl
;**********************************
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "./taylor_diagram.ncl"
;**********************************
begin
    ;**********************************
    data =
readAsciiTable("../results_all_vs_trmm_7-18_ratio.text",5,"float",0)
    ; Only dataset. One text file. Split the data up into
    model=data(:,0)
            ; Model points
    cc_1=data(0:63,1)
    ratio_1=data(0:63,3)
            ; Precipitation Observations
    cc_2=data(64:66,1)
    ratio_2=data(64:66,3)
            ; Reanalysis Data
    cc_3=data(67:68,1)
    ratio_3=data(67:68,3)
            ; Benchmark Model
    cc_4=data(69,1)
    ratio_4=data(69,3)

    ;**********************************
    ; Put the ratios and pattern correlations into
    ; arrays for plotting
    ;**********************************
    nDataSets = 4 ; number of datasets
    npts = dimsizes(data(:,0))
    ratio = new ((/nDataSets, npts/),typeof(ratio_1) )
    cc = new ((/nDataSets, npts/),typeof(cc_1) )

    ratio(0,0:63) = ratio_1
    cc(0,0:63) = cc_1

    ratio(1,64:66) = ratio_2
        cc(1,64:66) = cc_2

    ratio(2,67:68) = ratio_3
        cc(2,67:68) = cc_3

    ratio(3,69) = ratio_4
        cc(3,69) = cc_4
    ;**********************************
    ; create plot
    ;**********************************
    wks = gsn_open_wks("pdf","taylor_test")

    res = True ; diagram mods desired
    res@tiMainString = "Example" ; title
    res@Colors = (/"red","blue","green","orange"/) ;
marker colors
    res@Markers = (/2,7,6,15/) ; marker styles
    res@markerTxYOffset = 0.04 ; offset btwn marker &
label
    ;res@gsMarkerSizeF = 0.008 ; marker size
    ;res@txFontHeightF = 0.039 ; text size

    res@stnRad = (/ 0.5, 1.5 /) ; additional standard radii
    res@ccRays = (/ 0.6, 0.9 /) ; correllation rays
    res@ccRays_color = "LightGray" ; default is "black"

    res@centerDiffRMS = True ; RMS 'circles'
    res@centerDiffRMS_color = "LightGray" ; default is "black"/)

    plot = taylor_diagram(wks, ratio, cc, res)
end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

Received on Wed Jun 6 13:00:34 2012

This archive was generated by hypermail 2.1.8 : Wed Jun 06 2012 - 15:17:44 MDT