Regridding

There are numerous regridding functions available in NCL. In NCL version 6.1.0, there will be new regridding capabilities available via the use of software from the Earth System Modeling Framework (ESMF). This new software will do regridding on rectilinear, curvilinear, and unstructured grids, using bilinear, patch, or conservative interpolation.

Of the older regridding routines, some are unique (eg: use of spherical harmonics); conservative remapping (eg, area_conserve_remap and area_conserve_remap_Wrap). Conventional bilinear interpolation is available (eg: linint2 and linint2_Wrap). Some allow regridding from rectilinear grids to curvilinear grids, (eg: rgrid2rcm and rgrid2rcm_Wrap), and curvilinear grids to rectilinear grids (eg: rcm2rgrid and rcm2rgrid_Wrap).

Other functions exist to "fill" existing grids via extrapolation or a Poisson based algorithm. The Grid_Fill and Vertical Interpolation Application pages provide examples.

regrid_1.ncl: An example of using linint2, which interpolates from one grid to another using bilinear interpolation.

regrid_2.ncl: An example of using g2gsh, which interpolates from one gaussian grid to another using spherical harmonics.

Example: the actual interpolation is conducted using g2gsh_Wrap, a wrapper function that will assign all the appropriate meta data, including the gaussian latitudes, to resulting output.

regrid_3.ncl: An example of using g2fsh, which interpolates from one gaussian grid to a fixed grid using spherical harmonics.

regrid_4.ncl: An example of using f2fsh, which interpolates from one fixed grid to another using spherical harmonics.

regrid_5.ncl: An example of using f2gsh, which interpolates from a fixed grid to a gaussian grid using spherical harmonics.

Even though this example uses f2gsh_Wrap the coordinate variables could be created manually. See example one for the creation of the longitude array. The resulting gaussian latitude array can be created using latGau.

regrid_6.ncl: An example of using area_conserve_remap_Wrap, to perform an interpolation from a high resolution fixed (regular) grid to a lower resolution fixed grid. The interpolation is globally conservative.

regrid_7.ncl: An example of using area_conserve_remap_Wrap, to perform an interpolation from a high resolution fixed (regular) grid to a lower resolution Gaussian grid. The interpolation is globally conservative.

regrid_8.ncl: An example of using area_conserve_remap_Wrap, to perform an interpolation from a high resolution Gaussian grid to a lower resolution Gaussian grid. The interpolation is globally conservative.

regrid_9.ncl: An example of using area_conserve_remap_Wrap, to perform an interpolation from a high resolution Gaussian grid to a lower resolution fixed grid. The interpolation is globally conservative.

regrid_10.ncl: An example of using area_conserve_remap_Wrap, to perform an interpolation from a high resolution fixed (regular) grid with limited latitudinal extent to a lower resolution fixed grid with approximately the same latitudinal extent. The interpolation is globally conservative.

regrid_11.ncl: An example of using area_conserve_remap_Wrap, to perform an interpolation from a high resolution fixed (regular) grid with limited latitudinal extent to a lower resolution gaussian grid with approximately the same latitudinal extent. The interpolation is not conservative.

regrid_12.ncl: Another simple example of using area_conserve_remap_Wrap to perform an interpolation from a high resolution regular grid to a lower resolution grid. The variable being regridded is surface topography [elevation].

regrid_13.ncl: Another example of using area_conserve_remap_Wrap to perform an interpolation from a high resolution regular grid to a lower resolution grid. The variable being regridded is surface topography [elevation]. It is demonstrated how one might create a set of lat/lon values that outline Tibet. The user specified variable "zcrit" can be changed as needed. The value "zcrit=1500" corresponds to approximately 850 hPa. "zcrit" values of 2500, 3500 and 4500 would correspond to approxiamtely 750, 650 and 575 hPa, respectively.