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=480,ysize=480
end
;
; mv idl.ps ~/rmhd/other/anvar/cluster/structure_RM.ps
;
!p.charsize=1.4
!p.multi=[0,1,1]
;
; Some values for the initial point
;
t0=2517.55
urms0=0.123
Brms0=0.063
to=1/urms0
kf=1.5
;
; Read in data
;
file='histogram_t=2520.01_xx.sav'
restore,file=file
SF2_2520=SF2
Brms_2520=6.29e-2
norm_fact_2520=(Brms_2520/kf)^2.
norm_fact_2520=mean_2RM02*256.^2
;
file='histogram_t=2760.01_xx.sav'
restore,file=file
SF2_2760=SF2
Brms_2760=1.41e-2
norm_fact_2760=(Brms_2760/kf)^2.
norm_fact_2760=mean_2RM02*256.^2
;
file='histogram_t=3111.99_xx.sav'
restore,file=file
SF2_3111=SF2
Brms_3111=6.42e-3
norm_fact_3111=(Brms_3111/kf)^2.
norm_fact_3111=mean_2RM02*256.^2
;
; Define r vector
;
nx=256
r=indgen(nx/2)*!pi/(nx/2-1)
;
; Plot results
;
!y.title='!8c!6(!8R!6)'
!x.title='!8R!6'
;
plot,r,1-SF2_2520/norm_fact_2520,yr=[-0.25,1.05],ticklen=0.03
oplot,r,1-SF2_2760/norm_fact_2760,li=2
oplot,r,1-SF2_3111/norm_fact_3111,li=1
;
!x.title=''
!y.title=''
;
END