pro test_msis

;  msis00, date, sec, glat, glon, alt, t_alt, t_exo, nd, md, $
;          ft7p=<ft7p_value>, ft7a=<ft7a_value>, apd=<ap_value[s]>, $
;          /SI, /Ap_3hr, /HELP

  m = 20L
  date = lonarr(m)
  date[0] = 2002001
  delta_date = 1
  for i=1,m-1 do date[i] = date[i-1]+delta_date
  s = 0
  lat = -80.
  lon = 0.
  alt = 70.
  h = fltarr(m)
  he = fltarr(m)
  n = fltarr(m)
  o = fltarr(m)
  n2 = fltarr(m)
  o2 = fltarr(m)
  ar = fltarr(m)
  ao = fltarr(m)
  t = fltarr(m)

  for j=0,m-1 do begin
;   msis00, date[j], s, lat, lon, alt, t_alt, t_exo, nd, md, $
;     ft7p=80., ft7a=80., apd=4.
    msis00, date[j], s, lat, lon, alt, t_alt, t_exo, nd, md
    h[j] = nd.h
    he[j] = nd.he
    n[j] = nd.n
    o[j] = nd.o
    n2[j] = nd.n2
    o2[j] = nd.o2
    ar[j] = nd.ar
    ao[j] = nd.ao
    t[j] = t_alt
;   print,'j=',j,' date=',date[j],' h=',h[j],' he=',he[j],' n=',n[j],' o=',o[j],$
;     ' n2=',n2[j],' o2=',o2[j],' ar=',ar[j],' ao=',ao[j],' t=',t[j]
    
  endfor
   
end