fluxEddy

Calculates time averaged eddy flux quantities.

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.

Synopsis

    function fluxEddy(
        x  : float,
        y  : float
    )

Arguments

x
A float array with one or more dimensions. The rightmost dimension must be "time". Missing values should be indicated by x@_FillValue. If x@_FillValue is not set, then the NCL default (-999) will be assumed.
y
A float array with one or more dimensions. It must have the same dimensions as x. Missing values should be indicated by y@_FillValue. If y@_FillValue is not set, then the NCL default (-999) will be assumed.

Description

fluxEddy returns a float array of the same dimensionality as its input arrays with the last dimension omitted. The resulting quantities are computed in one pass. The approach used is let x=X+x' and y=Y+y' where X and Y are time averages, then x'y' = ave{x*y} - X*Y. Missing values are allowed.

Example 1

Let x and y be dimensioned nlat x mlon x ktimes (nlat=64, mlon=128, ktimes=90): 
    xpyp = fluxEddy (x,y)   ; xpyp=x'y' will be an nlat x mlon array

Example 2

Let u and v be dimensioned ntimes x klev x nlat x mlon with named dimensions "time" , "lev" , "lat" , "lon". Then:

    upvp = fluxEddy (u(lev|:,lat|:,lon|:,time|:), v(lev|:,lat|:,lon|:,time|:))
upvp will be a 3-dimensional array dimensioned klev x nlat x mlon.

Example 3

To calculate the eddy variance of a quantity z. (I.e. ave{z'z'}):

    zpzp = fluxEddy (z,z)
or

    zpzp = fluxEddy (z(lev|:,lat|:,lon|:,time|:), z(lev|:,lat|:,lon|:,time|:))
depending upon the dimension order.

Reference Manual Control Panel

NG4.1 Home, Index, Examples, Glossary, Feedback, Ref Contents, Ref WhereAmI?

$Revision: 1.3 $ $Date: 1999/03/22 20:49:09 $