root = '/home/bcarry/data/quick/hebe-hygeia/'


  readcol, root+'date', id, iso, jd, format='(I,A,D)', /Silent
  N = n_elements(id)


  for s=0,1 do begin
    for k=0, N-1 do begin

      idM = strtrim( string(id[k],format='(I)'),2 )
      idS = strtrim( string(s+1,format='(I)'),2 )

      ssoName =  idM+'_'+idS
      model = readVer(  root+ssoName+'.shape' )
      spin  = readSpin( root+ssoName+'.spin'  )

      vGeo = ssoPositions( idM, jd[k], type='aster', /geocentric   )
      vSun = ssoPositions( idM, jd[k], type='aster', /heliocentric )
      sep = subObserver_coord( -vGeo, spin, jd[k], /LTC )
      ssp = subObserver_coord( -vSun, spin, jd[k], /LTC )


      print, idM, idS, iso[k], sep.lon, sep.lat, ssp.lon, ssp.lat, $
             format='(I2,1x,I1,2x,A-22,1x,4(I3,2x))'

      image = shapeShow3d( model, spin, sep=sep, ssp=ssp )

      writefits, root+ssoName+'-'+iso[k]+'.fits', image

    endfor
  endfor





end