;
pro getmgii,ps=ps
;
; Get MgII time-series and find uars day of min and max
; (these days will be used as solar min and max)
;
print,'Calling get_l3_ts for MgII_279.6'
get_l3_ts,'MgII_279.6',data,status
print,'Status=',status
if status ne 0 then begin
  print,'Error ',status,' from get_l3_ts'
  return
endif
fdat = data.flux
uday = data.uday
;
fmax = max(fdat,imax)
daymax = uday(imax)
uday = uday(where(fdat ne 0))
fdat = fdat(where(fdat ne 0))
fmin = min(fdat,imin)
daymin = uday(imin)
;
; Flux data (getbrass.pro) not available for 1st 5 min days,
; so 6th min day = 1417 was used.
;
;fdat(imin) = 1.e36
;fmin = min(fdat,imin)
;daymin = uday(imin)
;
;fdat(imin) = 1.e36
;fmin = min(fdat,imin)
;daymin = uday(imin)
;
;fdat(imin) = 1.e36
;fmin = min(fdat,imin)
;daymin = uday(imin)
;
;fdat(imin) = 1.e36
;fmin = min(fdat,imin)
;daymin = uday(imin)
;
;fdat(imin) = 1.e36
;fmin = min(fdat,imin)
;daymin = uday(imin)
;
fdat = data.flux
uday = data.uday
fdat(where(fdat eq 0)) = 1.e36
ndays = n_elements(fdat)
print,'ndays=',ndays
print,'  max MgII =',fmax,' on day ',daymax
print,'  min MgII =',fmin,' on day ',daymin
;
if keyword_set(ps) then begin
  psfile = 'MgII.ps'
  print,'Will make ps file ',psfile
  set_plot,'PS'
  device,filename=psfile
endif
title='SOLSTICE ' + data.name
plot,uday,fdat,max_value=1.e35,yrange=[.115,.155],xrange=[0.,1600.],$
  /xstyle,/ystyle,title=title,ytitle='Core/Wing',xtitle='UARS Day'
if (keyword_set(ps)) then device,/close_file
return
end