Correction to previous response. I neglected the latitude array subscript.
xLong = ezfftb(ab(:,:,0:nLong-1), 0.0) ; total amplitude (1-3)
xShort = ezfftb(ab(:,:,nLong: ) , 0.0)
>>I have a easy question , I want to get long wave(1-3)
>>and short wave(4-) by spatial Fourier decomposition ,
>>do NCL have a function can do this?
>>x(lat ,lon)---> xwave
>>
>>Thanks for any suggestion!
>
>--
>Yes, of course. One approach is to use the FFT routines: ezfftf, ezfftb
>(Missing data are not allowed.)
>
>See Example 3
>
>--
> nLong = 3
> xLong = x
> xShort = x
>
> ab = ezfftf(x)
>
> xLong = ezfftb(ab(:,0:nLong-1), 0.0) ; total amplitude (1-3)
> xLong@long_name = "Long Waves"
>
> xShort = ezfftb(ab(:,nLong: ), 0.0) ; total amplitude (4-nlon/2)
> xShort@long_name = "Short Waves"
>
>--
>
>Youc could also try the spherical harmonic analysis/synthesis
>family of functions "sha*" and "shs*"
>
>good luck
>D
>
>_______________________________________________
>ncl-talk mailing list
>ncl-talk@ucar.edu
>http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
ncl-talk@ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2b29 : Thu Apr 22 2004 - 07:40:18 MDT