topcompute_fromirr_bilinear_weigaddr interpolation
Return value
weig, addr: 2D arrays, weig and addr are the weight and addresses used to
perform the interpolation:
dataout = total(weig*datain[addr], 1)
dataout = reform(dataout, jpia, jpja, /over)
Parameters
olonin
in
required
olat
in
required
omsk
in
required
alonin
in
required
alat
in
required
amsk
in
required
weig
addr
Version history
History
June 2006: Sebastien Masson (smasson@lodyc.jussieu.fr)
Known issues
Restrictions
- the input grid must be an "irregular 2D grid", defined as a grid made
of quadrilateral cells which corners positions are defined with olonin and olat
- We supposed the data are located on a sphere, with a periodicity along
the longitude
- to perform the bilinear interpolation within quadrilateral cells, we
first morph the cell into a square cell and then compute the bilinear
interpolation.
- if some corners of the cell are land points, their weight is set to 0
and the weight is redistributed on the remaining "water" corners
- points located out of the southern and northern boundaries or in cells
containing only land points are set the the same value as their closest neighbor l;