Date: Wed Sep 18 2013 - 14:58:32 MDT
Dear NCL
Attached is an image I made from the NCL Taylor plot examples.
Here is the part of the taylor_plot.ncl script that deals with the RMS circles.
My only problem is that I cannot get the RMS lines to draw correctly on the plot. How do I prevent them from extending past the outer-correlation arch on the plot?
Thank you.
Sincerely,
Erik
; ----------------------------------------------------------------
; Part 3:
; Concentric about 1.0 on XB axis
; I think this is correct. Still test mode.
; ----------------------------------------------------------------
if (rOpts .and. isatt(rOpts,"centerDiffRMS") \
.and. rOpts@centerDiffRMS) then
respl = True ; polyline mods desired
respl@xyLineThicknessF = 1.0 ; line thickness
respl@xyLineDashPattern = 2 ; short dash lines
respl@gsLineColor = "Black" ; line color
if (isatt(rOpts,"centerDiffRMS_color")) then
respl@gsLineColor = "LightGray"
end if
dx = 0.25
ncon = 8 ; 0.75, 0.50, 0.25, 0.0
npts = 100 ; arbitrary
ang = fspan(180,360,npts)*rad
dum9 = new(ncon,graphic)
;;;; IF YOU DESIRE RMS LABELS
rms_labels = (/0.0,0.25,0.50,0.75,1.0,1.25,1.50,1.75,2.0/)
do n=1,ncon
respl@gsLineLabelString = rms_labels(n-1) ;IF YOU DESIRE RMS LABELS
respl@gsLineLabelFontColor = "LightGray"
rr = n*dx ; radius from 1.0 [OBS] abscissa
xx = 1. + rr*cos(ang)
yy = fabs( rr*sin(ang) )
if (n.le.2) then
dum9(n-1) = gsn_add_polyline(wks,taylor,xx,yy,respl)
end if
if (n.eq.3) then
n3 = floattointeger( 0.61*npts )
dum9(n-1) = gsn_add_polyline(wks,taylor,xx(0:n3),yy(0:n3),respl)
dum9(n-1) = gsn_add_polyline(wks,taylor,xx,yy,respl)
end if
if (n.eq.4) then
n4 = floattointeger( 0.61*npts )
dum9(n-1) = gsn_add_polyline(wks,taylor,xx(0:n4),yy(0:n4),respl)
dum9(n-1) = gsn_add_polyline(wks,taylor,xx,yy,respl)
end if
if (n.eq.5) then
n5 = floattointeger( 0.22*npts )
dum9(n-1) = gsn_add_polyline(wks,taylor,xx(0:n5),yy(0:n5),respl)
dum9(n-1) = gsn_add_polyline(wks,taylor,xx,yy,respl)
end if
if (n.eq.6) then
n6 = floattointeger( 0.61*npts )
dum9(n-1) = gsn_add_polyline(wks,taylor,xx(0:n6),yy(0:n6),respl)
dum9(n-1) = gsn_add_polyline(wks,taylor,xx,yy,respl)
end if
if (n.eq.7) then
n7 = floattointeger( 0.61*npts )
dum9(n-1) = gsn_add_polyline(wks,taylor,xx(0:n7),yy(0:n7),respl)
dum9(n-1) = gsn_add_polyline(wks,taylor,xx,yy,respl)
end if
if (n.le.8) then
n8 = floattointeger( 0.71*npts )
dum9(n-1) = gsn_add_polyline(wks,taylor,xx(0:n8),yy(0:n8),respl)
dum9(n-1) = gsn_add_polyline(wks,taylor,xx,yy,respl)
end if
end do
delete(ang)
delete(xx)
delete(yy)
taylor@$unique_string("dum")$ = dum9
end if
_______________________________________________________
Erik U. Noble
Tel: +1 (212) 678-5512
Fax: +1 (212) 678-5552
E-mail: erik.noble@nasa.gov
WEB: http://www.giss.nasa.gov/staff/enoble.html
GISS: http://www.giss.nasa.gov
_______________________________________________________
* Open in Google Docs Viewer
* Open link in new tab
* Open link in new window
* Open link in new incognito window
*
* Download file
* Copy link address
* Edit PDF File on PDFescape.com
* Open in Google Docs Viewer
* Open link in new tab
* Open link in new window
* Open link in new incognito window
*
* Download file
* Copy link address
* Edit PDF File on PDFescape.com
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
