15 namespace GeographicLib {
19 void MagneticCircle::Field(
real lon,
bool diffp,
23 CircularEngine::cossin(lon, clam, slam);
24 real M[Geocentric::dim2_];
25 Geocentric::Rotation(_sphi, _cphi, slam, clam, M);
26 real BX0, BY0, BZ0, BX1, BY1, BZ1;
27 _circ0(clam, slam, BX0, BY0, BZ0);
28 _circ1(clam, slam, BX1, BY1, BZ1);
30 BX1 = (BX1 - BX0) / _dt0;
31 BY1 = (BY1 - BY0) / _dt0;
32 BZ1 = (BZ1 - BZ0) / _dt0;
38 Geocentric::Unrotate(M, BX1, BY1, BZ1, Bxt, Byt, Bzt);
43 Geocentric::Unrotate(M, BX0, BY0, BZ0, Bx, By, Bz);