<div dir="ltr"><div><div><div><div><div>Hi Guilherme,<br><br></div>In function get_areas_of_interest(), line 28:<br><br>return(features(ind(.not.ismissing(features))))<br><br></div>for whatever reason, "features" contains all missing values, and .not.ismissing(features) returns all False. The docs for ind() state that if there are no values that evaluate to true, it returns a missing value, which for integer types is -2147483647, and that's where the index out of range comes from :-/<br><br></div>So, I don't know why features is coming up all missing, but perhaps that gives and idea where to look.<br><br></div>Hope that helps...<br></div>Rick<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Apr 20, 2017 at 6:27 AM, Guilherme Martins <span dir="ltr"><<a href="mailto:jgmsantos@gmail.com" target="_blank">jgmsantos@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear Mary,<div><br></div><div>Thank you for all informations.</div><div><br></div><div>I'm finishing my script using your informations (shapefiles_21.ncl) but I found a problem.</div><div><br></div><div>When I run my script, I got that message:</div><div><br></div><div><div><div><div>fatal:Subscript out of range, error in subscript #0</div><div>fatal:An error occurred reading features</div><div>fatal:["Execute.c":8575]:<wbr>Execute: Error occurred at or near line 28 in file teste.ncl</div><div><br></div><div>fatal:["Execute.c":8575]:<wbr>Execute: Error occurred at or near line 132 in file teste.ncl</div></div></div></div><div><br></div><div><p style="margin:0px">I am using the "sigla" variable that has 5 informations (S, SE, N, NE, CO), and I can only pass 4 parameters in the "get_areas_of_interest" function. If I try to pass 5 parameters <span style="line-height:1.5">(S, SE, N, NE, CO)</span><span style="line-height:1.5">, the script returns the above error. T</span><span style="line-height:1.5">he goal is to use the 5 parameters in the function.</span></p></div><div><br></div><div>The script and data can download at:</div><div><br></div><div><a href="http://ftp-supercomputacao.inpe.br/public/jose.martins/teste.tar.gz" target="_blank">http://ftp-supercomputacao.<wbr>inpe.br/public/jose.martins/<wbr>teste.tar.gz</a><br></div><div><br></div><div>The shapefile informations (*.png) are in attachment.</div><span class="HOEnZb"><font color="#888888"><div><br></div><div><span style="line-height:1.5">Guilherme Martins.</span><br></div></font></span><div><div class="h5"><div><br></div><div><br></div><br><div class="gmail_quote"><div dir="ltr">Em seg, 17 de abr de 2017 às 18:02, Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>> escreveu:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-size:small">Guilherme,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Please see examples shapefiles_21.ncl at:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/shapefiles.shtml#ex21" target="_blank">http://www.ncl.ucar.edu/<wbr>Applications/shapefiles.shtml#<wbr>ex21</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">It explains why you see these "blocky" contours in NCL.</div><div class="gmail_default"><br></div><div class="gmail_default">Simply put, if you have four points that represent a box at the edge of a domain, and three of them have data one them and one is missing, then NCL will only draw contours in half of the box, because it won't extrapolate to that missing point.</div><div class="gmail_default"><br></div><div class="gmail_default">The real question is: do you simply want the graphics not to have gaps, or do you actually want there to be data in locations that fall just outside the geographical boundary? You can do things to force values just outside the geographical boundary, but this depends on how you are masking in the first place.</div><div class="gmail_default"><br></div><div class="gmail_default">If you are masking based on a shapefile, then see example shapefiles_18.ncl on the same page:</div><div class="gmail_default"><br></div><div class="gmail_default"><a href="http://www.ncl.ucar.edu/Applications/shapefiles.shtml#ex18" target="_blank">http://www.ncl.ucar.edu/<wbr>Applications/shapefiles.shtml#<wbr>ex18</a><br></div><div class="gmail_default"><br></div><div class="gmail_default">Note, however, that this mask is now returning values outside the range of interest.</div><div class="gmail_default"><br></div><div class="gmail_default">If you are simply trying to mask your data *graphically* by some geographical area, then instead of trying to use a mask function, you can draw the full data, and then draw polygons filled in with white in the areas you don't want to see graphically. That's what example 21 is showing.</div><div class="gmail_default"><br></div><div class="gmail_default">--Mary</div><div class="gmail_default"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote"></div></div><div class="gmail_extra"><div class="gmail_quote">On Mon, Apr 17, 2017 at 1:00 PM, Guilherme Martins <span dir="ltr"><<a href="mailto:jgmsantos@gmail.com" target="_blank">jgmsantos@gmail.com</a>></span> wrote:<br></div></div><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Dear users,<div><br></div><div>After searching extensively on the NCL website I was not able to find a solution to my problem. I only want values within the contour of the continent, but some edges appear in the figure that are stemming from the mask that I applied and I do not know how to remove this. The figure is attached for a better view of my problem.<br></div><div><br></div><div>Best regards,</div><div><br></div><div>Guilherme.</div></div><span class="m_-699587347489471141m_1770811509721896995HOEnZb"><font color="#888888"><div dir="ltr">-- <br></div><div data-smartmail="gmail_signature"><div dir="ltr"><div>------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>------</div><div><div>Instituto Nacional de Pesquisas Espaciais (INPE)<br></div><div>Centro de Previsão de Tempo e Estudos Climáticos (CPTEC)</div><div>Divisão de Satélites e Sistemas Ambientais (DSA)</div><div>Programa de Monitoramento de Queimadas</div><div>Telefone (INPE/CP): <a href="tel:+55%2012%203186-9205" value="+551231869205" target="_blank">+55 12 3186-9205</a> || Celular (TIM): <a href="tel:+55%2012%2098111-4292" value="+5512981114292" target="_blank">+55 12 98111-4292</a></div><div>E-mail: <a href="mailto:guilherme.martins@inpe.br" target="_blank">guilherme.martins@inpe.br</a> || <a href="mailto:jgmsantos@gmail.com" target="_blank">jgmsantos@gmail.com</a></div><div>Skype: guilherme.martins.</div><div>Homepage: <a href="https://sites.google.com/site/jgmsantos" target="_blank">https://sites.google.com/site/<wbr>jgmsantos</a></div><div>Currículo Lattes: <a href="http://lattes.cnpq.br/5997657584785803" target="_blank">http://lattes.cnpq.br/<wbr>5997657584785803</a></div></div><div><span style="font-family:"helvetica neue",helvetica,arial,sans-serif">------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>------</span><br></div></div></div>
</font></span><br></blockquote></div></div><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div></div></div></div><div class="HOEnZb"><div class="h5"><div dir="ltr">-- <br></div><div data-smartmail="gmail_signature"><div dir="ltr"><div>------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>------</div><div><div>Instituto Nacional de Pesquisas Espaciais (INPE)<br></div><div>Centro de Previsão de Tempo e Estudos Climáticos (CPTEC)</div><div>Divisão de Satélites e Sistemas Ambientais (DSA)</div><div>Programa de Monitoramento de Queimadas</div><div>Telefone (INPE/CP): <a href="tel:+55%2012%203186-9205" value="+551231869205" target="_blank">+55 12 3186-9205</a> || Celular (TIM): <a href="tel:+55%2012%2098111-4292" value="+5512981114292" target="_blank">+55 12 98111-4292</a></div><div>E-mail: <a href="mailto:guilherme.martins@inpe.br" target="_blank">guilherme.martins@inpe.br</a> || <a href="mailto:jgmsantos@gmail.com" target="_blank">jgmsantos@gmail.com</a></div><div>Skype: guilherme.martins.</div><div>Homepage: <a href="https://sites.google.com/site/jgmsantos" target="_blank">https://sites.google.com/site/<wbr>jgmsantos</a></div><div>Currículo Lattes: <a href="http://lattes.cnpq.br/5997657584785803" target="_blank">http://lattes.cnpq.br/<wbr>5997657584785803</a></div></div><div><span style="font-family:"helvetica neue",helvetica,arial,sans-serif">------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>------</span><br></div></div></div>
</div></div><br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>