<div dir="ltr">Forgot to Cc ncl-talk...I apologize!<br><br><div><div class="gmail_quote">---------- Forwarded message ----------<br>From: <b class="gmail_sendername">Carlos J. Valle-Diaz</b> <span dir="ltr">&lt;<a href="mailto:cj.vallediaz@gmail.com" target="_blank">cj.vallediaz@gmail.com</a>&gt;</span><br>Date: Wed, Oct 21, 2015 at 8:27 PM<br>Subject: Re: [ncl-talk] Help with creating a script for precipitation stored in shapefiles<br>To: Rick Brownrigg &lt;<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>&gt;<br><br><div dir="ltr"><div><div><div><div>Thank you very much Rick,<br><br></div>I&#39;ve been playing around with the script to understand the functions and resources used to make the plot and learn NCL from it. I managed to focus it on my area of interest (Puerto Rico). At the moment, I&#39;ve been trying to put a legend and text. I&#39;ve seen an example at the website &quot;polyg_8.ncl&quot;, but I have a question.<br><br></div>I know that the gsn_polymarker_ndc(wks, x, y, res) function lets me put the legend markers and the gsn_text_ndc(wks, text, x, y, res) function lets me put the text. I&#39;m using &quot;gsres&quot; as a resource for the polymarkers and &quot;txres&quot; as a resource for the text to go along with it. My question lies in the numeric values (x and y). In the example script I mentioned above, the code is set for an array set defined. In my case I would have to define it for the Globalvalues (precipitation) in my shapefile. In summary, how do I define x and y to make a legend at an even interval of my GlobalValues and to play with the location of both the markers and text? (See script attached)<br><br></div><div>Also, I&#39;ve been trying to include black circles in the loop that generates the colored polymarkers in the plot to make it more visible, but haven&#39;t been succesful. I&#39;ve tried the following sequence, for example, and marks an error.  What am I doing wrong?<br></div><div><br>   do i=0, dimsizes(f-&gt;Globvalue) - 1<br><br>      indx = doubletoint( (f-&gt;Globvalue(i)-globMin) / globRange * numColors) + 1  ; skip first color in table<br>      if (indx.ge.numColors) then<br>          indx = numColors - 1<br>      end if<br><br>       gsres@gsMarkerIndex = 4<br>       gsres@gsMarkerThicknessF = 2<br>       gsres@gsMarkerColor = &quot;Black&quot;<br>       gsn_polymarker(wks, plot, f-&gt;Lon(i), f-&gt;Lat(i), gsres)<br>       gsres@gsMarkerIndex = 1<br>       gsres@gsMarkerSizeF = 0.03<br>       gsres@gsMarkerColor = cmap(indx, :)<br>       gsn_polymarker(wks, plot, f-&gt;Lon(i), f-&gt;Lat(i), gsres)<br><br>   end do<br><br></div>Thanks again,<br><br></div>Carlos<br></div><br></div></div></div>