wgrib2: new_grid_interpolation
Introduction
The new_grid_interpolation option selects the type of
interpolation used by the new_grid option.
The possible values are bilinear, bicubic, nearest neighbor, and budget. (Spectral is planned.)
The new_grid_interpolation option must appear before the
new_grid option. If the new_grid_interpolation
option is not used, the interpolation defaults to bilinear.
Some fields may require different types of interpolation such as the soil type should
be found using a nearest neighbor interpolation. (A fractional soil type is meaningless.)
It is common practice to use if and fi options
to control the setting of the the interpolation type as shown by the following command.
wgrib2 IN.grb new_grid_winds earth new_grid_interpolation bilinear \
if ":(VGTYPSOTYP):" new_grid_interpolation neighbor fi \
new_grid latlon 0:360:1 90:181:1 OUT.grb
Spectral interpolation is specialized interpolation scheme. A global field
is transformed into spherical harmonics (user specified truncation), and
the grid point values are determined from the spherical harmonics. This
interpolation is unlike other interpolation scheme. For example, the
nearest neighbor is based on the 1 grid point, the bilinear internpolation
is based on 4 grid points. The spectral interpolation is based on all
the grid points. So the interpolation scheme will act as noise reduction
when projected to a fewer spherical harmonics than the number of grid points.
The limiting zonal wave number is specified and should not be prime as it makes
the Fast Fourier Transform into a slow
Fourier Transform. The limiting zonal wave number should be a product
of many 2s, 3s, 5s and other small primes. The limiting zonal wave number
may need to be compatible with the Fast Fourier Transform (FFT) code used.
Usage
new_grid_interpolation X
X = bilinear, bicubic, neighbor, budget, spectral(trun)(num)
spectra* is alpha
bilinear: linear interpolation in X and then Y.
neighbor: Values from nearest grid point.
budget: Preserves largescale means. Good for precipitation.
spectral: Converts a global field into spectral harmonics
and derives values from the spectral harmonics.
trun = 't', 'T', 'r' or 'R' for triangular or rhomboidal trunction
(num) = max zonal wave number
The input has to be a global field and contain valid values for all grid points.
If undefined values are found, bilinear interpolation is used.
See also: new_grid,
new_grid_winds
new_grid_ipopt
if
fi
