[ncl-talk] How to get explained variance from spectra
Dennis Shea
shea at ucar.edu
Thu Sep 16 20:39:43 MDT 2021
As notes in the *specx_anal*
<http://www.ncl.ucar.edu/Document/Functions/Built-in/specx_anal.shtml>
documentation:
-====
These spectra have been normalized so that the area under the curve:
(*spcx*(0)+*spcx*(*N*/2-1))*(df/2) + SUM{*spcx*(1:N/2-2)*df}
equals the variance of the detrended series, where df=(1/*N*)=frequency
spacing.
The units are variance/(unit frequency interval).
====
You know the frequencies corresponding to 5-to-10 days.
5-to-10day variance: var_510 = SUM{*spcx*('f5->f10')*df}
fractional_variance = var_510/total_variance
On Thu, Sep 16, 2021 at 2:39 AM Lyndz via ncl-talk <
ncl-talk at mailman.ucar.edu> wrote:
> Dear NCL-experts,
>
> I have a sample data that I attached here and below is the script that I
> used to plot the spectra.
>
> I found significant peaks between 5-10 (days).
>
> I would like to ask how to estimate the fractional variance of
> total rainfall variability explained by this 5-10 days ISO. Im not sure if
> this can be estimated from the spectra. If not, any suggestions on how to
> do this NCL?
>
>
> -Lyndz
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at mailman.ucar.edu
> List instructions, subscriber options, unsubscribe:
> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210916/c7f73dde/attachment.html>
More information about the ncl-talk
mailing list