pro pltmgii,ps=ps
;
; Plot MgII Core/Wing data from solstice_mgii.dat
; (obtained via anon ftp at rescha.colorado.edu)
;
filename = "solstice_mgii.dat"
;
; read_dat.pro (in cwd) obtained from Barry Knapp, 6/99:
;
mgii = read_dat(filename,status=istat)
; MGII            DOUBLE    = Array[3, 2414]
fdat = mgii(2,*)
fdat(where(fdat eq 0)) = 1.e36
fmax = max(fdat(where(fdat ne 1.e36)))
fmin = min(fdat(where(fdat ne 1.e36)))
print,'fmin=',fmin,' fmax=',fmax
ndays = n_elements(fdat)
print,'ndays=',ndays
udays = mgii(0,*)
title="SOLSTICE Mg-II core/wing ratio time series"
if keyword_set(ps) then begin
  psfile = 'MgII.ps'
  print,'Will make ps file ',psfile
  set_plot,'PS'
  device,filename=psfile
endif
plot,udays,fdat,max_value=1.e35,title=title,$
  ytitle='Core/Wing',yrange=[.1400,.1900],/ystyle,$
; ytitle='Core/Wing',/ystyle,$
  xtitle='UARS Day',/xstyle
if (keyword_set(ps)) then device,/close_file
end