NCL Home > Documentation > Functions > Interpolation

linint1

Interpolates from one series to another using piecewise linear interpolation across the rightmost dimension.

Prototype

	function linint1 (
		xi           : numeric,  
		fi           : numeric,  
		fiCyclic [1] : logical,  
		xo       [*] : numeric,  
		foOption     : integer   
	)

	return_val  :  float or double

Arguments

xi

An array that specifies the X coordinates of the fi array. It must be strictly monotonically increasing or decreasing, and can be unequally spaced. If xi is multi-dimensional, then its dimensions must be the same size as fi's dimensions. If it is one-dimensional, its length (call it nxi) must be the same as the rightmost (fastest varying) dimension of fi.

fi

An array of one or more dimensions. The rightmost dimension (nxi) is the dimension to be interpolated. If missing values are present, the attribute fi@_FillValue must be set appropriately.

fiCyclic

An option to indicate whether the rightmost dimension of fi is cyclic.

This should be set to True only if you have global data, but your longitude values don't quite wrap all the way around the globe. For example, if your longitude values go from, say, -179.75 to 179.75, or 0.5 to 359.5, then you would set this to True.

xo

A one-dimensional array that specifies the X coordinates of the return array. It must be strictly monotonically increasing or decreasing and may be unequally spaced.

foOption

Reserved for future use. It is currently not used, but set it to 0.

Return value

The returned value will have the same dimensions as fi, except for the rightmost dimension which will have the same dimension size as the length of xo. The return type will be double if fi is double, and float otherwise.

Description

The linint1 function uses piecewise linear interpolation to interpolate from one series to another. The series may be cyclic in the X direction.

If missing values are present, then linint1 will perform the piecewise linear interpolation at all points possible, but will return missing values at coordinates which could not be used.

If any of the the output coordinates xo are outside those of the input coordinates xi, the fo values at those coordinates will be set to missing (i.e. no extrapolation is performed).

Use linint1_n if the dimension to do the interpolation across is not the rightmost dimension. This function can be significantly faster than linint1.

Use the linint1_Wrap function if metadata retention is desired. The interface is identical.

See Also

linint1_n, linint1_Wrap, linint1_n_Wrap, linint2, linint2_points, linmsg, ESMF_regrid

Examples

Example 1

Assume fi is a 1D array, xi is a 1D array with values from 30 to 80 (they don't have to be equally-spaced), and that the rightmost dimension of fi is not to be treated as cyclic. Further assume that the output grid, xo, also has values from 30 to 80. Then, to interpolate fi to the grid represented by xo:

  fo = linint1 (xi,fi, False, xo, 0)

  ; Use linint1_Wrap if metadata retention is desired
  ; fo = linint1_Wrap (xi,fi, False, xo, 0)

fo will be 1D and be the same size as xo.

Example 2

Assume fi is dimensioned ntim x nlvl x nlat x mlon (ntim=50, nlvl=30, nlat=64, mlon=128), and that the rightmost dimension is to be treated as cyclic (the user should not add a cyclic point for the rightmost dimension). All times, levels, and latitudes will be interpolated and returned in a new array fo, dimensioned ntim x nlvl x nlat x 144:

  lon = (0., 2.8125, .... , 357,0125)

  LON = (0., 2.5, ... , 357.5)    ; length 144

  fo = linint1 (lonfi, fi, True, LON, 0)

  ; Use linint1_Wrap if metadata retention is desired
  ; fo = linint1_Wrap (lonfi, fi, True, LON, 0)

Example 3

Assume xi is dimensioned ntim x nlvl x nlat x mlon (ntim=100, nlvl=30, nlat=64, mlon=128) and has named dimensions "time", "lev", "lat", "lon" and coordinate variables of the same name. Further, assume the values of time range from 15 to 500.

To create new functional values at arbitrarily specified times, the following approach could be used:

  tNew = (/15., 15.5,16.8,19.0, ...488.23 /)  ; new times
  time = xi&time   ; for clarity

  ; Use linint1_n to avoid reordering:
  xo   = linint1_n (time, x, False, tNew, 0, 0)

  ; Use linint1_n_Wrap if metadata retention is desired
  ;  xo   = linint1_n_Wrap (time, x, False, tNew, 0, 0) 

The function will interpolate all levels, latitudes and longitudes to the user-specified times and return in a new array xo. If NTIM = dimsizes(tNew) (number of new time steps), then the returned array, xo, will be of size nlvl x nlat x mlon x NTIM.