pro plotdst
  filename = '/project/image/serv8/dst/dst200111.txt'
  dst = readdst(filename) 
  filename = '/project/image/serv8/imf/ace_level0_20011105.txt'
  imf = processimf0(filename)
;** Structure <821c99c>, 2 tags, length=7440, refs=1:
;   DST             INT       Array[744]
;   T               DOUBLE    Array[744]
;% Compiled module: READLEVEL0.
;** Structure <821f234>, 6 tags, length=40348, refs=1:
;   BX              FLOAT     Array[1441]
;   BY              FLOAT     Array[1441]
;   BZ              FLOAT     Array[1441]
;   VX              FLOAT     Array[1441]
;   NP              FLOAT     Array[1441]
;   T               DOUBLE    Array[1441]
  days_per_year = julday(1, 1, fix(dst.t)+1, 0, 0, 0)- $
    julday(1, 1, fix(dst.t), 0, 0, 0)
  days_per_year = double(days_per_year)
  t0 = 2001.+308./days_per_year[0]
  t1 = 2001.+310./days_per_year[0]
  xrange = [t0, t1]
;  position = [0.1, 0.1, 0.85, .5]
  plot, dst.t, dst.dst, ystyle = 8, psym = 10, $
    ytitle = 'D!Dst!N (nT)', charsize = 1.2, xrange = xrange, $
    xtickformat = 'year2hhmm', xticks = 12, $
    xminor = 12, xstyle = 1
;; Calculate pressure
;  p = 1.9e-27*imf.np*1e6*imf.vx^2*1e6*1.e9 ;; nPa
;  yrange = [0, 10]
;  axis, yax = 1, yrange = yrange, /save, ytitle = 'Pressure (nPa)', $
;    ystyle = 1, charsize = 1.2
;  oplot, imf.t, p, thick = 2
  yrange = [-30, 20]
  axis, yax = 1, yrange = yrange, /save, ytitle = 'IMF B!Dz!N (nT)', $
    ystyle = 1, charsize = 1.2
  oplot, imf.t, imf.bz, thick = 2
t0 = 2000.+277./days_per_year+06./(24.*days_per_year)
t1 = 2000.+277./days_per_year+09./(24.*days_per_year)
t2 = 2000.+277./days_per_year+17.5/(24.*days_per_year)
t3 = 2000.+278./days_per_year+22./(24.*days_per_year)
plots, [t0, t0], yrange, /data
plots, [t1, t1], yrange, /data
plots, [t2, t2], yrange, /data
plots, [t3, t3], yrange, /data
  return
end