HNBody
Version 1.2.1

A simple custom driver adding the mean lunar quadrupole. More...
Functions  
void  lunar_version (FILE *f, const char *prefix) 
void  lunar_kick (double(*dv)[3], double t, double(*x)[3], const double m[], int nHL, int n, double dt, const hnb_data_t *sys) 
double  lunar_energy (double t, double(*x)[3], double(*v)[3], const double m[], int nNL, const hnb_data_t *sys) 
int  main (int argc, char *argv[]) 
SRCLevel: ISO C
MTLevel: TBD
This custom driver includes the extra kick and energy due to the Sun's attraction on the mean EarthMoon quadrupole, as derived by Quinn, Duncan, & Tremaine (1991). If QTD_DRIVER
is defined, the expansion of the lunar orbit due to tidal dissipation is included and the standard correction factor f = 0.9473 is used. This fully duplicates the Quinn et al. calculation. Otherwise (and by default), the dissipation is neglected (it is useful mainly when EarthMoon rigid body dynamics are being calculated, and it also breaks energy conservation at the rate of 1e13 per Myr), and the correction f = 0.8525 is used. The latter was found by Varadi, Runnegar, & Ghil (2003) to minimize longterm errors in Earth's orbital elements (note that this turns out to be at the expense of shortterm accuracy, however.)