Kenneth P. Bowman writes:
> > Compute the "interpolation coordinates" from the original grid
> > j = VALUE_LOCATE(y_original, y)
> > yj = j + (y - y_original[j])/(y_original[j+1] - y_original[j])
To tell you the truth, I can't get this to work at all. :-(
IDL> lat = [-87.5, 50, 25, 0, 30, 45, 64, 87.5]
IDL> y = Scale_Vector(findgen(7), -87.5, 87.5)
IDL> j = Value_Locate(lat, y)
IDL> yj = j + (y - lat[j])/(lat[j+1] - lat[j])
% Program caused arithmetic error: Floating illegal operand
IDL> print, yj
0.000000 0.212121 0.424242 3.00000 3.97222 5.70175 -NaN
Cheers,
David
--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.dfanning.com/
Se****e ma de ni thui. ("Perhaps thou speakest truth.")
|
