<div dir="ltr"><div dir="ltr">Shouli Xuan,<br><br>Documentation for bw_bandpass_filter says it is optimized for narrow bandwidths. Low pass is far outside the intended purpose of this function. I looked at the reference paper (David Russell, 2005) and the source code. It looks like this function may be unstable if you try to use it as a low pass filter. Therefore you should find alternate code for low pass Butterworth filter, or else use a different filter.</div><div dir="ltr"><br></div><div>The source code applies fca and fcb in a symmetric way. This is why exchanging them gives the same results. However, this is completely accidental, and outside the specifications for this function. In general it is a mistake to use function parameters outside of their intended constraints <span style="color:rgb(51,51,51);font-family:verdana,sans-serif;font-size:13.3333px">(fcb > fca in this case).</span></div><div><br></div><div>I want to reinforce Dennis's suggestion. Compare any special filter code with the NCL Lanczos low pass filter. Check the filter responses in both time and frequency domains, as shown in examples filters_8.ncl and filters_9.ncl.</div><div><br></div><div>--Dave</div><div dir="ltr"><br></div><div dir="ltr"><br><div class="gmail_quote"><div dir="ltr">On Sun, Jan 13, 2019 at 11:46 AM Dennis Shea <<a href="mailto:shea@ucar.edu">shea@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"><div dir="ltr"><div>[1] DATA: The years have changed but the files are available at:<br></div><div><br></div><div><div dir="ltr"><div dir="ltr"><div>ftp <a href="http://ftp.cgd.ucar.edu" target="_blank">ftp.cgd.ucar.edu</a></div><div>anonymous</div><div>your_email</div><div>cd pub/shea/MJO</div><div>mget * .<br></div><div>quit</div><div>==========</div><div><div><br></div><div>Note: The file names have changed slightly from those used in the examples. For example:<br></div><div><pre><a href="http://olr.day.anomalies.1980-2005.nc" target="_blank">olr.day.anomalies.1980-2005.nc</a> ==> <a href="http://olr.day.anomalies.1980-2011.nc" target="_blank">olr.day.anomalies.1980-2011.nc</a></pre></div></div></div><div>Perhaps, some minor changes in the graphic resources will be needed. <br></div><div dir="ltr">==========</div><div dir="ltr">/ftp/pub/shea/MJO<br>-rw-r--r-- 1 shea cgdcas 458055500 Jan 11 10:08 <a href="http://MODEL.SST.HAD187001-198110.OI198111-201703.nc" target="_blank">MODEL.SST.HAD187001-198110.OI198111-201703.nc</a><br>-rw-r--r-- 1 shea cgdcas 30699392 Jan 11 10:18 <a href="http://OLR.12hr_2yrs.wheeler.nc" target="_blank">OLR.12hr_2yrs.wheeler.nc</a><br>-rw-r--r-- 1 shea cgdcas 858482293 Jan 11 09:40 olr.day.anomalies.1979-2011.nc.gz<br>-rw-r--r-- 1 shea cgdcas 304335976 Jan 11 10:10 <a href="http://olr.day.mean.nc" target="_blank">olr.day.mean.nc</a><br>-rw-r--r-- 1 shea cgdcas 615513004 Jan 11 10:24 <a href="http://pregpcp19962008.daily.nc" target="_blank">pregpcp19962008.daily.nc</a><br>-rw-r--r-- 1 shea cgdcas 350655700 Jan 11 10:23 <a href="http://TRMM.200006-200109.thin.nc" target="_blank">TRMM.200006-200109.thin.nc</a><br>-rw-r--r-- 1 shea cgdcas 929887809 Jan 11 09:41 uwnd.day.anomalies.1979-2011.200.nc.gz<br>-rw-r--r-- 1 shea cgdcas 934253770 Jan 11 09:41 uwnd.day.anomalies.1979-2011.850.nc.gz<br>-rw-r--r-- 1 shea cgdcas 939714739 Jan 11 09:41 vwnd.day.anomalies.1979-2011.200.nc.gz<br>-rw-r--r-- 1 shea cgdcas 939092767 Jan 11 09:41 vwnd.day.anomalies.1979-2011.850.nc.gz<br></div>====</div><div>The *gz files will have to be unpacked:</div><div> %> gzip -d *nc.gz</div><div>====</div><div><br></div><div>FILTERS:</div><div> <a href="http://www.ncl.ucar.edu/Applications/filter.shtml" target="_blank"><b>http://www.ncl.ucar.edu/Applications/filter.shtml</b></a></div><div><br></div><div>There are several examples of Butterworth filters. <br></div><div>These included comaprison with Lanczos filters.</div><div><br></div><div>I do not know why: ".. exchange values of ca and cb .." does not change the results.</div><div>I don't have time to look right now. <br></div><div>Perhaps the interface between NCL and the ubderlying fortran code does an exvhange.</div><div><br></div><div>re: making the band pass into a low-pass filter. Pergform xome comparisons withe a Lanczos low pass filter.</div><div><br></div><div>Good Luck<br></div><div><br></div></div></div></div><br><div class="gmail_quote"><div dir="ltr">On Sat, Jan 12, 2019 at 8:27 PM 15295538792 <<a href="mailto:15295538792@163.com" target="_blank">15295538792@163.com</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"><u></u>
<div style="border-width:0px;font-size:10.5pt;font-family:微软雅黑;color:rgb(0,0,0);margin:12px;line-height:1.5"><u></u>
<div>Dear NCL experts,</div>
<div> </div>
<div>I'm writing to ask if you could send me the data file named
"<a href="http://uwnd.day.850.anomalies.1980-2005.nc" target="_blank">uwnd.day.850.anomalies.1980-2005.nc</a>" which is used in the script
"bfband_2.ncl". I can't find it in the online data files. In addition, I'd like
to consult you about the function "bw_bandpass_filter" which applies a
Butterworth bandpass filter to time series. The online documentation says that
the cut-off paramaters ca and cb of this function should meet a condition of
ca>cb (fcb>fca). However, when I exchanged the values of ca and cb, the
results remained the same. How could this happen? In fact, I want to use this
bandpass filter as a low pass filter which need only one cut-off frequency. So
which parameter should I set to arbitrarily large (ca or cb)? Is there anything
wrong in my script? Please help. Thanking in anticipation. (PS: Input and
output data files as well as ncl script I used has been attached to this
email.)</div>
<div> </div>
<div> </div>
<div>Sincerely,<br>Shouli Xuan<br><a href="mailto:Email%3A15295538792@163.com" target="_blank">Email:15295538792@163.com</a></div>
<div> </div>
<div style="font-size:10pt;font-family:Verdana;color:rgb(192,192,192)">
<div align="left">2019-01-13</div>
<hr id="gmail-m_-8398776638075760935gmail-m_5506818083859252814SignNameHR" style="border-color:rgb(192,192,192) currentcolor currentcolor;border-style:solid none none;border-width:1px 0px 0px;height:1px;width:122px" align="left">
<span id="gmail-m_-8398776638075760935gmail-m_5506818083859252814_FlashSignName">15295538792</span></div></div></blockquote></div>
</blockquote></div></div></div>