Skip Navigation Links www.nws.noaa.gov 
NOAA logo - Click to go to the NOAA home page National Weather Service   NWS logo - Click to go to the NWS home page
Climate Prediction Center
 
 

 
About Us
   Our Mission
   Who We Are

Contact Us
   CPC Information
   CPC Web Team

 
HOME > Monitoring_and_Data > Oceanic and Atmospheric Data > Reanalysis: Atmospheric Data > wgrib2-new_grid_interpolation
 

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 ":(VGTYP|SOTYP):" -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 large-scale 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

NOAA/ National Weather Service
National Centers for Environmental Prediction
Climate Prediction Center
5830 University Research Court
College Park, Maryland 20740
Climate Prediction Center Web Team
Page last modified: July 2, 2010
Disclaimer Privacy Policy