ep='2016-10-10'
  nbd=13
  step=30


  jd = dblarr(nbd)
  jd[0] = date_conv(ep,'JULIAN')
  iso = strarr(nbd)

  planet = replicate({x:0.,y:0.,r:0.,t:0.},8,nbd)

  for kEp=0, nbd-1 do begin
    jd[kEp]  = jd[0]+ kEp*step
    iso[kEp] = strMid(date_conv(jd[kEp],'FITS'),0,10)
  endfor

  for kP=0, 7 do for kEp=0, nbd-1 do begin

    sso = ssoPositions( kP+1, jd[kEp], /planet, rPlane=2, /heliocentric )

    planet[kP,kEp].x = sso[0]
    planet[kP,kEp].y = sso[1]

    planet[kP,kEp].r = sqrt( sso[0]*sso[0] + sso[1]*sso[1] )
    planet[kP,kEp].t = ( atan( sso[1], sso[0] )/!DTOR + 720 ) mod 360

  endfor

  forprint, iso, textout='planet.csv', format='(A10,8(", ",I4))', $
            planet[0,*].t, planet[1,*].t, planet[2,*].t, planet[3,*].t, $
            planet[4,*].t, planet[5,*].t, planet[6,*].t, planet[7,*].t, $
            comment='Epoque, Mercure, Venus, Terre, Mars, Jupiter, Saturne, Uranus, Neptune'


end