Anda membuat kesalahan mendasar yang sama dengan yang dilakukan Anton Gromov dalam pertanyaannya di situs jaringan Space Exploration StackExchange. Anda menggunakan barycenter tata surya alih-alih Matahari sebagai asal bingkai. Seandainya Anda menggunakan Matahari, perbedaan yang nyata akan turun hampir dua kali lipat.
Kelemahan yang jauh lebih kecil adalah bahwa Anda tampaknya menggunakan tengah malam daripada siang hari. Zaman J2000.0 adalah 12 siang Waktu Terestrial pada 1 Januari 2000. Ini tidak banyak mengubah banyak hal, tetapi tetap penting. Setengah hari offset itu telah menggigit saya lebih dari sekali.
Koordinat z posisi Bumi tidak akan turun ke nol seandainya Anda membuat koreksi di atas (berpusat pada matahari alih-alih barycenter tata surya, dan 12 siang daripada tengah malam). Itu tidak akan jatuh ke nol bahkan jika Anda telah menggunakan barycenter Bumi-Bulan daripada pusat Bumi.
Salah satu alasannya tidak akan jatuh ke nol adalah bahwa orbit Bumi tentang Matahari bukanlah elips. Orbit Bumi bahkan tidak terletak di pesawat! Penyebab utama dari gerakan non-planar ini adalah bahwa orbit Bulan tentang Bumi condong sekitar 5.15 ° terhadap ekliptika, membuat Bumi melonjak naik turun sehubungan dengan ekliptika 1 . Karena orbit Bumi tidak benar-benar elips atau bahkan planar, bidang ekliptika pada beberapa zaman adalah bidang rata-rata waktu yang berpusat pada waktu zaman yang membuat posisi bumi sehubungan dengan Matahari (atau posisi Matahari terhadap Bumi) rata-rata ke nol.
Alasan yang lebih kecil adalah bahwa J2000.0 berarti ekliptika dan rata-rata kerangka ekuinoks didefinisikan lebih dari 36 tahun yang lalu. HORIZONS menggunakan DE431, yang dirilis 6 tahun lalu. Kemampuan JPL untuk memodelkan tata surya meningkat secara signifikan dalam tiga dekade terakhir.
1 Sebenarnya, itu membuat Matahari tampak naik turun sehubungan dengan ekliptika. Konsep bidang ekliptika berasal dari zaman Yunani kuno, dan masih mempertahankan sisa definisi ini hingga hari ini.