Interpolation/
compute_fromreg_imoms3_weigaddr.pro
compute the weight and address need to interpolate data from a
"regular grid" to any grid using the imoms3 method
compute_fromreg_imoms3_weigaddr
Interpolation
compute_fromreg_imoms3_weigaddr, alonin, alatin, olonin, olat, weig, addr, NONORTHERNLINE=scalar 0 or 1, NOSOUTHERNLINE=scalar 0 or 1
Parameters
alonin
in
required
type: 2d array
longitude of the input data
alatin
in
required
type: 2d array
latitude of the input data
olonin
in
required
type: 2d array
longitude of the output data
olat
in
required
type: 2d array
latitude of the output data
weig
out
type: 2d array
(see ADDR)
addr
out
type: 2d array
2D arrays, weig and addr are the weight and addresses used to
perform the interpolation:
dataout = total(weig*datain[addr], 1)
dataout = reform(dataout, jpio, jpjo, /over)
Keywords
NONORTHERNLINE
type: scalar 0 or 1
default: 0
put 1 if you don't want to take into
account the northern line of the input data when performing the interpolation.
NOSOUTHERNLINE
type: scalar 0 or 1
default: 0
put 1 if you don't want to take into
account the southern line of the input data when performing the interpolation.
Version history
Version
$Id: compute_fromreg_imoms3_weigaddr.pro 372 2008-08-08 12:31:53Z pinsard $
History
November 2005: Sebastien Masson (smasson@lodyc.jussieu.fr)
March 2006: works for rectangular grids
Known issues
Restrictions
- the input grid must be a "regular/rectangular grid", defined as a grid for
which each longitude lines have the same latitude and each latitude columns
have the same longitude.
- We supposed the data are located on a sphere, with a periodicity along
the longitude.
- points located between the first/last 2 lines are interpolated
using a imoms3 interpolation along the longitudinal direction and linear
interpolation along the latitudinal direction
- points located out of the southern and northern boundaries are interpolated
using a imoms3 interpolation only along the longitudinal direction.
Produced by IDLdoc 2.0.
