Spectral Analysis

spec_1.ncl: Reads in a time series from a netCDF file, calculates the spectrum and plots the results.

specx_anal: Calculates the spectra of a time series with options on smoothing, demeaning, detrending, and tapering.

dimsizes: calculates the dimension sizes of any array.

spec_2.ncl: Reads in two time series and calculates their spectra, cospectra, coherence, quadrature spectra, and phase and creates a panel plot of all values.

specxy_anal: Calculates the cospectra, coherence, quadrature spectra, and phase as well as the spectra of two time series with options on smoothing, demeaning, detrending, and tapering.

gsn_panel, is the plot interface that panels plots together. More examples of panel plots.

spec_3.ncl: Calculates red noise confidence intervals and creates several plot variations.

specx_ci is the Shea Utility that calculates a red noise confidence interval. This function returns four curves: the calculated spectrum, the "red noise" curve and the curves indicating the upper and lower confidence bounds.

spec_4.ncl: Plots spectra and "red noise" confidence intervals on a log scale and draws the bandwidth.

trYLog = True
trYMinF = 0.10
trYMaxF = 30.0, Changes the Y-axis to log and manually sets the lower and upper limits.

gsn_polyline and gsn_text are used to add the bandwidth line and label to the figure.