I didn't expect a response from the author himself. Thank you for your time!
I expect that I'm just miss-reading the documentation, but I sure can't see where I'm going wrong.
Right now I'm using values from 8 November 00:00:00 UT in the ephemeris to check my code with. My expectation is that if I feed OE's from the ephemeris and a MJD from a line in the ephemeris into Plante I should get the ra/dec for that line. For reference, the ephemeris line I'm using is:
Date__(UT)__HR:MN:SC.fff Date_________JDUT R.A.___(ICRF/J2000.0)___DEC R.A.__(airls-apparent)__DEC
2011-Nov-08 00:00:00.000 2455873.500000000 * 18 22 52.6104 -19 16 57.681 18 23 34.6940 -19 16 30.326
> What time scale are you using for the date argument?
The date is 55873.0 MJD (2455873.500000000 - 2400000.5) from above.
> What sort of RA,Dec are you comparing the results with?
The ones from the ephemeris line (above).
> Have you got refraction enabled in the JPL predictions?
No, because the way I read the Plante call documentation I adjust for refraction after the Plante call.
The OE's from the ephemeris are:
EPOCH= 2454648.5 ! 2008-Jul-01.00 (CT)
EC= 2.566687011667356E-01 QR= 2.980651303494526E+01 TP= 2448174.216430557892
OM= 1.103600311387948E+02 W = 1.161890224972339E+02 IN= 1.709558287650895E+01
A = 4.009855777865089E+01 MA= 25.130588774518851 ADIST= 5.039060252235653E+01
The values (from above) passed in are:
date: (55873.00000000000000000000)
epoch: (54648.00000000000000000000)
orbinc: ( 17.095582876509 +17:05:44.9835543 (0.298374) )
anode: ( 110.360031138795 +110:21:36.11209966 (1.926146) )
perih: ( 116.189022497234 +116:11:20.48099004 (2.027881) )
aorq: (40.09855777865089265788)
e: (0.25666870116673562441)
aorl: ( 25.130588774519 +25:07:50.11958827 (0.438612) )
dm: 0.000000
The result is:
RA: 275.927687261573 +18:23:42.64494278 (4.815847)
Dec -19.279067405550 -19:16:44.64265998 (-0.336483)
----
Starlink User Support list
For list configuration, including subscribing to and unsubscribing from the list, see
https://www.jiscmail.ac.uk/cgi-bin/webadmin?A0=STARLINK
|