<div dir="ltr"><div><div><div><div>Dear All,<br><br></div>I am trying to understand how to do the wavelet analysis and started using an example code and the data available online. But I got an error message as &quot;segmentation fault (core dumped)&quot; when I ran the program. I have given the code that I used as below for you to review. Also I have attached the data I used in the email. <br><br></div>So, could you please suggest me how I solve this problem? <br><br></div>Thanks in advance.<br><br>THE CODE is here and the data is attached in the email,<br><br> load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;<br> load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;<br> load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl&quot;<br> load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl&quot;<br> load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl&quot;<br><br>;--------------------------------------------------------------------<br>; CREATE SAME WAVELET FIGURE AS TORRENCE &amp; COMPO using seasonal Nino3<br>; from 1871-1997<br>begin<br>ninoseas = asciiread(&quot;C:\diptic\sst\sst_nino3.dat&quot;,-1,&quot;float&quot;)<br>  ninoseas!0 = &quot;time&quot;<br>  ntime = dimsizes(ninoseas)<br>  timeo = fspan(1871.25,1996.,ntime)<br>  print(ntime)<br>  ninoseas&amp;time = timeo<br>  ninomam = dim_avg_Wrap(ninoseas(0::4))<br>  ninojja = dim_avg_Wrap(ninoseas(1::4))<br>  ninoson = dim_avg_Wrap(ninoseas(2::4))<br>  ninodjf = dim_avg_Wrap(ninoseas(3::4))<br>  ninoseas(0::4) = ninoseas(0::4) - ninomam<br>  ninoseas(1::4) = ninoseas(1::4) - ninojja<br>  ninoseas(2::4) = ninoseas(2::4) - ninoson<br>  ninoseas(3::4) = ninoseas(3::4) - ninodjf<br>  time  = timeo<br>  N     = dimsizes(time)<br>;_FillValue=1e+36<br>;printMinMax(ninoseas,False)<br>;return<br>;************************************<br>; compute wavelet<br>;************************************<br>  mother  = 0<br>  param   = 6.0<br>  dt      = 0.25    ;timesteps in units of years<br>  s0      = dt<br>  dj      = 0.25<br>  jtot    = 1+floattointeger(((log10(N*dt/s0))/dj)/log10(2.))<br>  npad    = N<br>  nadof   = 0.0<br>  noise   = 1<br>  siglvl  = 0.05<br>  isigtest= 0<br>;print(nadof)<br>w =wavelet(ninoseas,mother,dt,param,s0,dj,jtot,npad,noise,isigtest,siglvl,nadof)<br>;printVarSummary(w)<br>return<br>end<br></div>segmentation fault (core dumped)<br><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Apr 18, 2017 at 3:29 AM,  <span dir="ltr">&lt;<a href="mailto:ncl-talk-request@ucar.edu" target="_blank">ncl-talk-request@ucar.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send ncl-talk mailing list submissions to<br>
        <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<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>
or, via email, send a message with subject or body &#39;help&#39; to<br>
        <a href="mailto:ncl-talk-request@ucar.edu">ncl-talk-request@ucar.edu</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:ncl-talk-owner@ucar.edu">ncl-talk-owner@ucar.edu</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than &quot;Re: Contents of ncl-talk digest...&quot;<br>
<br>
<br>
Today&#39;s Topics:<br>
<br>
   1. Re: Problem with gsn_add_text (Stavros Dafis)<br>
   2. Problem with contour plot (Guilherme Martins)<br>
   3. Re: Problem with gsn_add_text (Rick Brownrigg)<br>
   4. Re: Problem with contour plot (Mary Haley)<br>
   5. how to draw this kind of plot (wen)<br>
