Question on the low pass filter?

Dear NCL user:

      Hi, I want to do a low pass filter (filtering signal with period less than 10a) for a 100a time series.

After see some tips on the forum , I still have some qestions.

Using "ezfftf" and "ezfftb" :
 
      x = random_uniform(-1,1,100) ; using a random series as example
      fc = ezfftf (x)
      fc (:,5:49) = 0.0 ; Is it right for filtering the period longer than 10a?
      xband = ezzftb(fc) ; the result

Using filwgts_lancos :
 
    I want to know how to define the nwt and fca to get a similar low pass filter as listed above ?

Thanks!

Lin
      
       
   
 
    

LEO_ARIES
2006-12-26

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Dec 25 2006 - 11:17:12 MST