<div dir="ltr"><div><div>Hi,<br><br></div>It appears to work for me; if I uncomment the printVarSummary(w), I get the output below. Do the values look reasonable? Which version of NCL are you using? Perhaps an upgrade to v6.4.0 would solve the issue<br><br></div>Rick<br><div><br>$ ncl foo.ncl<br> Copyright (C) 1995-2017 - All Rights Reserved<br> University Corporation for Atmospheric Research<br> NCAR Command Language Version 6.4.1<br> The use of this software is governed by a License Agreement.<br> See <a href="http://www.ncl.ucar.edu/">http://www.ncl.ucar.edu/</a> for more details.<br><br><br>Variable: ntime<br>Type: integer<br>Total Size: 4 bytes<br> 1 values<br>Number of Dimensions: 1<br>Dimensions and sizes: [1]<br>Coordinates: <br>(0) 504<br><br>Variable: w<br>Type: float<br>Total Size: 145152 bytes<br> 36288 values<br>Number of Dimensions: 3<br>Dimensions and sizes: [2] x [36] x [504]<br>Coordinates: <br>Number Of Attributes: 15<br> psi0 : 0.7511256<br> cdelta : 0.776<br> lag1 : 0.7689273<br> stdev : 0.7335985<br> mean : -4.20585e-09<br> r1 : 0.7689273<br> phase : <ARRAY of 18144 elements><br> power : <ARRAY of 18144 elements><br> gws : <ARRAY of 36 elements><br> signif : <ARRAY of 36 elements><br> fft_theor : <ARRAY of 36 elements><br> dof : <ARRAY of 36 elements><br> coi : <ARRAY of 504 elements><br> period : <ARRAY of 36 elements><br> scale : <ARRAY of 36 elements><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Apr 18, 2017 at 9:39 AM, Dipti Sharma <span dir="ltr"><<a href="mailto:sdipti596@gmail.com" target="_blank">sdipti596@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"><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 "segmentation fault (core dumped)" 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 "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_code.ncl"<br> load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/gsn_csm.ncl"<br> load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/contributed.<wbr>ncl"<br> load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/shea_util.ncl"<br> load "$NCARG_ROOT/lib/ncarg/<wbr>nclscripts/csm/shea_util.ncl"<br><br>;-----------------------------<wbr>------------------------------<wbr>---------<br>; CREATE SAME WAVELET FIGURE AS TORRENCE & COMPO using seasonal Nino3<br>; from 1871-1997<br>begin<br>ninoseas = asciiread("C:\diptic\sst\sst_<wbr>nino3.dat",-1,"float")<br> ninoseas!0 = "time"<br> ntime = dimsizes(ninoseas)<br> timeo = fspan(1871.25,1996.,ntime)<br> print(ntime)<br> ninoseas&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>;*****************************<wbr>*******<br>; compute wavelet<br>;*****************************<wbr>*******<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/<wbr>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,<wbr>param,s0,dj,jtot,npad,noise,<wbr>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"><<a href="mailto:ncl-talk-request@ucar.edu" target="_blank">ncl-talk-request@ucar.edu</a>></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" target="_blank">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/mailma<wbr>n/listinfo/ncl-talk</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:ncl-talk-request@ucar.edu" target="_blank">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" target="_blank">ncl-talk-owner@ucar.edu</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of ncl-talk digest..."<br>
<br>
<br>
Today'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 <<a href="mailto:sdafis@cc.uoi.gr" target="_blank">sdafis@cc.uoi.gr</a>><br>
Subject: Re: [ncl-talk] Problem with gsn_add_text<br>
To: Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>><br>
Cc: "<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Message-ID: <<a href="mailto:1492454911.58f50dffc9203@webmail.uoi.gr" target="_blank">1492454911.58f50dffc9203@webm<wbr>ail.uoi.gr</a>><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>
<a href="tel:+33%209%2081%2094%2022%2012" value="+33981942212" target="_blank">(+33)9 81 94 22 12</a>, 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 <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>>:<br>
<br>
> Hi Stavros,<br>
><br>
> I don't know for certain, but I suspect you need to assign the _FillValue<br>
> to the string array "points" that you constructed from "t2" (I don't see<br>
> how the attributes from t2 would carry to points via the sprintf() ).<br>
><br>
> Hope that helps...<br>
> Rick<br>
><br>
> On Sun, Apr 16, 2017 at 12:56 PM, Stavros Dafis <<a href="mailto:sdafis@cc.uoi.gr" target="_blank">sdafis@cc.uoi.gr</a>> wrote:<br>
><br>
> > Hello everyone, I am facing a problem with the function gsn_add_text. In<br>
> > the<br>
> > description of the function, it is mentioned that no missing values will be<br>
> > plotted:<br>
> ><br>
> > &#8220;This function attaches a list of text strings to the given plot.<br>
> > The text<br>
> > strings will be drawn only when the plot is drawn. If a missing value is<br>
> > encountered in the text strings or x/y values, then no text string will be<br>
> > drawn at this pair.&#8221;<br>
> ><br>
> > My script reads a netCDF file with 3 variables t2(lat,lon), lat(lat),<br>
> > lon(lon),<br>
> > so in order to attach values on the plot using this function, I have to<br>
> > make<br>
> > new 2-D variables for lat/lon and define the string to be plotted on the<br>
> > map.<br>
> > Values are assigned just fine on the plot, over the land, but missing<br>
> > values<br>
> > are also printed over the sea. I tried to use landsea masking but my data<br>
> > have<br>
> > high horizontal resolution (5km). Any suggestions?<br>
> ><br>
> > I have the same problem with WRF data, with already 2-D lat/lon arrays.<br>
> ><br>
> > This is the part of the code which prints the values:<br>
> > .<br>
> > .<br>
> > .<br>
> > contour = gsn_csm_contour_map(wks,t2,res<wbr>)<br>
> ><br>
> > t2@missing_value = -1<br>
> > t2@_FillValue = -1<br>
> > (default missing values = -999)<br>
> ><br>
> > lat2 = new(dimsizes(t2),typeof(lat),-<wbr>1)<br>
> > lon2 = new(dimsizes(t2),typeof(lon),-<wbr>1)<br>
> ><br>
> > nb = dimsizes(t2)<br>
> ><br>
> > do nl=0,nb(0)-1,1<br>
> > do ml=0,nb(1)-1,1<br>
> ><br>
> > lat2(nl,ml) = lat(nl)<br>
> > lon2(nl,ml) = lon(ml)<br>
> ><br>
> > end do<br>
> > end do<br>
> ><br>
> > points = sprintf("%3.0f", t2)<br>
> > ;points = ""+t2(:,:)+""<br>
> ><br>
> > values = gsn_add_text(wks,contour,point<wbr>s,lon2(:,:),lat2(:,:),True)<br>
> ><br>
> ><br>
> ><br>
> > --<br>
> > Stavros NTAFIS (DAFIS)<br>
> > ------------------------------<wbr>-----------------<br>
> > Physicist - Meteorologist, M.Sc.<br>
> > Research Associate, National Observatory of Athens<br>
> > <a href="tel:+33%209%2081%2094%2022%2012" value="+33981942212" target="_blank">(+33)9 81 94 22 12</a>, 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>
> ><br>
> > ______________________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
> ><br>
><br>
<br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Mon, 17 Apr 2017 19:00:02 +0000<br>
From: Guilherme Martins <<a href="mailto:jgmsantos@gmail.com" target="_blank">jgmsantos@gmail.com</a>><br>
Subject: [ncl-talk] Problem with contour plot<br>
To: NCL <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>
<CAH+t_hMaRUDHM0wMkLC9cz2pFmMY<wbr>GXgGC+X+bcxDLo=+<a href="mailto:8fgWQg@mail.gmail.com" target="_blank">8fgWQg@mail.<wbr>gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<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" target="_blank">+55 12 3186-9205</a> || Celular (TIM): <a href="tel:%2B55%2012%2098111-4292" value="+5512981114292" target="_blank">+55 12 98111-4292</a><br>
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><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/59976575<wbr>84785803</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/piperm<wbr>ail/ncl-talk/attachments/<wbr>20170417/c3215b0e/attachment-<wbr>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/piperm<wbr>ail/ncl-talk/attachments/<wbr>20170417/c3215b0e/attachment-<wbr>0001.png</a><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Mon, 17 Apr 2017 14:32:43 -0600<br>
From: Rick Brownrigg <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Problem with gsn_add_text<br>
To: Stavros Dafis <<a href="mailto:sdafis@cc.uoi.gr" target="_blank">sdafis@cc.uoi.gr</a>><br>
Cc: "<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>
<CAGKRhbHX1v-1sJ3G=<a href="mailto:QrXGA%2B0wSJ2KKsBB_7GqbCnp9krh9fkkQ@mail.gmail.com" target="_blank">QrXGA+0wSJ2<wbr>KKsBB_7GqbCnp9krh9fkkQ@mail.<wbr>gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Hi,<br>
<br>
You used -1 as _FillValue for variable "t2". 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 "points" from t2 as "%3.0f", so that any data values<br>
in t2 that were -1 are now " -1", and thus it seems your _FillValue should<br>
be " -1". Did you do this, and do the types between data and _FillValue<br>
agree everywhere?<br>
<br>
If that'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 <<a href="mailto:sdafis@cc.uoi.gr" target="_blank">sdafis@cc.uoi.gr</a>> wrote:<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<br>
> 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>
> <a href="tel:+33%209%2081%2094%2022%2012" value="+33981942212" target="_blank">(+33)9 81 94 22 12</a>, 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 <<a href="mailto:brownrig@ucar.edu" target="_blank">brownrig@ucar.edu</a>>:<br>
><br>
> > Hi Stavros,<br>
> ><br>
> > I don't know for certain, but I suspect you need to assign the _FillValue<br>
> > to the string array "points" that you constructed from "t2" (I don't see<br>
> > how the attributes from t2 would carry to points via the sprintf() ).<br>
> ><br>
> > Hope that helps...<br>
> > Rick<br>
> ><br>
> > On Sun, Apr 16, 2017 at 12:56 PM, Stavros Dafis <<a href="mailto:sdafis@cc.uoi.gr" target="_blank">sdafis@cc.uoi.gr</a>><br>
> wrote:<br>
> ><br>
> > > Hello everyone, I am facing a problem with the function gsn_add_text.<br>
> In<br>
> > > the<br>
> > > description of the function, it is mentioned that no missing values<br>
> will be<br>
> > > plotted:<br>
> > ><br>
> > > &#8220;This function attaches a list of text strings to the given plot.<br>
> > > The text<br>
> > > strings will be drawn only when the plot is drawn. If a missing value<br>
> is<br>
> > > encountered in the text strings or x/y values, then no text string<br>
> will be<br>
> > > drawn at this pair.&#8221;<br>
> > ><br>
> > > My script reads a netCDF file with 3 variables t2(lat,lon), lat(lat),<br>
> > > lon(lon),<br>
> > > so in order to attach values on the plot using this function, I have to<br>
> > > make<br>
> > > new 2-D variables for lat/lon and define the string to be plotted on<br>
> the<br>
> > > map.<br>
> > > Values are assigned just fine on the plot, over the land, but missing<br>
> > > values<br>
> > > are also printed over the sea. I tried to use landsea masking but my<br>
> data<br>
> > > have<br>
> > > high horizontal resolution (5km). Any suggestions?<br>
> > ><br>
> > > I have the same problem with WRF data, with already 2-D lat/lon arrays.<br>
> > ><br>
> > > This is the part of the code which prints the values:<br>
> > > .<br>
> > > .<br>
> > > .<br>
> > > contour = gsn_csm_contour_map(wks,t2,res<wbr>)<br>
> > ><br>
> > > t2@missing_value = -1<br>
> > > t2@_FillValue = -1<br>
> > > (default missing values = -999)<br>
> > ><br>
> > > lat2 = new(dimsizes(t2),typeof(lat),-<wbr>1)<br>
> > > lon2 = new(dimsizes(t2),typeof(lon),-<wbr>1)<br>
> > ><br>
> > > nb = dimsizes(t2)<br>
> > ><br>
> > > do nl=0,nb(0)-1,1<br>
> > > do ml=0,nb(1)-1,1<br>
> > ><br>
> > > lat2(nl,ml) = lat(nl)<br>
> > > lon2(nl,ml) = lon(ml)<br>
> > ><br>
> > > end do<br>
> > > end do<br>
> > ><br>
> > > points = sprintf("%3.0f", t2)<br>
> > > ;points = ""+t2(:,:)+""<br>
> > ><br>
> > > values = gsn_add_text(wks,contour,point<wbr>s,lon2(:,:),lat2(:,:),True)<br>
> > ><br>
> > ><br>
> > ><br>
> > > --<br>
> > > Stavros NTAFIS (DAFIS)<br>
> > > ------------------------------<wbr>-----------------<br>
> > > Physicist - Meteorologist, M.Sc.<br>
> > > Research Associate, National Observatory of Athens<br>
> > > <a href="tel:+33%209%2081%2094%2022%2012" value="+33981942212" target="_blank">(+33)9 81 94 22 12</a>, 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>
> > ><br>
> > > ______________________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
> > ><br>
> ><br>
><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/piperm<wbr>ail/ncl-talk/attachments/<wbr>20170417/01bd2242/attachment-<wbr>0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Mon, 17 Apr 2017 15:02:08 -0600<br>
From: Mary Haley <<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>><br>
Subject: Re: [ncl-talk] Problem with contour plot<br>
To: Guilherme Martins <<a href="mailto:jgmsantos@gmail.com" target="_blank">jgmsantos@gmail.com</a>><br>
Cc: NCL <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Message-ID:<br>
<<a href="mailto:CACNN_CJz1yvpMXMABZ8FhSEov1T0qN63rzsNBBVXBhLos3UKjg@mail.gmail.com" target="_blank">CACNN_CJz1yvpMXMABZ8FhSEov1T0<wbr>qN63rzsNBBVXBhLos3UKjg@mail.<wbr>gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<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/Applic<wbr>ations/shapefiles.shtml#ex21</a><br>
<br>
It explains why you see these "blocky" 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'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/Applic<wbr>ations/shapefiles.shtml#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't want to see graphically. That's what example 21 is showing.<br>
<br>
--Mary<br>
<br>
<br>
On Mon, Apr 17, 2017 at 1:00 PM, Guilherme Martins <<a href="mailto:jgmsantos@gmail.com" target="_blank">jgmsantos@gmail.com</a>><br>
wrote:<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>------------------------------<br>
> ------------------------------<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" target="_blank">+55 12 3186-9205</a> <+55%2012%203186-9205> || Celular<br>
> (TIM): <a href="tel:%2B55%2012%2098111-4292" value="+5512981114292" target="_blank">+55 12 98111-4292</a> <+55%2012%2098111-4292><br>
> 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><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/59976575<wbr>84785803</a><br>
> ------------------------------<wbr>------------------------------<br>
> ------------------------------<wbr>------<br>
><br>
> ______________________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
><br>
><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/piperm<wbr>ail/ncl-talk/attachments/<wbr>20170417/96b2b947/attachment-<wbr>0001.html</a><br>
<br>
------------------------------<br>
<br>
Message: 5<br>
Date: Tue, 18 Apr 2017 15:29:04 +0800 (CST)<br>
From: wen <<a href="mailto:wenguanhuan@163.com" target="_blank">wenguanhuan@163.com</a>><br>
Subject: [ncl-talk] how to draw this kind of plot<br>
To: ncl-list <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
Message-ID: <<a href="mailto:5de3ed86.9f7f.15b7ff6fb6f.Coremail.wenguanhuan@163.com" target="_blank">5de3ed86.9f7f.15b7ff6fb6f.Cor<wbr>email.wenguanhuan@163.com</a>><br>
Content-Type: text/plain; charset="gbk"<br>
<br>
Hi all,<br>
I am trying to use TRMM 3B42 data to draw this kind of plot. I haven'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/piperm<wbr>ail/ncl-talk/attachments/<wbr>20170418/b09be439/attachment.<wbr>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/piperm<wbr>ail/ncl-talk/attachments/<wbr>20170418/b09be439/attachment.<wbr>png</a><br>
<br>
------------------------------<br>
<br>
______________________________<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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br>
<br>
End of ncl-talk Digest, Vol 161, Issue 23<br>
******************************<wbr>***********<br>
</blockquote></div><br></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>