<div dir="ltr"><div><div><div><div><div><div>Looking at scripts can be a tedious and time consuming task. <br></div>--<br><br>Your code <br><br> infile = addfile("/home/amadou/NCL/<a href="http://all_new.nc">all_new.nc</a>", "r")<br> temp = infile->air_temperature<br> Temp = temp(256392:258575,7) ; for 5 years (2012_JFMAMJ)<br> time = infile->time<br> Time = time(256392:258575)<br><br> index = ind(Temp(:) .ne. Temp@_FillValue) ; indices of Temp without missing values<br> printVarSummary(index)<br> Time0 = time(256444:258575)<br> Temp1 = new((/2132/), float)<br> Temp1 = Temp(index)<br><br></div>[snip]<br><br></div>[1] You can **not** use the following to check for missing values. These can *only* be checked via the 'ismissing' function.<br><br> index = ind(Temp(:) .ne. Temp@_FillValue) ; indices of Temp without missing values<br><br></div> Use<br><br></div> index = ind(.not.ismissing(Temp)) ; indices of Temp without missing values<br><br> <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/ismissing.shtml">https://www.ncl.ucar.edu/Document/Functions/Built-in/ismissing.shtml</a><br> <a href="https://www.ncl.ucar.edu/Document/Functions/Built-in/ind.shtml">https://www.ncl.ucar.edu/Document/Functions/Built-in/ind.shtml</a><br><br></div>[2] <br>As I mentioned to you in a previous ncl-talk response, you can not just eliminate missing values and perform a wavelet analysis on the 'good' values. Omitting values will result in 'numbers' but they are not the 'correct numbers'. Omitting values results in spectral gaps. If you have (say) hourly values and very short spans with missing values you could perform a linear interpolation. This will affect the spectra but likely not in any substantive manner.<br><div><div><div><br>[3] <br></div><div>You have two different sets of index values: 256392:258575 and 256444:258575 ... this is likely the source of your plotting issue.<br><br></div><div>Good programming practice suggest that you use variables rather than hard wired constants.<br><br></div><div> tStrt = 256392<br></div><div> tLast = 258575<br><br></div><div>and use these<br><br> temp = infile->air_temperature<br> Temp = temp(tStrt:tLast,7) ; for 5 years (2012_JFMAMJ)<br> time = infile->time<br> Time = time(tStrt:tLast)<br><br></div><div>Actually .... better<br><br></div><div> Temp = infile->air_temperature(tStrt:tLast,7) <br></div><div> printVarSummary(Temp)<br><br><br></div><div><br></div><div><br></div><div><br><div><br><div><br><br></div></div></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jun 3, 2015 at 6:00 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">Dear NCL users,<br><br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">I tried to do a panel plot between "<span style="color:rgb(0,0,255)">Wavelet Power Spectrum</span>" and "<span style="color:rgb(0,0,255)">Timesires</span>" from hourly temperature data of April-May-June. But this data contained 61 hours of missing values that I cancelled. When I plotted my data, I remarked that the x-axis
on the wavelet power spectrum does not match the timeserie one. <br><br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">What is wrong with my script?<br><br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">Please find attached my script and my plot.<br><br></div><div class="gmail_default" style="font-family:courier new,monospace;font-size:small">Best <br></div>-- <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>_______________________________________________<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" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>