
exp_tapershC
Performs tapering (filtering) of the spherical harmonic coefficients.
Prototype
function exp_tapershC ( ab : numeric, N : numeric, r : integer ) return_val [dimsizes(ab)] : float or double
Arguments
abspherical harmonic coefficients returned to the user by the NCL functions shagC/shaeC.
Nmode at which the taper weight will equal exp(-1) [=0.3678795 ]. (N < maximum-wave-number possible)
rpower of exponent. This determines the rate at which the coefficients decrease (see below).
Return value
The results are returned in an array that is of size (dimsizes(ab)).
Description
exp_tapershC performs tapering on the spherical harmonic coefficients.
The effect is to perform (spatial) isotropic smoothing by reducing the amplitude of the coefficients at higher modes. Most frequently, this is done for graphical purposes. The formula used is:
S(n) = EXP{-[ n(n+1)/N(N+1)]^r }n = total wavenumber, N,r defined above. This is equation (9) in the following reference:
Reference:
Spatial Smoothing on the Sphereexp_tapershC is the function version of exp_tapersh.
P. D. Sardeshmukh and B. J. Hoskins
Monthly Weather Review, December 1984, pp 2424-2529.
See Also
exp_tapersh_wgts, exp_tapersh, shagC, shaeC
Examples
Example 1
a and b are calculated from a function (say shagC/shaeC). g is a data array dimensioned ntim x nlat x nlon.
ab = shaeC(g) T = 62 ab = tri_trunC (ab,T) ;triangularly truncate the spherical harmonic coefficients ab = exp_tapershC (ab, 45, 2) ;taper the spherical harmonic coefficients g_mod = shseC(ab,nlon) ;perform the spherical harmonic synthesis on the truncated coefficients, store results in h