def __calc_equatorialCoordinates(self):
''' calculates position in equatorial coords from ecliptic '''
import math
Jcty = self._TSMgr.time_get_julianCenturiesJ2000()
Lam_deg = self._moonCoordinates["ecl_Lon"] # Lambda
Bet_deg = self._moonCoordinates["ecl_Lat"] # Beta
Eps_deg = MH.eps_deg(Jcty) # Epsilon
sinLam = MH.sin_deg(Lam_deg)
cosLam = MH.cos_deg(Lam_deg)
sinBet = MH.sin_deg(Bet_deg)
cosBet = MH.cos_deg(Bet_deg)
sinEps = MH.sin_deg(Eps_deg)
cosEps = MH.cos_deg(Eps_deg)
# [o. Montenbruck] "Grundlagen der Ephemeridenrechnung" S.14
# cos(del)*cos(alp) = cos(bet)*cos(lam)
# cos(del)*sin(alp) = cos(eps)*cos(bet)*sin(lam) - sin(eps)*sin(bet)
# sin(del) = sin(eps)*cos(bet)*sin(lam) + cos(eps)*sin(bet)
Del_rad = math.asin( sinEps*cosBet*sinLam + cosEps*sinBet )
Del_deg = math.degrees(Del_rad)
cosDel = math.cos(Del_rad)
cosAlp = (cosBet*cosLam) / cosDel
sinAlp = (cosEps*cosBet*sinLam - sinEps*sinBet) / cosDel
Alp_deg = MH.atan2_deg(sinAlp,cosAlp)%360
self._moonCoordinates["equa_Lon"] = Alp_deg # Alpha / RA
self._moonCoordinates["equa_Lat"] = Del_deg # Delta / DE
def __init__(self, TSMgrObj: TSMgr.TimeSpaceMgr, planet: str):
''' Create a planet by its name,
with reference to a TimeSpaceMgr object
(preferably one for all planets). '''
#print("PLANET: "+str(planet))
# default values
self._successPlanet = False
self._planetName = None
self._TSMgr = None
if (planet.lower() in Planetlist.Planetlist
and planet.lower() != "earth"):
#wenn planet sich in der liste bekannter planeten befindet
self._planetName = planet.lower()
self._TSMgr = TSMgrObj # reference to the global TSMgr
Jcty = self._TSMgr.time_get_julianCenturiesJ2000()
self._keplerPlanet = Planetlist.getCurrentKeplerElem(
self._planetName, Jcty) #get kepler elements
self._keplerEarth = Planetlist.getCurrentKeplerElem(
"earth", Jcty) #get kepler elements
# astronomical almanach 2018
# schiefe der eklitptic epsilon = epsJ2000 + epsDot * T
# T = julian centuries since J2000
# epsJ2000 = 84381.406 arcsec
# epsDot = -46.836769 arcsec per julian century
# self._eps = (84381.406 -
# 46.836769*self._TSMgr.time_get_julianCenturiesJ2000()
# )/3600
self._eps = MH.eps_deg(Jcty)
# astronomic unit in km
self._AU = 149597870.7
self._successPlanet = True