<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: "Missing values (_FillValue) are not allowed."<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 'rea;' data.<br><br>---<br></div>Punch line: I don't think you can get what you want ... given your data.<br><br>--<br><br></div>Your statement:<br><br>"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"><<a href="mailto:mpapin24@gmail.com" target="_blank">mpapin24@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 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'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"><<a href="mailto:shea@ucar.edu" target="_blank">shea@ucar.edu</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>[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 "Temp" look like? printVarSummary(Temp)<br><br>[3] Are there missing values?<br><br></div>[4] If you want the "power spectrum", 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"><<a href="mailto:mpapin24@gmail.com" target="_blank">mpapin24@gmail.com</a>></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:"bookman old style","new york",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">"<i>The time is always right to do right": <b>Nelson Mandela</b></i></font></p><p><i><font color="#ff9900">"Character is like a tree and reputation like a shadow. The shadow is what we think of it; the tree is the real thing" : <b>Abraham Lincoln</b></font></i></p><p><i><font color="#ff9900">"Do what you can, with what you have, where you are" <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:'bookman old style','new york',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">"<i>The time is always right to do right": <b>Nelson Mandela</b></i></font></p><p><i><font color="#ff9900">"Character is like a tree and reputation like a shadow. The shadow is what we think of it; the tree is the real thing" : <b>Abraham Lincoln</b></font></i></p><p><i><font color="#ff9900">"Do what you can, with what you have, where you are" <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>