def calc_geocentric_alignments(planet_name, d):

    sun = calc_orbital_elements('Sun', d)
    planet = calc_orbital_elements(planet_name, d)

    lonsun = sun.get('v') + sun.get('w')

    xs = sun.get('r') * math.cos(lonsun)
    ys = sun.get('r') * math.sin(lonsun)

    xh = planet.get('xh')
    yh = planet.get('yh')

    xg = xh + xs
    yg = yh + ys

    helio_degree = math.degrees(math.atan2( xh, yh ))
    geo_degree = math.degrees(math.atan2( xg, yg ))


    helio_degree =  90 - helio_degree
    geo_degree = 90 + geo_degree

    if helio_degree < 0:
        helio_degree = helio_degree + 360
    if geo_degree < 0:
        geo_degree = geo_degree + 360

    return geo_degree
dataset/ETHPy150Open lhartikk/AstroBuild/astro_build.py/calc_geocentric_alignments/operator `-` -> `+`