Re: helio(bary)centeric correction (Thorstensen) Dear Dr. Thorstensen, Since I haven't had a conventional analytical formula for obtaining the barycenter, I have written some code in order to adapt a free software using DE200 (NIPE, written by a Japanese amateur Novas) to my purpose. XYZ coordinates for the earth and the barycenter are given for some epochs close to the DE200 initial epoch (I will later make numerical integrations for more recent epochs on a faster machine). I hope Dr. Thorstensen can directly compare these results with his computation. Note JDs are given (in natural framework of DE200) in TDT. Comparison with my old Newcomb-based program has yielded differences to an order of a few 10^-7 AU for the coordinate of the earth. JD(TDT) X2000.0 Y2000.0 Z2000.0 ---------------------------------------------------------------- 2440404.00000 0.174605559 -0.918888826 -0.398457537 center of grav -0.004503521 -0.000791360 -0.000276299 2440404.20000 0.177936955 -0.918353811 -0.398225610 center of grav -0.004503445 -0.000792393 -0.000276744 2440404.40000 0.181266392 -0.917808505 -0.397989225 center of grav -0.004503369 -0.000793427 -0.000277189 2440404.60000 0.184593836 -0.917252912 -0.397748385 center of grav -0.004503293 -0.000794461 -0.000277634 2440404.80000 0.187919251 -0.916687036 -0.397503091 center of grav -0.004503217 -0.000795494 -0.000278079 2440405.00000 0.191242603 -0.916110881 -0.397253344 center of grav -0.004503141 -0.000796527 -0.000278524 2440405.20000 0.194563856 -0.915524452 -0.396999145 center of grav -0.004503064 -0.000797561 -0.000278969 2440405.40000 0.197882976 -0.914927751 -0.396740497 center of grav -0.004502987 -0.000798594 -0.000279414 2440405.60000 0.201199926 -0.914320783 -0.396477399 center of grav -0.004502909 -0.000799627 -0.000279859 2440405.80000 0.204514673 -0.913703551 -0.396209855 center of grav -0.004502832 -0.000800660 -0.000280304