<br>
<br>
------------------------------<wbr>------------------------------<wbr>----------<br>
<br>
Message: 1<br>
Date: Mon, 17 Apr 2017 21:48:31 +0300<br>
From: Stavros Dafis &lt;<a href="mailto:sdafis@cc.uoi.gr">sdafis@cc.uoi.gr</a>&gt;<br>
Subject: Re: [ncl-talk] Problem with gsn_add_text<br>
To: Rick Brownrigg &lt;<a href="mailto:brownrig@ucar.edu">brownrig@ucar.edu</a>&gt;<br>
Cc: &quot;<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>&quot; &lt;<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>&gt;<br>
Message-ID: &lt;<a href="mailto:1492454911.58f50dffc9203@webmail.uoi.gr">1492454911.58f50dffc9203@<wbr>webmail.uoi.gr</a>&gt;<br>
Content-Type: text/plain; charset=ISO-8859-7<br>
<br>
Dear Rick, thank you for your suggestion, I have tried to assign both the<br>
FillValue and the missing value (@missing_value) but nothing seems to work.<br>
<br>
I forgot to mention that I am using NCL 6.2.0, 6.3.0 and 6.4.0, on different<br>
machines.<br>
<br>
<br>
--<br>
Stavros NTAFIS (DAFIS)<br>
------------------------------<wbr>-----------------<br>
Physicist - Meteorologist, M.Sc.<br>
Research Associate, National Observatory of Athens<br>
(+33)9 81 94 22 12, Mobile: +30 697 04 20 242<br>
---------------<br>
Weather charts:<br>
<a href="http://www.meteo.gr" rel="noreferrer" target="_blank">http://www.meteo.gr</a><br>
<a href="http://www.meteo.gr/meteomaps/" rel="noreferrer" target="_blank">http://www.meteo.gr/meteomaps/</a><br>
<br>
<br>
Quoting Rick Brownrigg &lt;<a href="mailto:brownrig@ucar.edu">brownrig@ucar.edu</a>&gt;:<br>
<br>
&gt; Hi Stavros,<br>
&gt;<br>
&gt; I don&#39;t know for certain, but I suspect you need to assign the _FillValue<br>
&gt; to the string array &quot;points&quot; that you constructed from &quot;t2&quot; (I don&#39;t see<br>
&gt; how the attributes from t2 would carry to points via the sprintf() ).<br>
&gt;<br>
&gt; Hope that helps...<br>
&gt; Rick<br>
&gt;<br>
&gt; On Sun, Apr 16, 2017 at 12:56 PM, Stavros Dafis &lt;<a href="mailto:sdafis@cc.uoi.gr">sdafis@cc.uoi.gr</a>&gt; wrote:<br>
&gt;<br>
&gt; &gt; Hello everyone, I am facing a problem with the function gsn_add_text. In<br>
&gt; &gt; the<br>
&gt; &gt; description of the function, it is mentioned that no missing values will be<br>
&gt; &gt; plotted:<br>
&gt; &gt;<br>
&gt; &gt; &amp;#8220;This function attaches a list of text strings to the given plot.<br>
&gt; &gt; The text<br>
&gt; &gt; strings will be drawn only when the plot is drawn. If a missing value is<br>
&gt; &gt; encountered in the text strings or x/y values, then no text string will be<br>
&gt; &gt; drawn at this pair.&amp;#8221;<br>
&gt; &gt;<br>
&gt; &gt; My script reads a netCDF file with 3 variables t2(lat,lon), lat(lat),<br>
&gt; &gt; lon(lon),<br>
&gt; &gt; so in order to attach values on the plot using this function, I have to<br>
&gt; &gt; make<br>
&gt; &gt; new 2-D variables for lat/lon and define the string to be plotted on the<br>
&gt; &gt; map.<br>
&gt; &gt; Values are assigned just fine on the plot, over the land, but missing<br>
&gt; &gt; values<br>
&gt; &gt; are also printed over the sea. I tried to use landsea masking but my data<br>
&gt; &gt; have<br>
&gt; &gt; high horizontal resolution (5km). Any suggestions?<br>
&gt; &gt;<br>
&gt; &gt; I have the same problem with WRF data, with already 2-D lat/lon arrays.<br>
&gt; &gt;<br>
&gt; &gt;  This is the part of the code which prints the values:<br>
&gt; &gt; .<br>
&gt; &gt; .<br>
&gt; &gt; .<br>
&gt; &gt; contour = gsn_csm_contour_map(wks,t2,<wbr>res)<br>
&gt; &gt;<br>
&gt; &gt; t2@missing_value = -1<br>
&gt; &gt; t2@_FillValue = -1<br>
&gt; &gt; (default missing values = -999)<br>
&gt; &gt;<br>
&gt; &gt; lat2 = new(dimsizes(t2),typeof(lat),-<wbr>1)<br>
&gt; &gt; lon2 = new(dimsizes(t2),typeof(lon),-<wbr>1)<br>
&gt; &gt;<br>
&gt; &gt; nb = dimsizes(t2)<br>
&gt; &gt;<br>
&gt; &gt; do nl=0,nb(0)-1,1<br>
&gt; &gt; do ml=0,nb(1)-1,1<br>
&gt; &gt;<br>
&gt; &gt; lat2(nl,ml) = lat(nl)<br>
&gt; &gt; lon2(nl,ml) = lon(ml)<br>
&gt; &gt;<br>
&gt; &gt; end do<br>
&gt; &gt; end do<br>
&gt; &gt;<br>
&gt; &gt; points = sprintf(&quot;%3.0f&quot;, t2)<br>
&gt; &gt; ;points = &quot;&quot;+t2(:,:)+&quot;&quot;<br>
&gt; &gt;<br>
&gt; &gt; values = gsn_add_text(wks,contour,<wbr>points,lon2(:,:),lat2(:,:),<wbr>True)<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; --<br>
&gt; &gt; Stavros NTAFIS (DAFIS)<br>
&gt; &gt; ------------------------------<wbr>-----------------<br>
&gt; &gt; Physicist - Meteorologist, M.Sc.<br>
&gt; &gt; Research Associate, National Observatory of Athens<br>
&gt; &gt; (+33)9 81 94 22 12, Mobile: +30 697 04 20 242<br>
&gt; &gt; ---------------<br>
&gt; &gt; Weather charts:<br>
&gt; &gt; <a href="http://www.meteo.gr" rel="noreferrer" target="_blank">http://www.meteo.gr</a><br>
&gt; &gt; <a href="http://www.meteo.gr/meteomaps/" rel="noreferrer" target="_blank">http://www.meteo.gr/meteomaps/</a><br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; ______________________________<wbr>_________________<br>
&gt; &gt; ncl-talk mailing list<br>
&gt; &gt; <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
&gt; &gt; List instructions, subscriber options, unsubscribe:<br>
&gt; &gt; <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>
&gt; &gt;<br>
&gt;<br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Mon, 17 Apr 2017 19:00:02 +0000<br>
From: Guilherme Martins &lt;<a href="mailto:jgmsantos@gmail.com">jgmsantos@gmail.com</a>&gt;<br>
Subject: [ncl-talk] Problem with contour plot<br>
To: NCL &lt;<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>&gt;<br>
Message-ID:<br>
        &lt;CAH+t_<wbr>hMaRUDHM0wMkLC9cz2pFmMYGXgGC+<wbr>X+bcxDLo=+<a href="mailto:8fgWQg@mail.gmail.com">8fgWQg@mail.gmail.<wbr>com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;utf-8&quot;<br>
