<div><div>Dear NCL developers,</div><div><br></div><div>I use the script "taylor_diagram.ncl" to draw Taylor diagram. I extend the outer curve from the original 1.65 to 1.75 to fit my data, but it seems that the parameters w.r.t. RMS semi-circles need also be changed accordingly (as is shown below, the semi-circles are incomplete).</div><div>The issue is, I have no idea how to change it.</div><div><br></div><div><img src="cid:2E05F4FE@0FC9496C.9F046F55"></div><div><br></div><div>Below is a snippet of the script, which I guess is related to the RMS semi-circles</div><div>Any hints about how to change the parameters would be sincerely appreciated!</div><div><br></div><div>Best,</div><div>Qi</div><div><br></div><div>dx &nbsp; = 0.25</div><div>&nbsp; &nbsp; &nbsp; ncon = 4 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; 0.75, 0.50, 0.25, 0.0</div><div>&nbsp; &nbsp; &nbsp; npts = 100 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; ; arbitrary</div><div>&nbsp; &nbsp; &nbsp; ang &nbsp;= fspan(180,360,npts)*rad</div><div><br></div><div>&nbsp; &nbsp; &nbsp; dum9 = new(ncon,graphic)</div><div><br></div><div>&nbsp; &nbsp; &nbsp; do n=1,ncon&nbsp;</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;rr &nbsp;= n*dx &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;; radius from 1.0 [OBS] abscissa</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;xx &nbsp;= 1. + rr*cos(ang)</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;yy &nbsp;= fabs( rr*sin(ang) )</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;if (n.le.2) then</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;dum9(n-1) = gsn_add_polyline(wks,taylor,xx,yy,respl)</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;end if</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;if (n.eq.3) then</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;n3 = floattointeger( 0.77*npts )&nbsp;</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;dum9(n-1) = gsn_add_polyline(wks,taylor,xx(0:n3),yy(0:n3),respl)</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;end if</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;if (n.eq.4) then</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;n4 = floattointeger( 0.61*npts )&nbsp;</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;dum9(n-1) = gsn_add_polyline(wks,taylor,xx(0:n4),yy(0:n4),respl)</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;end if</div><div>&nbsp; &nbsp; &nbsp; end do</div></div>