x      = output coordinates of earth or earth-moon barycenter
relative to sun determined from n-body tape
output units are au and au/day
input quantities are
jd     = julian day number
fract  = coordinate time fraction of day
nvel   = 0 x(1-3) position output
nvel   = 1 x(4-6) velocity output
ntype  =-1 moon relative to earth
ntype  = 0 earth-moon barycenter relative to sun
ntype  = 1 earth relative to sun