<div dir="ltr"><div>Oi Marilia<br><br></div>[SNIP]<br><div><br>;*****************************<wbr>*******************<br>; create BandPass Filter<br>;*****************************<wbr>*******************<br> ihp = 2 <wbr> ; bpf=>band pass filter => (BPF)<br> nWgt = 201<br> sigma = 1.0 <wbr> ; Lanczos sigma<br> fca = 1./100.<br> fcb = 1./20.<br> wgt = filwgts_lanczos (nWgt, ihp, fca, fcb, sigma )<br><br>;*****************************<wbr>******************************<br>; Read daily U anomalies<br>;*****************************<wbr>******************************<br><br> fili = "uwnd.day."+plev+".<a href="http://anomalies.1980-2005.nc" target="_blank">anomalies.<wbr>1980-2005.nc</a>"<br> f = addfile (diri+fili , "r")<br> TIME = f->time <wbr> ; days since ...<br> YMD = ut_calendar(TIME, -2) ; entire (time,6)<br><br> iStrt = ind(YMD.eq.ymdStrt) <wbr> ; index start<br> iLast = ind(YMD.eq.ymdLast) <wbr> ; index last <br> delete(TIME)<br> delete(YMD )<br><br> x = f->U_anom(iStrt:iLast,:,:) <wbr> ; U_anom_sm<br> printVarSummary(x)<br> printMinMax(x, True)<br><br></div><div> [SNIP] <br><br>;*****************************<wbr>*******************<br>; Apply the Lanczos band pass filter weights to the original anomalies (time,lat,lon)<br>;*****************************<wbr>*******************<br><br> xBPF = wgt_runave_n_Wrap (x, wgt, 0, 0) ; BPF (Band Pass Filtered)<br></div><div>[SNIP]<br><br>==========<br><br></div><div>If you do a <br></div><div> printVarSummary(xBPF)<br></div><div> printMinMax(xBPF,0)<br><br></div><div>You will 'see' the xBPF contains the Lanczos filtered daily values.<br><br>========<br> <br><br></div><div>D<br></div><div><br><br><br><br><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Sep 12, 2016 at 12:42 PM, <span dir="ltr"><<a href="mailto:marilia.gregorio@cptec.inpe.br" target="_blank">marilia.gregorio@cptec.inpe.br</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello, how are you?<br>
<br>
I'm trying to use the mjoclivar_3.ncl script, but would like to generate<br>
a field filtered by lanczos of method for every day, you guys know how I<br>
do it?<br>
<br>
Thank you.<br>
______________________________<wbr>_________________<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" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
</blockquote></div><br></div>