Module LOTPS in CMLIBPasses smooth function through two-dimensional scattered data and returns an array of interpolated values on user specified grid. Based on local thin-plate splines.
CALL LOTPS(MODE, NPPR, NPI, XI, YI, FI, NXO, XO, NYO, YO, IWK, NIWK, IJWKU, WK, NWK, FO, KER)| E2b | Interpolation of scattered, non-gridded multivariate data |