<br>
Dear users,<br>
<br>
After searching extensively on the NCL website I was not able to find a<br>
solution to my problem. I only want values within the contour of the<br>
continent, but some edges appear in the figure that are stemming from the<br>
mask that I applied and I do not know how to remove this. The figure is<br>
attached for a better view of my problem.<br>
<br>
Best regards,<br>
<br>
Guilherme.<br>
--<br>
------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>------<br>
Instituto Nacional de Pesquisas Espaciais (INPE)<br>
Centro de Previs?o de Tempo e Estudos Clim?ticos (CPTEC)<br>
Divis?o de Sat?lites e Sistemas Ambientais (DSA)<br>
Programa de Monitoramento de Queimadas<br>
Telefone (INPE/CP): <a href="tel:%2B55%2012%203186-9205" value="+551231869205">+55 12 3186-9205</a> || Celular (TIM): <a href="tel:%2B55%2012%2098111-4292" value="+5512981114292">+55 12 98111-4292</a><br>
E-mail: <a href="mailto:guilherme.martins@inpe.br">guilherme.martins@inpe.br</a> || <a href="mailto:jgmsantos@gmail.com">jgmsantos@gmail.com</a><br>
Skype: guilherme.martins.<br>
Homepage: <a href="https://sites.google.com/site/jgmsantos" rel="noreferrer" target="_blank">https://sites.google.com/site/<wbr>jgmsantos</a><br>
Curr?culo Lattes: <a href="http://lattes.cnpq.br/5997657584785803" rel="noreferrer" target="_blank">http://lattes.cnpq.br/<wbr>5997657584785803</a><br>
------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>------<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170417/c3215b0e/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20170417/c3215b0e/<wbr>attachment-0001.html</a><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: figure.png<br>
Type: image/png<br>
Size: 58495 bytes<br>
Desc: not available<br>
Url : <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170417/c3215b0e/attachment-0001.png" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20170417/c3215b0e/<wbr>attachment-0001.png</a><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Mon, 17 Apr 2017 14:32:43 -0600<br>
From: Rick Brownrigg &lt;<a href="mailto:brownrig@ucar.edu">brownrig@ucar.edu</a>&gt;<br>
Subject: Re: [ncl-talk] Problem with gsn_add_text<br>
To: Stavros Dafis &lt;<a href="mailto:sdafis@cc.uoi.gr">sdafis@cc.uoi.gr</a>&gt;<br>
Cc: &quot;<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>&quot; &lt;<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>&gt;<br>
Message-ID:<br>
        &lt;CAGKRhbHX1v-1sJ3G=<a href="mailto:QrXGA%2B0wSJ2KKsBB_7GqbCnp9krh9fkkQ@mail.gmail.com">QrXGA+<wbr>0wSJ2KKsBB_7GqbCnp9krh9fkkQ@<wbr>mail.gmail.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;utf-8&quot;<br>
