cwd,run
fo2="(a,' & ',a,' & ',a,' & ',a,' & ',a,' \\% ', a)"
fo2="(a,' & ',a,' & ',a,' & ',a,' \\% ', a)"
fo='(5e10.1)'
nt=n_elements(tt)
tA=xiM/vA
default,CM,20.
t0=CM*tA-tt
print,'CM=',CM
for it=0,nt-1 do begin
  ;print,tt[it],vA[it],xiM[it],xiM[it]/vA[it],tt[it]*vA[it]/xiM[it],fo=fo
  ;print,latex(tt[it],-1),latex(vA[it],1),latex_float(reform(xiM[it]),1),latex_float(reform(t0[it]),1),latex_float(reform((tt[it]+t0[it])/tA[it]),1),dir,fo=fo2
  print,latex(tt[it],-1),latex(vA[it],1),latex_float(reform(xiM[it]),1),latex(reform(t0[it]),1),dir,fo=fo2
endfor
END