<div dir="ltr"><div><div><div><div><div><div><div><div><div><div><div>[1] If your data has missing values, you can not use either specx_anal or the wavelet tool.<br><br>[2] <a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/dtrend.shtml">http://www.ncl.ucar.edu/Document/Functions/Built-in/dtrend.shtml</a><br><br></div>     states: &quot;Missing values (_FillValue) are not allowed.&quot;<br><br></div>[3] If there are **isolated** missing values you could likely do a simple linear interpolation<br></div>     to fill them. However, **any** interpolation will affect the spectra.<br><br></div>     If you have many _FillValue .... I would not even try to fill them in.<br></div>     The resulting spectra will reflect the effects of the interpolation as much as the &#39;rea;&#39; data.<br><br>---<br></div>Punch line: I don&#39;t think you can get what you want ... given your data.<br><br>--<br><br></div>Your statement:<br><br>&quot;My script is from this link: <br><a href="http://paos.colorado.edu/research/wavelets/" target="_blank">http://paos.colorado.edu/research/wavelets/</a><br><br></div>There is no NCL script at that web site.<br><br></div>In fact, NCL call the fortran code presented at the website.<br><br>---<br></div>The interactive code there is IDL ... not NCL<br><br></div>Regards<br><div><div><div><div><div><div>     <br><div><div><div><div><div><div><div>     <br></div></div></div></div></div></div></div></div></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, May 23, 2015 at 9:54 AM, Amadou Coulibaly <span dir="ltr">&lt;<a href="mailto:mpapin24@gmail.com" target="_blank">mpapin24@gmail.com</a>&gt;</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 class="gmail_default" style="font-family:courier new,monospace;font-size:small">Hi Dennis,<br><br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">My script is from this link: <br><a href="http://paos.colorado.edu/research/wavelets/" target="_blank">http://paos.colorado.edu/research/wavelets/</a><br><br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">My Temp look like this:<br><span style="color:rgb(0,0,255)"><br>Variable: Temp<br>Type: float<br>Total Size: 1051968 bytes<br>            262992 values<br>Number of Dimensions: 1<br>Dimensions and sizes:    [time | 262992]<br>Coordinates: <br>            time: [410227200..1356994800]<br>Number Of Attributes: 7<br>  coordinates :    lon lat<br>  grid_type :    unstructured<br>  comment :    Air temperature is the bulk temperature of the air, not the surf<br>ace (skin) temperature.<br>  shortname :    ta<br>  code :      11<br>  units :    K<br>  _FillValue :    9.96921e+36<br>(0)    299.15<br>(1)    9.96921e+36<br>(2)    9.96921e+36<br>(3)    296.35</span><br><br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">It contains a lot of missing data. That&#39;s why I used the following command:<br>; Remove mean and linear trend (recommended, not required)<br>; Missing values are not allowed. <br>;************************************  <br> Temp = dtrend(Temp,False)<br><br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">After that I printed Temp, but look stange for me with negative values:<br><br><span style="color:rgb(0,0,255)">Variable: Temp<br>Type: float<br>Total Size: 1051968 bytes<br>            262992 values<br>Number of Dimensions: 1<br>Dimensions and sizes:    [time | 262992]<br>Coordinates: <br>            time: [410227200..1356994800]<br>Number Of Attributes: 7<br>  coordinates :    lon lat<br>  grid_type :    unstructured<br>  comment :    Air temperature is the bulk temperature of the air, not the surf<br>ace (skin) temperature.<br>  shortname :    ta<br>  code :      11<br>  units :    K<br>  _FillValue :    9.96921e+36<br>(0)    -8.256236e+36<br>(1)    1.713005e+36<br>(2)    1.713035e+36<br>(3)    -8.256145e+36<br></span><br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">I think the problem is from that command trying to remove mean and linear trend. Because without that I can get the wavelet function (w), but it contains only zero as values.<br><br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">How to handle this issue?<br><br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">Best<br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small"><br></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On 23 May 2015 at 14:17, Dennis Shea <span dir="ltr">&lt;<a href="mailto:shea@ucar.edu" target="_blank">shea@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"><div dir="ltr"><div><div><div>[1] I have no idea why the script should seg fault. Maybe your parameter settings are not appropriate for the input data.<br><br></div>[2] What does &quot;Temp&quot; look like? printVarSummary(Temp)<br><br>[3] Are there missing values?<br><br></div>[4] If you want the &quot;power spectrum&quot;, why not look at <br><br><a href="http://www.ncl.ucar.edu/Document/Functions/Built-in/specx_anal.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Built-in/specx_anal.shtml</a><br><a href="http://www.ncl.ucar.edu/Applications/spec.shtml" target="_blank">http://www.ncl.ucar.edu/Applications/spec.shtml</a><br><br></div>Good luck<br></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div>On Sat, May 23, 2015 at 8:00 AM, Amadou Coulibaly <span dir="ltr">&lt;<a href="mailto:mpapin24@gmail.com" target="_blank">mpapin24@gmail.com</a>&gt;</span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div><div dir="ltr"><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">Dear ncl users,<br><br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">When I tried ti run my script (wavelet1.ncl), I am getting a strange error which is: <br><span style="color:rgb(255,0,0)"><span style="background-color:rgb(243,243,243)">Segmentation fault (core dumped)</span></span><br><br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">I am using <span style="color:rgb(255,153,0)">ncl.6.2.1</span> and my system is: <span style="color:rgb(255,153,0)">Linux amadou-Satellite-C855 3.19.0-16-generic #16-Ubuntu SMP Thu Apr 30 16:09:58 UTC 2015 x86_64 x86_64 x86_64 GNU/Linux<br></span><br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">My dataset contains hourly data (1983-2012) of wind(speed/direction), air_temperature, ... for 46 sites. I just wanted to plot the power spectrum of only one site and for the all 30 years. I think the error comes from line 45 when I tried to compute the wavelet function (w). Please find attached my script.<br><br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">How can I solve this issue?<br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small"><br clear="all"></div><br>-- <br><div><div dir="ltr"><b style="font-size:small;font-family:tahoma,sans-serif"><font color="#000066"> </font><font color="#000066">COULIBALY   </font><font color="#000066">AMADOU   </font></b><br><span style="color:rgb(0,0,102);font-family:tahoma,sans-serif;font-size:small">PhD  Student  on  West  African  Climate  System (WACS)</span><br><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">FUTA</font><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066"> - </font><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">F</font><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066">ederal  </font><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">U</font><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066">niversity  of  </font><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">T</font><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066">echnology of </font><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">A</font><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066">kure, Nigeria</font><div><font face="tahoma, sans-serif" color="#000066" size="2"><br></font><div><span style="color:rgb(76,17,48)"><span style="font-family:arial,helvetica,sans-serif"><i style="font-weight:bold">Visiting Student - </i></span><i style="font-weight:bold"><span style="font-family:arial,helvetica,sans-serif"><i><font size="1"><span style="font-size:13px">University of Cologne, Germany</span></font></i></span><br></i><i><span style="font-family:arial,helvetica,sans-serif">Institute of Geophysics and Meteorology</span></i><span style="font-family:arial,helvetica,sans-serif"><br></span></span><div style="letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px"><span style="color:rgb(76,17,48)"><span style="font-family:arial,helvetica,sans-serif"><font size="1">Pohligstr. 3 / Office 3.102</font></span></span></div><div style="letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px"><span style="color:rgb(76,17,48)"><span style="font-family:arial,helvetica,sans-serif"><font size="1">D-50969 Köln</font></span></span></div></div></div><div style="color:rgb(80,0,80)"><i style="color:rgb(34,34,34);font-family:&quot;bookman old style&quot;,&quot;new york&quot;,times,serif"><font size="1"><span style="color:rgb(0,0,127)"><span style="font-size:13px"><br></span></span></font></i></div><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)"><u>Project</u>: </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">WASCAL</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)"> (</span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">W</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">est </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">A</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">frican </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">S</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">cience </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">S</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">ervice </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">C</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">entre on </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">C</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">limate </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">C</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">hange and </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">A</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">dapted </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">L</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">and Use) </span><br><font style="font-family:tahoma,sans-serif;font-size:small;color:rgb(0,0,102)" color="#000066">Phone:<a href="tel:%28%2B234%29%20810%20795%202836" value="+2348107952836" target="_blank">(+234) 810 795 2836</a> </font><font style="font-family:tahoma,sans-serif;font-size:small" color="#ff0000">/</font><font style="font-family:tahoma,sans-serif;font-size:small;color:rgb(0,0,102)" color="#000066"><a href="tel:%28%2B223%29%2065%2067%2027%2057" value="+22365672757" target="_blank">(+223) 65 67 27 57</a> </font><font style="font-family:tahoma,sans-serif;font-size:small" color="#ff0000">/</font><font style="font-family:tahoma,sans-serif;font-size:small;color:rgb(0,0,102)" color="#000066"><a href="tel:%28%2B226%29%2064%2057%2037%2027" value="+22664573727" target="_blank">(+226) 64 57 37 27</a> </font><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066"><span style="color:rgb(34,34,34)"></span><font color="#ff0000">/</font></font><font style="font-family:tahoma,sans-serif;font-size:small;color:rgb(0,0,102)" color="#000066"><a href="tel:%2B49%2015218352574" value="+4915218352574" target="_blank">+49 15218352574</a><br></font><div><p><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066"></font></p><div style="text-align:left"><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066"><font color="#000066">E-mail: </font><a href="mailto:mpapin24@gmail.com" style="color:rgb(0,0,255)" target="_blank">mpapin24@gmail.com</a><span style="color:rgb(34,34,34)"> </span><font color="#ff0000">/ </font><span style="color:rgb(0,0,255)"><a href="mailto:coulibalya68@yahoo.com" target="_blank">coulibalya68@yahoo.com</a> </span></font></div><p></p></div><div><font face="verdana,sans-serif" size="2"><p><font color="#ff9900">&quot;<i>The time is always right to do right&quot;: <b>Nelson Mandela</b></i></font></p><p><i><font color="#ff9900">&quot;Character is like a tree and reputation like a shadow.  The shadow is what we think of it; the tree is the real thing&quot; : <b>Abraham Lincoln</b></font></i></p><p><i><font color="#ff9900">&quot;Do what you can, with what you have, where you are&quot; <b>Theodore Roosevelt</b></font></i></p></font><p><b><font color="#000066" size="2"><i> </i></font></b></p><p><b><font color="#000066" size="4"> </font></b></p><div> </div></div></div></div>
</div>
<br></div></div>_______________________________________________<br>
ncl-talk mailing list<br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>
</blockquote></div><br><br clear="all"><br>-- <br><div><div dir="ltr"><b style="font-size:small;font-family:tahoma,sans-serif"><font color="#000066"> </font><font color="#000066">COULIBALY   </font><font color="#000066">AMADOU   </font></b><br><span style="color:rgb(0,0,102);font-family:tahoma,sans-serif;font-size:small">PhD  Student  on  West  African  Climate  System (WACS)</span><br><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">FUTA</font><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066"> - </font><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">F</font><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066">ederal  </font><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">U</font><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066">niversity  of  </font><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">T</font><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066">echnology of </font><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">A</font><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066">kure, Nigeria</font><div><font face="tahoma, sans-serif" color="#000066" size="2"><br></font><div><span style="color:rgb(76,17,48)"><span style="font-family:arial,helvetica,sans-serif"><i style="font-weight:bold">Visiting Student - </i></span><i style="font-weight:bold"><span style="font-family:arial,helvetica,sans-serif"><i><font size="1"><span style="font-size:13px">University of Cologne, Germany</span></font></i></span><br></i><i><span style="font-family:arial,helvetica,sans-serif">Institute of Geophysics and Meteorology</span></i><span style="font-family:arial,helvetica,sans-serif"><br></span></span><div style="letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px"><span style="color:rgb(76,17,48)"><span style="font-family:arial,helvetica,sans-serif"><font size="1">Pohligstr. 3 / Office 3.102</font></span></span></div><div style="letter-spacing:normal;text-transform:none;white-space:normal;word-spacing:0px"><span style="color:rgb(76,17,48)"><span style="font-family:arial,helvetica,sans-serif"><font size="1">D-50969 Köln</font></span></span></div></div></div><div style="color:rgb(80,0,80)"><i style="color:rgb(34,34,34);font-family:&#39;bookman old style&#39;,&#39;new york&#39;,times,serif"><font size="1"><span style="color:rgb(0,0,127)"><span style="font-size:13px"><br></span></span></font></i></div><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)"><u>Project</u>: </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">WASCAL</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)"> (</span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">W</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">est </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">A</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">frican </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">S</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">cience </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">S</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">ervice </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">C</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">entre on </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">C</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">limate </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">C</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">hange and </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">A</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">dapted </span><font style="font-size:small;font-family:tahoma,sans-serif" color="#0000ff">L</font><span style="font-size:small;font-family:tahoma,sans-serif;color:rgb(0,0,102)">and Use) </span><br><font style="font-family:tahoma,sans-serif;font-size:small;color:rgb(0,0,102)" color="#000066">Phone:<a href="tel:%28%2B234%29%20810%20795%202836" value="+2348107952836" target="_blank">(+234) 810 795 2836</a> </font><font style="font-family:tahoma,sans-serif;font-size:small" color="#ff0000">/</font><font style="font-family:tahoma,sans-serif;font-size:small;color:rgb(0,0,102)" color="#000066"><a href="tel:%28%2B223%29%2065%2067%2027%2057" value="+22365672757" target="_blank">(+223) 65 67 27 57</a> </font><font style="font-family:tahoma,sans-serif;font-size:small" color="#ff0000">/</font><font style="font-family:tahoma,sans-serif;font-size:small;color:rgb(0,0,102)" color="#000066"><a href="tel:%28%2B226%29%2064%2057%2037%2027" value="+22664573727" target="_blank">(+226) 64 57 37 27</a> </font><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066"><span style="color:rgb(34,34,34)"></span><font color="#ff0000">/</font></font><font style="font-family:tahoma,sans-serif;font-size:small;color:rgb(0,0,102)" color="#000066"><a href="tel:%2B49%2015218352574" value="+4915218352574" target="_blank">+49 15218352574</a><br></font><div><p><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066"></font></p><div style="text-align:left"><font style="font-size:small;font-family:tahoma,sans-serif" color="#000066"><font color="#000066">E-mail: </font><a href="mailto:mpapin24@gmail.com" style="color:rgb(0,0,255)" target="_blank">mpapin24@gmail.com</a><span style="color:rgb(34,34,34)"> </span><font color="#ff0000">/ </font><span style="color:rgb(0,0,255)"><a href="mailto:coulibalya68@yahoo.com" target="_blank">coulibalya68@yahoo.com</a> </span></font></div><p></p></div><div><font face="verdana,sans-serif" size="2"><p><font color="#ff9900">&quot;<i>The time is always right to do right&quot;: <b>Nelson Mandela</b></i></font></p><p><i><font color="#ff9900">&quot;Character is like a tree and reputation like a shadow.  The shadow is what we think of it; the tree is the real thing&quot; : <b>Abraham Lincoln</b></font></i></p><p><i><font color="#ff9900">&quot;Do what you can, with what you have, where you are&quot; <b>Theodore Roosevelt</b></font></i></p></font><p><b><font color="#000066" size="2"><i> </i></font></b></p><p><b><font color="#000066" size="4"> </font></b></p><div> </div></div></div></div>
</div>
</div></div></blockquote></div><br></div>