[ncl-talk] Lanczos Filtering

Kiswendsida Hyacinthe GUIGMA karongseba at gmail.com
Tue Dec 5 09:14:08 MST 2017


Dear NCL users,I am very new on NCL and I am trying to filter my rainfall data using Lanczos Filter (bandpass 10-90 day filtering). My question concerns the processing of the weights. I would like to know 
1-is there any relationship between the number of weights to choose and the cutoff frequencies?2-is there any practical criterion to choose the number of weights? I know that this number influences the loss of data.Below is my code:
----------------------------------------------------------------------------------------------------------------------------------------------------------------------;import the file
gpcp=addfile("gpcp_Sahel.nc","r")
rr=gpcp->precip
;filter parameters
ihp=2 ;bandpass
sigma=1.0 ; Lanczos sigma
nWgt    = 81  ; number of weights -> means 40 data lost at the beginning and 40 data lost at the end.                      
fca     = 1./90.            ; start freq ;highpass= 90 days
fcb     = 1./10.                        ; last  freq; lowpass= 10 days
wgt     = filwgts_lanczos (nWgt, ihp, fca, fcb, sigma )
printVarSummary(wgt)
xBPF    = wgt_runave_n_Wrap ( rr, wgt, 0, 0)      ; 10-90 day bp performed across time dimension(0) of rr (meaning the time dimension)
 ---------------------------------------------------------------------------------------------------------------------------------------------------------------------




Cheers,

GUIGMA
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20171205/17200199/attachment.html>


More information about the ncl-talk mailing list