Function lupnt::SolidEarthTideDisplacement

Function Documentation

Vec3d lupnt::SolidEarthTideDisplacement(const Vec3d &r_station, const std::vector<std::pair<Vec3d, double>> &tide_bodies, double GM_earth, double earth_radius_m, double h2 = 0.6078, double l2 = 0.0847)

Degree-2 in-phase solid Earth tide station displacement (IERS Conventions form).

The crust displacement is sum_j (GM_j/GM_E)(R_E^4/R_j^3){ h2 r_hat[(3(R_hat_j.r_hat)^2 - / 1)/2] + 3 l2 (R_hat_j.r_hat)[R_hat_j - (R_hat_j.r_hat) r_hat] } summed over the tide-raising bodies (Moon and Sun). Because it is built from frame-invariant dot products of unit vectors, evaluating it with all vectors expressed in an inertial frame yields the displacement directly in that frame.

Parameters:

  • r_station – Geocentric station position [m].

  • tide_bodies – List of (geocentric body position [m], GM [m^3/s^2]) (Moon, Sun).

  • GM_earth – Earth gravitational parameter [m^3/s^2].

  • earth_radius_m – Earth equatorial radius [m].

  • h2 – Degree-2 vertical (radial) Love number.

  • l2 – Degree-2 horizontal Shida number.

Returns:

Station displacement vector [m], same frame as the inputs.