ksearthshadow.cpp

/*
    SPDX-FileCopyrightText: 2018 Valentin Boettcher <valentin@boettcher.cf>
 
    SPDX-License-Identifier: GPL-2.0-or-later
*/
#include "ksearthshadow.h"
#include "kssun.h"
#include "ksmoon.h"
 
KSEarthShadow::KSEarthShadow(const KSMoon *moon, const KSSun *sun, const KSPlanet * earth)
    : KSPlanetBase (i18n("Earth Shadow")), m_sun { sun }, m_moon { moon }, m_earth {earth}
{
}
 
bool KSEarthShadow::shouldUpdate()
{
    return (m_moon->illum() > 0.8);
}
 
bool KSEarthShadow::isInEclipse()
{
    double dist = angularDistanceTo(m_moon).Degrees() * 60;
    return (dist - m_moon->angSize() / 2) <= m_penumbra_ang;
}
 
KSEarthShadow::ECLIPSE_TYPE KSEarthShadow::getEclipseType()
{
    double half_a_moon = m_moon->angSize() / 2;
    double dist = angularDistanceTo(m_moon).Degrees() * 60; // arcminutes
    if (dist <= (m_penumbra_ang - half_a_moon))
    {
        if(dist <= (m_umbra_ang - half_a_moon))
            return FULL_UMBRA;
        else
            return FULL_PENUMBRA;
    }
    else if ((dist - half_a_moon) <= m_penumbra_ang)
    {
        return PARTIAL;
    }
 
    return NONE;
}
 
bool KSEarthShadow::findGeocentricPosition(const KSNumbers *, const KSPlanetBase *)
{
    updateCoords();
    return true; //TODO: not nice!
}
 
 
void KSEarthShadow::updateCoords(const KSNumbers *, bool, const CachingDms *, const CachingDms *, bool )
{
    updateCoords();
}
 
//TODO: Abort if Null
void KSEarthShadow::updateCoords()
{
    // flip the sun around to get the 'shadow coordinates'
    dms t_ra(m_sun->ra().Degrees() + 180);
    t_ra.reduceToRange(dms::ZERO_TO_2PI);
    dms t_dec(-1 * (m_sun->dec().Degrees()));
 
    set(t_ra, t_dec);
    Rearth = m_moon->rearth();
}
 
 
//NOTE: This can easily be generalized to any three bodies.
void KSEarthShadow::calculateShadowRadius()
{
    double d_sun = m_sun->rearth() * AU_KM;
    double d_moon = m_moon->rearth() * AU_KM;
    double r_sun = m_sun->physicalSize() / 2;
    double r_earth = m_earth->physicalSize() / 2;
 
    double umbraRad = 1.01 * r_earth + (r_earth - r_sun) / d_sun * d_moon;
    double penumbraRad = 1.01 * r_earth + (r_sun + r_earth) / d_sun * d_moon;
 
    m_umbra_ang = asin(umbraRad / d_moon) * 60. * 180. / dms::PI;
    m_penumbra_ang = asin(penumbraRad / d_moon) * 60. * 180. / dms::PI;
 
    return;
}