mirror
/
XEphem
mirror of https://github.com/XEphem/XEphem.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

Add quick test showing Pluto jump in year 2247

This commit is contained in:
Brandon Rhodes 2022-01-24 16:08:40 -05:00
parent 06729a304a
commit 5317eae6eb
2 changed files with 37 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
tests/output.expected
View File

@ -13,3 +13,7 @@ RA +5.0000000000000000 -> +4.9986250865903390
Dec -5.0000000000000000 -> +1.2835978054231891
RA +6.0000000000000000 -> +6.0032362075042789
Dec -6.0000000000000000 -> +0.2845854258344847
Pluto discontinuity when Chapront (1995) ends on 2247/10/1
RA = 16:52:53.80
RA = 11:52:36.32

33
tests/test.c
View File

@ -1,6 +1,8 @@
#include "astro.h"
int main() {
char out[1024];
printf("Precession J2000 -> J2015\n");
double ra, pra;
double dec, pdec;
@ -11,5 +13,36 @@ int main() {
printf("RA %+.16f -> %+.16f\n", ra, pra);
printf("Dec %+.16f -> %+.16f\n", dec, pdec);
}
printf("\nPluto discontinuity when Chapront (1995) ends on 2247/10/1\n");
double mjd0, mjd1;
cal_mjd(5, 16, 2247, &mjd0);
cal_mjd(9, 27, 2248, &mjd1);
Now n;
n.n_lat = n.n_lng = n.n_tz = n.n_temp = n.n_pressure = n.n_elev
= n.n_dip = 0.0;
n.n_epoch = J2000;
Obj o;
o.o_type = PLANET;
o.pl_code = PLUTO;
o.pl_moon = X_PLANET;
n.n_mjd = mjd0;
obj_cir(&n, &o);
fs_sexa(out, radhr(o.s_gaera), 2, 360000);
printf("RA = %s\n", out);
n.n_mjd = mjd1;
obj_cir(&n, &o);
fs_sexa(out, radhr(o.s_gaera), 2, 360000);
printf("RA = %s\n", out);
/* From Elwood's email:
5/16/2247 00:00 UTC 4.1.0 RA 16:52:53.90 Dec -10:36:26.5
3.6.7 RA 16:52:53.90 Dec -10:36:26.5
9/27/2248 00:00 UTC 4.1.0 RA 11:52:36.30 Dec 16:13:59.6
3.6.7 RA 16:46:16.01 Dec -11:05:43.7
*/
return 0;
}