~georg-zotti/stellarium/gz_AtmosphereTweaks

                Planet("art", 0, 0, Vec3f(0,0,0), 0, 0, "", "", "", Q_NULLPTR, Q_NULLPTR, Q_NULLPTR, false, true, false, true, "artificial"),

                dest(Q_NULLPTR), orig_name(orig->getEnglishName()), orig_name_i18n(orig->getNameI18n())

// Since V0.14, we follow Meeus, Astr. Alg. 2nd ed, Ch.11., but used offset rho in a wrong way. (offset angle phi in distance rho.)

        if (getHomePlanet()->getRadius()==0.0) // the transitional ArtificialPlanet or SpaceShipObserver has this

        const double a=getHomePlanet()->getRadius();

        const double bByA = getHomePlanet()->getOneMinusOblateness(); // b/a;

        const double latRad=currentLocation.latitude*(M_PI/180.0);

        const double u = atan( bByA * tan(latRad));

        const double altFix = currentLocation.altitude/(1000.0*AU*a);

 

        const double rhoSinPhiPrime= bByA * sin(u) + altFix*sin(latRad);

        //double rhoCosPhiPrime= bByA * cos(u) + altFix*cos(latRad); // WARNING! bByA is not in the book!!! THIS IS A TEST!

        const double rhoCosPhiPrime=        cos(u) + altFix*cos(latRad);

 

        const double rho = sqrt(rhoSinPhiPrime*rhoSinPhiPrime+rhoCosPhiPrime*rhoCosPhiPrime);

// Used to approximate solution with assuming a spherical planet.

// Since V0.14, following Meeus, Astr. Alg. 2nd ed, Ch.11.

// Since V0.16, we can produce the usual offset values plus geocentric latitude phi'.

Vec3d StelObserver::getTopographicOffsetFromCenter(void) const

{

        if (getHomePlanet()->getRadius()==0.0) // the transitional ArtificialPlanet or SpaceShipObserver has this

                return currentLocation.altitude/(1000*AU);

 

        const double a=getHomePlanet()->getRadius();

        const double bByA = getHomePlanet()->getOneMinusOblateness(); // b/a;

 

        if (fabs(currentLocation.latitude)>=89.9) // avoid tan(90) issues.

                return a * bByA;

 

        const double latRad=currentLocation.latitude*(M_PI/180.0);

        const double u = atan( bByA * tan(latRad));

        // qDebug() << "getDistanceFromCenter: a=" << a*AU << "b/a=" << bByA << "b=" << bByA*a *AU  << "latRad=" << latRad << "u=" << u;

        Q_ASSERT(fabs(u)<= fabs(latRad));

        const double altFix = currentLocation.altitude/(1000.0*AU*a);

 

        const double rhoSinPhiPrime= bByA * sin(u) + altFix*sin(latRad);

        const double rhoCosPhiPrime=        cos(u) + altFix*cos(latRad);

 

        const double rho = sqrt(rhoSinPhiPrime*rhoSinPhiPrime+rhoCosPhiPrime*rhoCosPhiPrime);

        //return rho*a;

        double phiPrime=asin(rhoSinPhiPrime/rho);

        return Vec3d(rhoCosPhiPrime*a, rhoSinPhiPrime*a, phiPrime);

}