21 real slam, clam, M[Geocentric::dim2_];
23 real Wres = W(slam, clam, gx, gy, gz);
24 Geocentric::Rotation(_sphi, _cphi, slam, clam, M);
25 Geocentric::Unrotate(M, gx, gy, gz, gx, gy, gz);
31 real slam, clam, M[Geocentric::dim2_];
33 real Tres = InternalT(slam, clam, deltax, deltay, deltaz,
true,
true);
34 Geocentric::Rotation(_sphi, _cphi, slam, clam, M);
35 Geocentric::Unrotate(M, deltax, deltay, deltaz, deltax, deltay, deltaz);
40 if ((_caps & GEOID_HEIGHT) != GEOID_HEIGHT)
42 real slam, clam, dummy;
44 real T = InternalT(slam, clam, dummy, dummy, dummy,
false,
false);
45 real correction = _corrmult * _correction(slam, clam);
46 return T/_gamma0 + correction;
50 real& Dg01, real& xi, real& eta)
52 if ((_caps & SPHERICAL_ANOMALY) != SPHERICAL_ANOMALY) {
59 deltax, deltay, deltaz,
60 T = InternalT(slam, clam, deltax, deltay, deltaz,
true,
false);
62 real MC[Geocentric::dim2_];
63 Geocentric::Rotation(_spsi, _cpsi, slam, clam, MC);
64 Geocentric::Unrotate(MC, deltax, deltay, deltaz, deltax, deltay, deltaz);
66 Dg01 = - deltaz - 2 * T * _invR;
73 real Wres = V(slam, clam, gX, gY, gZ) + _frot * _Px / 2;
82 if ((_caps & GRAVITY) != GRAVITY) {
87 Vres = _gravitational(slam, clam, GX, GY, GZ),
88 f = _GMmodel / _amodel;
98 bool gradp,
bool correct)
const {
100 if ((_caps & DISTURBANCE) != DISTURBANCE) {
105 if ((_caps & DISTURBING_POTENTIAL) != DISTURBING_POTENTIAL)
111 ? _disturbing(slam, clam, deltaX, deltaY, deltaZ)
112 : _disturbing(slam, clam));
113 T = (T / _amodel - (correct ? _dzonal0 : 0) * _invR) * _GMmodel;
115 real f = _GMmodel / _amodel;
120 real r3 = _GMmodel * _dzonal0 * _invR * _invR * _invR;
121 deltaX += _Px * clam * r3;
122 deltaY += _Px * slam * r3;
GeographicLib::Math::real real
static void sincosd(T x, T &sinx, T &cosx)
void SphericalAnomaly(real lon, real &Dg01, real &xi, real &eta) const
Math::real Disturbance(real lon, real &deltax, real &deltay, real &deltaz) const
Math::real Gravity(real lon, real &gx, real &gy, real &gz) const
Namespace for GeographicLib.
Header for GeographicLib::Geocentric class.
Math::real GeoidHeight(real lon) const
Header for GeographicLib::GravityCircle class.