NAME: compute_fromreg_bilinear_weigaddr
PURPOSE: compute the weight and address neede to interpolate data from a
"regular grid" to any grid using the bilinear method
CATEGORY:interpolation
CALLING SEQUENCE:
compute_fromreg_bilinear_weigaddr, alon, alat, olon, olat, weig, addr
INPUTS:
lonin and latin: longitude/latitude of the input data
lonout and latout: longitude/latitude of the output data
KEYWORD PARAMETERS:
/NONORTHERNLINE and /NOSOUTHERNLINE: activate if you don't whant to take into
account the northen/southern line of the input data when perfoming the
interpolation.
OUTPUTS:
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, jpio, jpjo, /over)
COMMON BLOCKS: none
SIDE EFFECTS: ?
RESTRICTIONS:
- the input grid must be a "regular grid", defined as a grid for which each
lontitudes 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 out of the southern and northern boundaries are interpolated
using a linear interpolation only along the longitudinal direction.
EXAMPLE:
MODIFICATION HISTORY: