NCL Home>
Application examples>
Data Analysis ||
Data files for some examples
Example pages containing:
tips |
resources | 
functions/procedures
NCL: Simple Fourier Analysis of Climate Data
Global atmospheric data are periodic in longitude (0-360) and climatological 
data are periodic in time (here, 12 months). The Fourier Analysis in the 
following examples uses a climatological data set derived from ERA-Interim 
data spanning 1989-2005.
The fourier_info, ezfftf and 
ezfftb can be used to perform variations of 
Fourier Analysis. Using these functions on a variable with longitude
as the rightmost dimension performs spatial analysis. To examine temporal
harmonics, the input series must be 
reordered
so the dimension 'time' is the rightmost dimension.
 
fanal_1.ncl: 
A variable (here, gropotential height at 500 hPa) dimensioned (time,lat,lon)
is examined to determine the amplitude, phase (1st maximum) and percent variance explained. Harmonics 1 and 2 are displayed.
 
fanal_2.ncl: 
Two similar variables TREFHT (temperature at 2m) and TSKIN (temperature
at actual surface) are compared via harmonic analysis.
 
fanal_3.ncl: 
A forward fast Fourier transform (
ezfftf) performs
a 'Fourier Analysis'. Selected coefficents are set to zero to isolate
different waves. A backward fast Fourier transform (
ezfftb)
is used to perform a 'Fourier Synthesis'.
 
fanal_4.ncl: 
Create three simple sine waves (blue, red, green) and combine the waves 
(superposition; black).
Use 
fourier_info on the combined series 
to derive the amplitudes, phases and per-cent
variances explained by each harmonic. Use a polymarker to mark the 
derived locations of the harmonic phase. The (edited) printed output is:
     Variable: finfo
     Type: float
     Total Size: 36 bytes
                 9 values
     Number of Dimensions: 2
     Dimensions and sizes:	[3] x [3]       ; [3] x [nhx]
     Coordinates: 
     
     (0,0)	20         <=== Amplitudes
     (0,1)	10
     (0,2)	 5
     
     (1,0)	72         <=== Phases: location of first maximum
     (1,1)	36
     (1,2)	12
     
     (2,0)	76.19048   <=== % variance
     (2,1)	19.04762
     (2,2)	4.761904
 
fanal_5.ncl: 
Similar to 
fanal_3 and 
fanal_4 except that spherical
harmonic analysis and synthesis are also used.