;  Plot EM vs time


restgenx, file='ospex_results_nov3_multitherm.geny', results

;emission measure in cm^-3
a0 = results.spex_summ_params[15,*]
a0 = reform(a0)
aone = results.spex_summ_params[16,*]
aone = reform(aone)
a2 = results.spex_summ_params[17,*]
a2 = reform(a2)
a3 = results.spex_summ_params[18,*]
a3 = reform(a3)

;calculate emission measure from DEM
em = (2./ (a3-1.)) * a0 * ((aone/2.)^(1.-a3) - (a2/2.)^(1.-a3))
;em = em * 10.^49


;start and end time in sec
start_time = results.spex_summ_time_interval[0,*]
start_time = reform(start_time)
end_time = results.spex_summ_time_interval[1,*]
end_time = reform(end_time)
tarr = (start_time + end_time)/2.

; Time range to plot
s = 3 & e = n_elements(tarr) - 1
stime = anytim(tarr(s),/vms)
etime = anytim(tarr(e),/vms)

utb = min(tarr)
; Plot emission measure
utplot,tarr(s:e)-utb,em(s:e),utb, $
   /ylog,ytitle='Emission Measure in units of 10!u49!n cm!u-3!n', $
    xtitle = ' ', title = 'Plot of Emission Measure vs. TIme'
xyouts,0.3,0.9,'Start Time = ' + stime, color=0,/normal
xyouts,0.3,0.85,'End Time = ' + etime, color=0,/normal

end