<br>
Hi,<br>
<br>
You used -1 as _FillValue for variable &quot;t2&quot;. Note it is an integer value --<br>
what type is the array t2 -- float or int?  In any case, you created the<br>
formatted string array &quot;points&quot; from t2 as &quot;%3.0f&quot;, so that any data values<br>
in t2 that were -1 are now &quot; -1&quot;, and thus it seems your _FillValue should<br>
be &quot; -1&quot;. Did you do this, and do the types between data and _FillValue<br>
agree everywhere?<br>
<br>
If that&#39;s not the issue, I think we may need to see script and dataset to<br>
figure out what the issue is.<br>
<br>
HTH...<br>
Rick<br>
<br>
On Mon, Apr 17, 2017 at 12:48 PM, Stavros Dafis &lt;<a href="mailto:sdafis@cc.uoi.gr">sdafis@cc.uoi.gr</a>&gt; wrote:<br>
<br>
&gt; Dear Rick, thank you for your suggestion, I have tried to assign both the<br>
&gt; FillValue and the missing value (@missing_value) but nothing seems to work.<br>
&gt;<br>
&gt; I forgot to mention that I am using NCL 6.2.0, 6.3.0 and 6.4.0, on<br>
&gt; different<br>
&gt; machines.<br>
&gt;<br>
&gt;<br>
&gt; --<br>
&gt; Stavros NTAFIS (DAFIS)<br>
&gt; ------------------------------<wbr>-----------------<br>
&gt; Physicist - Meteorologist, M.Sc.<br>
&gt; Research Associate, National Observatory of Athens<br>
&gt; (+33)9 81 94 22 12, Mobile: +30 697 04 20 242<br>
&gt; ---------------<br>
&gt; Weather charts:<br>
&gt; <a href="http://www.meteo.gr" rel="noreferrer" target="_blank">http://www.meteo.gr</a><br>
&gt; <a href="http://www.meteo.gr/meteomaps/" rel="noreferrer" target="_blank">http://www.meteo.gr/meteomaps/</a><br>
&gt;<br>
&gt;<br>
&gt; Quoting Rick Brownrigg &lt;<a href="mailto:brownrig@ucar.edu">brownrig@ucar.edu</a>&gt;:<br>
&gt;<br>
&gt; &gt; Hi Stavros,<br>
&gt; &gt;<br>
&gt; &gt; I don&#39;t know for certain, but I suspect you need to assign the _FillValue<br>
&gt; &gt; to the string array &quot;points&quot; that you constructed from &quot;t2&quot; (I don&#39;t see<br>
&gt; &gt; how the attributes from t2 would carry to points via the sprintf() ).<br>
&gt; &gt;<br>
&gt; &gt; Hope that helps...<br>
&gt; &gt; Rick<br>
&gt; &gt;<br>
&gt; &gt; On Sun, Apr 16, 2017 at 12:56 PM, Stavros Dafis &lt;<a href="mailto:sdafis@cc.uoi.gr">sdafis@cc.uoi.gr</a>&gt;<br>
&gt; wrote:<br>
&gt; &gt;<br>
&gt; &gt; &gt; Hello everyone, I am facing a problem with the function gsn_add_text.<br>
&gt; In<br>
&gt; &gt; &gt; the<br>
&gt; &gt; &gt; description of the function, it is mentioned that no missing values<br>
&gt; will be<br>
&gt; &gt; &gt; plotted:<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; &amp;#8220;This function attaches a list of text strings to the given plot.<br>
&gt; &gt; &gt; The text<br>
&gt; &gt; &gt; strings will be drawn only when the plot is drawn. If a missing value<br>
&gt; is<br>
&gt; &gt; &gt; encountered in the text strings or x/y values, then no text string<br>
&gt; will be<br>
&gt; &gt; &gt; drawn at this pair.&amp;#8221;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; My script reads a netCDF file with 3 variables t2(lat,lon), lat(lat),<br>
&gt; &gt; &gt; lon(lon),<br>
&gt; &gt; &gt; so in order to attach values on the plot using this function, I have to<br>
&gt; &gt; &gt; make<br>
&gt; &gt; &gt; new 2-D variables for lat/lon and define the string to be plotted on<br>
&gt; the<br>
&gt; &gt; &gt; map.<br>
&gt; &gt; &gt; Values are assigned just fine on the plot, over the land, but missing<br>
&gt; &gt; &gt; values<br>
&gt; &gt; &gt; are also printed over the sea. I tried to use landsea masking but my<br>
&gt; data<br>
&gt; &gt; &gt; have<br>
&gt; &gt; &gt; high horizontal resolution (5km). Any suggestions?<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; I have the same problem with WRF data, with already 2-D lat/lon arrays.<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;  This is the part of the code which prints the values:<br>
&gt; &gt; &gt; .<br>
&gt; &gt; &gt; .<br>
&gt; &gt; &gt; .<br>
&gt; &gt; &gt; contour = gsn_csm_contour_map(wks,t2,<wbr>res)<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; t2@missing_value = -1<br>
&gt; &gt; &gt; t2@_FillValue = -1<br>
&gt; &gt; &gt; (default missing values = -999)<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; lat2 = new(dimsizes(t2),typeof(lat),-<wbr>1)<br>
&gt; &gt; &gt; lon2 = new(dimsizes(t2),typeof(lon),-<wbr>1)<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; nb = dimsizes(t2)<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; do nl=0,nb(0)-1,1<br>
&gt; &gt; &gt; do ml=0,nb(1)-1,1<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; lat2(nl,ml) = lat(nl)<br>
&gt; &gt; &gt; lon2(nl,ml) = lon(ml)<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; end do<br>
&gt; &gt; &gt; end do<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; points = sprintf(&quot;%3.0f&quot;, t2)<br>
&gt; &gt; &gt; ;points = &quot;&quot;+t2(:,:)+&quot;&quot;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; values = gsn_add_text(wks,contour,<wbr>points,lon2(:,:),lat2(:,:),<wbr>True)<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; --<br>
&gt; &gt; &gt; Stavros NTAFIS (DAFIS)<br>
&gt; &gt; &gt; ------------------------------<wbr>-----------------<br>
&gt; &gt; &gt; Physicist - Meteorologist, M.Sc.<br>
&gt; &gt; &gt; Research Associate, National Observatory of Athens<br>
&gt; &gt; &gt; (+33)9 81 94 22 12, Mobile: +30 697 04 20 242<br>
&gt; &gt; &gt; ---------------<br>
&gt; &gt; &gt; Weather charts:<br>
&gt; &gt; &gt; <a href="http://www.meteo.gr" rel="noreferrer" target="_blank">http://www.meteo.gr</a><br>
&gt; &gt; &gt; <a href="http://www.meteo.gr/meteomaps/" rel="noreferrer" target="_blank">http://www.meteo.gr/meteomaps/</a><br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; ______________________________<wbr>_________________<br>
&gt; &gt; &gt; ncl-talk mailing list<br>
&gt; &gt; &gt; <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
&gt; &gt; &gt; List instructions, subscriber options, unsubscribe:<br>
&gt; &gt; &gt; <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>
&gt; &gt; &gt;<br>
&gt; &gt;<br>
&gt;<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170417/01bd2242/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20170417/01bd2242/<wbr>attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Mon, 17 Apr 2017 15:02:08 -0600<br>
From: Mary Haley &lt;<a href="mailto:haley@ucar.edu">haley@ucar.edu</a>&gt;<br>
Subject: Re: [ncl-talk] Problem with contour plot<br>
To: Guilherme Martins &lt;<a href="mailto:jgmsantos@gmail.com">jgmsantos@gmail.com</a>&gt;<br>
Cc: NCL &lt;<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>&gt;<br>
Message-ID:<br>
        &lt;<a href="mailto:CACNN_CJz1yvpMXMABZ8FhSEov1T0qN63rzsNBBVXBhLos3UKjg@mail.gmail.com">CACNN_<wbr>CJz1yvpMXMABZ8FhSEov1T0qN63rzs<wbr>NBBVXBhLos3UKjg@mail.gmail.com</a><wbr>&gt;<br>
