pro eclipse_plotter, eclipse, date=date
if n_elements(date) eq 0 then date = '24-oct-95'
wdef,0
ss  = where(eclipse.daylight eq 1)
utplot,eclipse,eclipse.dd,yr=[0,7000],$
  ytit='Distance between centers, arc sec',$
  tit = 'Yohkoh eclipses for '+date,xstyle=1
outplot,eclipse(ss),eclipse(ss).dd,psym=2
outplot,eclipse,eclipse.r_sun+eclipse.r_moon,lines=1
fem_grid 
wdef,1
nn = findgen(1001)*2*!pi/1000
xx = 961*cos(nn)
yy = 961*sin(nn)
plot,eclipse.xx,eclipse.yy,xr=[-3000,3000],yr=[-3000,3000],$
  tit = 'Heliocentric motion (R_o, 2*R_o)'
oplot,eclipse(ss).xx,eclipse(ss).yy,psym=2
oplot,xx,yy
oplot,2*xx,2*yy,lines=1
;stop
end