eotvos(Lat,Long,H,d​s)

A function to calculate the Eotvos effect using the Harlan 1968 formula

A function to calculate the Eotvos effect using the Harlan (1968) formula
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Ref: https://en.wikipedia.org/wiki/E%C3%B6tv%C3%B6s_effect
Ref: Harlan, R.B. Eotvos Corrections for Airborne Gravimetry 1968. J. Geophysics Res.. 14,17.
Relative to the GRS80 reference ellipsoid
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The input values required are vectors
- Latitude (in decimal degrees)
- Longitude (in decimal degrees)
- Height (in metres)
- ds (the time interval between measurements in seconds)
The Output is the Eotvos effect in mGal.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Example usage:
e.g.
Latitude=[-43;-43.001;-43.002];
Longitude=[120;120.001;120.002];
H=[100;100;100];
ds=10;
Eotvoseffect=eotvos(Latitude,Longitude,H,ds)
%% Output
Eotvoseffect =
0
89.951763118312542
0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
N.B. The first and last of the output values are zero. This is because a derivative is taken, using a central difference, to estimate the velocity.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%