Content-Type: text/plain; charset=&quot;utf-8&quot;<br>
<br>
Guilherme,<br>
<br>
Please see examples shapefiles_21.ncl at:<br>
<br>
<a href="http://www.ncl.ucar.edu/Applications/shapefiles.shtml#ex21" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/<wbr>Applications/shapefiles.shtml#<wbr>ex21</a><br>
<br>
It explains why you see these &quot;blocky&quot; contours in NCL.<br>
<br>
Simply put, if you have four points that represent a box at the edge of a<br>
domain, and three of them have data one them and one is missing, then NCL<br>
will only draw contours in half of the box, because it won&#39;t extrapolate to<br>
that missing point.<br>
<br>
The real question is: do you simply want the graphics not to have gaps, or<br>
do you actually want there to be data in locations that fall just outside<br>
the geographical boundary?  You can do things to force values just outside<br>
the geographical boundary, but this depends on how you are masking in the<br>
first place.<br>
<br>
If you are masking based on a shapefile, then see example shapefiles_18.ncl<br>
on the same page:<br>
<br>
<a href="http://www.ncl.ucar.edu/Applications/shapefiles.shtml#ex18" rel="noreferrer" target="_blank">http://www.ncl.ucar.edu/<wbr>Applications/shapefiles.shtml#<wbr>ex18</a><br>
<br>
Note, however, that this mask is now returning values outside the range of<br>
interest.<br>
<br>
If you are simply trying to mask your data *graphically* by some<br>
geographical area, then instead of trying to use a mask function, you can<br>
draw the full data, and then draw polygons filled in with white in the<br>
areas you don&#39;t want to see graphically.  That&#39;s what example 21 is showing.<br>
<br>
--Mary<br>
<br>
<br>
On Mon, Apr 17, 2017 at 1:00 PM, Guilherme Martins &lt;<a href="mailto:jgmsantos@gmail.com">jgmsantos@gmail.com</a>&gt;<br>
wrote:<br>
<br>
&gt; Dear users,<br>
&gt;<br>
&gt; After searching extensively on the NCL website I was not able to find a<br>
&gt; solution to my problem. I only want values within the contour of the<br>
&gt; continent, but some edges appear in the figure that are stemming from the<br>
&gt; mask that I applied and I do not know how to remove this. The figure is<br>
&gt; attached for a better view of my problem.<br>
&gt;<br>
&gt; Best regards,<br>
&gt;<br>
&gt; Guilherme.<br>
&gt; --<br>
&gt; ------------------------------<wbr>------------------------------<br>
&gt; ------------------------------<wbr>------<br>
&gt; Instituto Nacional de Pesquisas Espaciais (INPE)<br>
&gt; Centro de Previs?o de Tempo e Estudos Clim?ticos (CPTEC)<br>
&gt; Divis?o de Sat?lites e Sistemas Ambientais (DSA)<br>
&gt; Programa de Monitoramento de Queimadas<br>
&gt; Telefone (INPE/CP): <a href="tel:%2B55%2012%203186-9205" value="+551231869205">+55 12 3186-9205</a> &lt;+55%2012%203186-9205&gt; || Celular<br>
&gt; (TIM): <a href="tel:%2B55%2012%2098111-4292" value="+5512981114292">+55 12 98111-4292</a> &lt;+55%2012%2098111-4292&gt;<br>
&gt; E-mail: <a href="mailto:guilherme.martins@inpe.br">guilherme.martins@inpe.br</a> || <a href="mailto:jgmsantos@gmail.com">jgmsantos@gmail.com</a><br>
&gt; Skype: guilherme.martins.<br>
&gt; Homepage: <a href="https://sites.google.com/site/jgmsantos" rel="noreferrer" target="_blank">https://sites.google.com/site/<wbr>jgmsantos</a><br>
&gt; Curr?culo Lattes: <a href="http://lattes.cnpq.br/5997657584785803" rel="noreferrer" target="_blank">http://lattes.cnpq.br/<wbr>5997657584785803</a><br>
&gt; ------------------------------<wbr>------------------------------<br>
&gt; ------------------------------<wbr>------<br>
&gt;<br>
&gt; ______________________________<wbr>_________________<br>
&gt; ncl-talk mailing list<br>
&gt; <a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
&gt; List instructions, subscriber options, unsubscribe:<br>
&gt; <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>
&gt;<br>
&gt;<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170417/96b2b947/attachment-0001.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20170417/96b2b947/<wbr>attachment-0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Tue, 18 Apr 2017 15:29:04 +0800 (CST)<br>
From: wen &lt;<a href="mailto:wenguanhuan@163.com">wenguanhuan@163.com</a>&gt;<br>
Subject: [ncl-talk] how to draw this kind of plot<br>
To: ncl-list &lt;<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>&gt;<br>
Message-ID: &lt;<a href="mailto:5de3ed86.9f7f.15b7ff6fb6f.Coremail.wenguanhuan@163.com">5de3ed86.9f7f.15b7ff6fb6f.<wbr>Coremail.wenguanhuan@163.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;gbk&quot;<br>
<br>
Hi all,<br>
  I am trying to use TRMM 3B42 data  to draw this kind of plot. I haven&#39;t found a example on the website to draw it.  Do anyone know  how to make it. Which function should I<br>
use. Thank you!<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
--<br>
<br>
??<br>
???<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170418/b09be439/attachment.html" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20170418/b09be439/<wbr>attachment.html</a><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: ??1.png<br>
Type: image/png<br>
Size: 249741 bytes<br>
Desc: not available<br>
Url : <a href="http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170418/b09be439/attachment.png" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>pipermail/ncl-talk/<wbr>attachments/20170418/b09be439/<wbr>attachment.png</a><br>
<br>
------------------------------<br>
<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>
<br>
End of ncl-talk Digest, Vol 161, Issue 23<br>
******************************<wbr>***********<br>
</blockquote></div><br></div></div>