Cálculo de las posiciones de los nodos ascendentes y descendentes de la Luna a partir de archivos de efemérides

Aquí hay una manera de hacerlo completamente en línea si no desea descargar las bibliotecas CSPICE. Visite https://naif.jpl.nasa.gov/naif/webgeocalc.html -> haga clic en WebGeocalc y haga clic en buscador de posición:

Llena los datos así:

Tenga en cuenta que una latitud planetocéntrica de 0 significa que la latitud de la eclíptica es 0, ya que estamos usando el marco de referencia ECLIPDATE (eclíptica de la fecha).

Presiona calcular:

Los resultados son todas las veces que la luna cruza la eclíptica entre 1970 y 2040. Esto incluye tanto los nodos ascendentes como los nodos descendentes.

Por supuesto, no quiere los tiempos, quiere las posiciones reales, así que haga clic en "Guardar todos los intervalos":

Regrese al menú de cálculo y presione "vector de estado":

Complétalo de la siguiente manera:

Tenga en cuenta que es posible que deba arrastrar su lista de intervalos guardados al cuadro "Lista de intervalos", o puede que se complete automáticamente.

Después de completar, presione calcular:

Y ahí lo tienes: la lista de posiciones lunares para cuando la luna está en sus nodos ascendentes y descendentes.

Observe que la ascensión recta del nodo se mueve algo linealmente:

pero no exactamente.

Si tuviera que trazar la longitud de la eclíptica del nodo, vería un ajuste más lineal debido al ciclo constante entre los meses dracónicos y siderales de la luna.

¡ La respuesta usando SPICE es excelente!

Un atajo es usar la interfaz JPL Horizons en modo ELEMENTOS ORBITAL, descargar una tabla y usar OMo$\Omega$, la longitud del nodo ascendente en grados.

  Symbol meaning:

    JDTDB    Julian Day Number, Barycentric Dynamical Time
      EC     Eccentricity, e
      QR     Periapsis distance, q (km)
      IN     Inclination w.r.t X-Y plane, i (degrees)
      OM     Longitude of Ascending Node, OMEGA, (degrees)
      W      Argument of Perifocus, w (degrees)
      Tp     Time of periapsis (Julian Day Number)
      N      Mean motion, n (degrees/sec)
      MA     Mean anomaly, M (degrees)
      TA     True anomaly, nu (degrees)
      A      Semi-major axis, a (km)
      AD     Apoapsis distance (km)
      PR     Sidereal orbit period (sec)

#!/usr/bin/env python
# -*- coding: utf-8 -*-

import numpy as np
import matplotlib.pyplot as plt

with open('datamoon.txt', 'r') as infile:
    linez = infile.readlines()

lines = [line.split(',') for line in linez]

Omega = np.array([float(line[5]) for line in lines])

JD = np.array([float(line[0]) for line in lines])
JD -= JD[0]

plt.plot(JD/365.2564 + 10, Omega)
plt.xlim(10, 30)
plt.xlabel('years since 2000')
plt.ylabel('geocentric lunar ascending node Ω (deg)')
plt.show()

Los almanaques nativos indios proporcionan la posición de los nodos. El nodo ascendente es Raahu; el descendente es Kethu. https://srirangaminfo.com/vakya-panchangam-srirangam.php