LSCLogarithmic_Grav​2T

A function to grid the disturbing potential from residual gravity anomaly data in 3-D space using least squares collocation.

This function can be used to grid residual disturbing potential (SolHa) from residual gravity anomaly data (FA - in mGal) anywhere in 3 dimensional space using the Forsberg 1987 logarithmic covariance function.

See
Forsberg, R. A New Covariance Model for Inertial Gravimetry and Gradiometry.
1987. JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 92, NO. B2, PAGES 1305-1310.

The required input arguments are the known values, X,Y,Z locations of the gravity observations, the residual gravity signal FA, the parameters of the covariance function C0, D and T, a vector of estimated variances N (one for each observation O(X,Y,Z,G)) and the locations Xi,Yi,Zi where the signal is desired to be interpolated to.

The full variance covariance of the interpolated residual disturbing potential is also output.