Calculates cross spectra quantities of a series.
Note: This function is only available in version 4.1.1 of NCL. If your site is licensed for version 4.1, then you can get version 4.1.1 for free. To get version 4.1.1 of NCAR Graphics software, please contact your site representative. If you don't know who your site representative is, then send email to ncarginf@ucar.edu or call (303) 497-1201.
function specxy_anal( x : float, y : float, iopt : integer, jave : integer, pct : float )
jave=0 : do no smoothing spcx contains raw spectra estimates (periodogram) jave>0 : average jave periodogram estimates together utilizing modified daniell smoothing (good stability but may lead to large bias ). all weights are 1/jave except wgt(1) and wgt(jave) which are 1/(2*jave). this is the recommended option. It is this number which has the most impact on the degrees of freedom.
These spectra have been normalized so that SUM{spcx(n)*df}= variance of the detrended series where df=frequency spacing. The units are variance/(unit frequency interval).
iopt = 1 ; remove least squares linear trends from each ; series prior to tapering and computing spectra. jave = 7 ; Average 7 periodogram estimates using modified Daniel pct = 0.1 ; taper 10% of the data dof = specxy_anal (x,y,iopt,jave,pct)
NG4.1 Home, Index, Examples, Glossary, Feedback, Ref Contents, Ref WhereAmI?