if !d.name eq 'PS' then begin
  device,xsize=15,ysize=13,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=480,ysize=480
end
;
; mv idl.ps structure_RM.ps
;
!p.charsize=1.6
!p.multi=[0,1,1]
;
; Some values for the initial point
;
t0=70.
urms0=0.135
Brms0=5e-3
to=1/urms0
kf=1.5
;
; Read in data
;
file='histogram_t=  35.00_xx.sav'
restore,file=file
SF2_35=SF2;mean_2RM02
Brms_35=1.93e-3
norm_fact_35=(Brms_35/kf)^2
norm_fact_35=mean_2RM02*256.^2
;
file='histogram_t=  70.00_xx.sav'
restore,file=file
SF2_70=SF2;mean_2RM02
Brms_70=4.71e-3
norm_fact_70=(Brms_70/kf)^2
norm_fact_70=mean_2RM02*256.^2
;
file='histogram_t= 105.00_xx.sav'
restore,file=file
SF2_105=SF2;mean_2RM02
Brms_105=3.87e-3
norm_fact_105=(Brms_105/kf)^2
norm_fact_105=mean_2RM02*256.^2
;
; Define r vector
;
nx=256
r=indgen(nx/2)*!pi/(nx/2-1)
;
; Plot results
;
!y.title='!6<[RM(!8x+r!6)-RM(!8x!6)]!u2!n>/<2[RM(!8x!6)]!u2!n>'
!x.title='!8r!6'
;
plot,r,SF2_35/norm_fact_35,yr=[0,2.]
oplot,r,SF2_70/norm_fact_70,li=2
oplot,r,SF2_105/norm_fact_105,li=3
;
; Legends
;
xx=0.2 & dx=0.3
yy=1.6 & dy=-0.15
leg_text='!8t!6='+string((35-t0)/to,FORMAT='(F5.2)')+'t!dto!n'
legend,xx,dx,yy+dy*0,0,leg_text
leg_text='!8t!6='+string((70-t0)/to,FORMAT='(F5.2)')+'t!dto!n'
legend,xx,dx,yy+dy*1,2,leg_text
leg_text='!8t!6='+string((105-t0)/to,FORMAT='(F5.2)')+'t!dto!n'
legend,xx,dx,yy+dy*2,3,leg_text
;
!p.title=''
;
END