pro xspacecraft_pv, spacecraft_id, missionDayNumber
;
; Plot uncertainties in r, drdt
;
; B. Knapp, 2003-11-11
;
  case spacecraft_id of
    21701: jdn = ud2jd(missionDayNumber)
    27651: jdn = sd2jd(missionDayNumber)
     else: stop,'Unknown spacecraft_id'
  endcase
;
  spacecraft_pv, spacecraft_id, jdn, pv, coorid=1, $
    tle_path='/tools/knapp/dat/tle/', $
;   tle_path='/mizar/proj/sorce/science_data_system/data/', $
;   tle_path='/home/betelgeuse/knapp/norad_tle/', $
    pv_uncertainty=dpv, interpolation_fraction=f, status=s
;
  get_dev, dev, /portrait
  ans=''
  !p.multi=[0,1,2]
  ang = 360*(findgen(21)/20)/!radeg
  usersym, cos(ang), sin(ang), /fill
  symsiz = 0.1

  plot,missionDayNumber,pv[2,*],psym=8, symsize=symsiz, /ynozero, $
    title='Radius (km)', xtitle='Mission Day'
  plot,missionDayNumber,dpv[2,*],psym=8, symsize=symsiz, $
    title='Radius uncertainty (km)', xtitle='Mission Day'
  if dev eq 0 then read,'CR..',ans
  plot,missionDayNumber,pv[5,*]*1000.d0, psym=8, symsize=symsiz, $
    title='dr/dt (m/s)', xtitle='Mission Day'
  plot,missionDayNumber,dpv[5,*]*1000.d0,psym=8, symsize=symsiz, $
    title='dr/dt uncertainty (m/s)', xtitle='Mission Day'
  send_plots,dev
  !p.multi=[0,0,0]
;
; stop
  return
end