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=NONORTHERNLINE, NOSOUTHERNLINE=NOSOUTHERNLINE
Parameters
alonin
in
required
longitude of the input data
alatin
in
required
latitude of the input data
olonin
in
required
longitude of the output data
olat
in
required
latitude of the output data
weig
out
addr
out
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
NOSOUTHERNLINE
activate if you don't want to take into account the northen/southern line
of the input data when performing the interpolation.
Version history
Version
$Id: compute_fromreg_imoms3_weigaddr.pro 163 2006-08-29 12:59:46Z navarro $
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 longitudes lines have the same latitude and each latitudes 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 on Wed Sep 13 16:32:13 2006.