I'm using https://github.com/Bill-Gray/jpl_eph for a personal project to read JPL DE files, namely DE430 and the more recent DE440. This generally works great and when tested against a reference implementation the output shows the correct vectors for location and velocity.
I am puzzled as to what time scale I should use for the DE files or more specifically, how I can convert the time scale that I have available (UTC) to the required TDB.
I read SOFA, explanations in the docs of Skyfield and https://gssc.esa.int/navipedia/index.php/Transformations_between_Time_Systems and have built a beginner's understanding of the time scales and their differences.
Question 1: My understanding is that DE430 and later version require TDB as input variable. Is my current working assumption correct that I can therefore take a time in UTC, add the appropriate number of leap-seconds that have been inserted since 1972 to arrive at TAI, add 32.184s to arrive at TT and use the approximate formula shown here to arrive at an approximation of TDB?
I would then convert this to Julian days and use as input to DE430?
Question 2: How do I do this for dates before 1972? I should mention that like any beginner, I have stumbled over $${\Delta}T = TT - UT$$ and historical records of it, but this seems to be the difference between UT (or UT1) rather than UTC. I think this is where I have the biggest confusion. Can I use records of ${\Delta}T$ available at the IERS to arrive at the correct TT (and use that to calculate TDB)? If so, how?
Question 3: Setting JPL DE aside for a moment, I've also came across other ephemeris like Brown's theory of the motion of the Moon, VSOP87 (and newer) or Newcombe's table of the sun and was wondering what time scale they require as input? I read an explanation given for VSOP87 and it seems to require TT, but what about the Improved Lunar Ephemeris?
TT~TDB
as close enough as per the SLALIB (precursor to SOFA). If you want the full precision, look at DE440t which has the TT-TDB transformation in it and you can iterate until convergence (start with just TT, see what TT-TDB that gives you for the date, add to your date and feed back in to get a refined TT-TDB value). This then only needs a leapsecond table for the UTC-TT part – astrosnapper Aug 20 '22 at 21:30astropy
module, which is quickly becoming the standard used by working astronomers. It provides a Time class for working with all sorts of time standards and an interface for JPL ephemerides. – Paul T. Aug 29 '22 at 16:28