<div dir="ltr"><div>ncl-talk has no idea what you have done. The series has <b> *pentad*</b> [5-day] temporal sampling. The plots have the frequency axis labeled 'cycles/month'. <br></div><div><br></div><div>ncl 0> f = addfile("<a href="http://rainfall_pc_00000.nc">rainfall_pc_00000.nc</a>","r")<br>ncl 1> ymd = cd_calendar(f->time,-2)<br>ncl 2> print(ymd)</div><div><br></div><div>(0)     19810303<br>(1)     19810308<br>(2)     19810313<br>(3)     19810318<br>(4)     19810323<br>(5)     19810328<br>(6)     19810402<br>(7)     19810407<br>(8)     19810412<br>(9)     19810417<br>(10)    19810422<br>(11)    19810427<br>(12)    19810502<br>(13)    19810507</div><div>.</div><div>.</div><div>It is <b>your responsibility </b>[**not ncl-talk's**]  to design a filter that addresses the pentad time sampling to create a (say) 20-60 day filter. <br></div><div>ncl-talk's experience is that dealing with this type of issue leads to many emails. We can no longer do that.</div><div><br></div><div>Good Luck<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Sep 15, 2020 at 6:43 AM Edwin Marco via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Dear ncl users,<div><br></div><div>I am trying to prove whether I have correctly filtered my dataset by plotting the power spectrum of the filtered data set.</div><div><br></div><div>My approach: </div><div>firstly I have plotted the spectrum of the unfiltered data set (<a href="http://rainfall_pc_00000.nc" target="_blank">rainfall_pc_00000.nc</a>) using attached script (plot_spectrum_unfiltered_data.ncl).</div><div><br></div><div>Secondly, I decided to filter my dataset at 20-60 days.</div><div><br></div><div>Thirdly; Using the Butterworth filtering technique, filtered the dataset (<a href="http://rainfall_pc_00000.nc" target="_blank">rainfall_pc_00000.nc</a>) at 20-60 days, and plot the spectrum of the filtered dataset.</div><div><br></div><div>Problem: After plotting the spectrum of the filtered dataset, I intended to see all peaks out of 60-20days bandpass to disappear however it came the opposite way.</div><div><br></div><div>I have attached scripts, plots, and datasets for easy follow-up. </div><div><br></div><div>Please help me know if my approach to proof whether I have correctly filtered the datasets is correctly and then where I am getting wrong and plot opposite results, contrary to the theory.</div><div><br></div><div>Thanks </div><div><br></div><div>Edwin</div><div><br></div><div><br></div><div> </div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>