NCL Home > Documentation > Functions > WRF, Interpolation


Interpolates to a specified pressure/height level.


	function wrf_interp_3d_z (
		v3d     : numeric,  
		vert    : numeric,  
		loc [1] : numeric   

	return_val  :  float or double



Input array to interpolate. This array must have at least 3 dimensions, whose rightmost dimensions are nz x ny x nx.


An array of the same dimensionality as v3d. Array must be either pressure (hPa) or height (m), depending on which vertical coordinate is used for interpolation. The units must be the same as loc's units.


A scalar. Either pressure (hPa) or height (m), depending on which vertical coordinate is used for interpolation. Data in v3d will be interpolated to this level. If loc is pressure (height), then vert must be pressure (height), and must be the same units.

Return value

The return array will be the same dimensionality as v3d and vert, except with no nz dimension.


This function interpolates terrain following data to either a given pressure or height level. If you have WRF data, then we recommend that you use the wrf_user_intrp3d function instead (see example 2 below).

If v3d has attributes "description" or "units", they will be returned with the output variable. The "description" attribute will be set to " " if v3d doesn't contain this attribute.

As of version 5.1.0, the return variable will contain all but the rightmost third named dimensions (if any) of v3d. Otherwise, the rightmost two dimensions will be named "south_north" and "west_east".

Note: in NCL V6.0.0, the default missing value was changed from -999999 to the default missing for an NCL float or double (9.96921e+36 or 9.969209968386869e+36, depending on output type).

See Also

wrf_interp_1d, wrf_interp_2d_xy, wrf_user_intrp3d, wrf_user_getvar

See the full list of WRF functions.


The function wrf_user_intrp3d, (available in the $NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl script) can also be used to interpolate in height/pressure.

Example 1

    a = addfile("","r")
    tk  = wrf_user_getvar(a,"tk",0)      ; temperature in Kelvin

  ; Interpolate tk to 850hPa
    pf = pf * 0.01                ; Convert to hPa
    var2d = wrf_interp_3d_z(tk,pf,850.) 
Example 2

  a = addfile("","r")

  time = 1
  tc = wrf_user_getvar(a,"tc",time)   ; T in deg C  
  z  = wrf_user_getvar(a,"z",time)    ; z on mass points

  ; Horizontally interpolate to height coordinates ("z")
  height = 2000.   ; 2km 
  tc_plane = wrf_user_intrp3d(tc,z,"h",height,0.,False)