if !d.name eq 'PS' then begin
  device,xsize=14,ysize=10,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
if !d.name eq 'X' then begin
  window,0,xsize=640*9/10,ysize=480*9/10
end
;
; mv idl.ps ~/rmhd/other/anvar/cluster/histogram2.ps
;
!p.multi=[0,1,1]
;
; Some values for the initial point
;
t0=2517.55
urms0=0.123
Brms0=0.063
to=1/urms0
;
; Read in data
;
file='histogram_t=2520.01_xx.sav'
restore,file=file
xvar_2520_xx=xvar
;yvar_2520_xx=yvar
binwidth=xvar(1)-xvar(0)
norm_fact=total(yvar)*binwidth
yvar_2520_xx=yvar/norm_fact
Brms_2520=6.29e-2
time_2520=time-t0
;
file='histogram_t=2580.00_xx.sav'
restore,file=file
Brms_2580=3.82e-2
xvar_2580_xx=xvar
;yvar_2580_xx=yvar
binwidth=xvar(1)-xvar(0)
norm_fact=total(yvar)*binwidth
yvar_2580_xx=yvar/norm_fact
time_2580=time-t0
;
file='histogram_t=2640.00_xx.sav'
restore,file=file
xvar_2640_xx=xvar
;yvar_2640_xx=yvar
binwidth=xvar(1)-xvar(0)
norm_fact=total(yvar)*binwidth
yvar_2640_xx=yvar/norm_fact
Brms_2640=2.41e-2
time_2640=time-t0
;
file='histogram_t=2700.00_xx.sav'
restore,file=file
xvar_2700_xx=xvar
;yvar_2700_xx=yvar
binwidth=xvar(1)-xvar(0)
norm_fact=total(yvar)*binwidth
yvar_2700_xx=yvar/norm_fact
Brms_2700=1.78e-2
time_2700=time-t0
;
file='histogram_t=2760.01_xx.sav'
restore,file=file
xvar_2760_xx=xvar
;yvar_2760_xx=yvar
binwidth=xvar(1)-xvar(0)
norm_fact=total(yvar)*binwidth
yvar_2760_xx=yvar/norm_fact
Brms_2760=1.41e-2
time_2760=time-t0
;
file='histogram_t=2820.01_xx.sav'
restore,file=file
xvar_2820_xx=xvar
;yvar_2820_xx=yvar
binwidth=xvar(1)-xvar(0)
norm_fact=total(yvar)*binwidth
yvar_2820_xx=yvar/norm_fact
Brms_2820=1.18e-2
time_2820=time-t0
;
file='histogram_t=2880.00_xx.sav'
restore,file=file
xvar_2880_xx=xvar
;yvar_2880_xx=yvar
binwidth=xvar(1)-xvar(0)
norm_fact=total(yvar)*binwidth
yvar_2880_xx=yvar/norm_fact
Brms_2880=1.01e-2
time_2880=time-t0
;
file='histogram_t=2940.00_xx.sav'
restore,file=file
xvar_2940_xx=xvar
;yvar_2940_xx=yvar
binwidth=xvar(1)-xvar(0)
norm_fact=total(yvar)*binwidth
yvar_2940_xx=yvar/norm_fact
Brms_2940=8.87e-3
time_2940=time-t0
;
file='histogram_t=3000.00_xx.sav'
restore,file=file
xvar_3000_xx=xvar
;yvar_3000_xx=yvar
binwidth=xvar(1)-xvar(0)
norm_fact=total(yvar)*binwidth
yvar_3000_xx=yvar/norm_fact
Brms_3000=7.82e-3
time_3000=time-t0
;
; Plot results
;
!x.range=[-13,13]
!y.range=[1e-4,2.2e-1]
!p.charsize=1.4
!x.title='RM/(!4q!6!d0!n!8B!6!drms!n/!8k!d!6f!n)'
!y.title='!8N/N!6!d0!d'
plot,xvar_2580_xx,yvar_2580_xx,background=255,col=1,li=2,$
  /ylog
oplot,xvar_2520_xx,yvar_2520_xx,col=1,li=0
oplot,xvar_2940_xx,yvar_2940_xx,col=1,li=1
;
!p.title=''
!x.range=''
!y.range=''
!p.charsize=''
!x.title=''
!y.title=''
!p.title=''
;
END