NASA Logo
Ocean Color Science Software

ocssw V2022
nutate.c
Go to the documentation of this file.
1 //
2 // Created by alex on 11/25/23.
3 //
4 #include "sun2000.h"
5 #include "nav.h"
6 
7 void nutate(double t, double xls, double gs, double xlm, double omega, double *dpsi, double *eps) {
8  // c Nutation in Longitude
9  *dpsi = -17.1996e0 * sin(omega / radeg) + 0.2062e0 * sin(2.0e0 * omega / radeg) -
10  1.3187e0 * sin(2.0e0 * xls / radeg) + 0.1426e0 * sin(gs / radeg) -
11  0.2274e0 * sin(2.0e0 * xlm / radeg);
12  // c Mean Obliquity of the Ecliptic
13  double epsm = 23.439291e0 - 3.560e-7 * t;
14  // Nutation in Obliquity
15  double deps = 9.2025e0 * cos(omega / radeg) + 0.5736e0 * cos(2.0e0 * xls / radeg);
16  // True Obliquity of the Ecliptic
17  *eps = epsm + deps / 3600.e0;
18  *dpsi = *dpsi / 3600.e0;
19 }
data_t t[NROOTS+1]
Definition: decode_rs.h:77
void nutate(double t, double xls, double gs, double xlm, double omega, double *dpsi, double *eps)
This subroutine computes the nutation in longitude and the obliquity of the ecliptic corrected for nu...
Definition: nutate.c:7
double dpsi
Definition: sun2000.c:6
double eps
Definition: gha2000.c:3
data_t omega[NROOTS+1]
Definition: decode_rs.h:77