;$Id: ppvar_power_comp.pro,v 1.2 2014/03/03 05:54:03 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=2 & !p.thick=2 & !x.thick=2 & !y.thick=2
end
;
;  mv idl.ps ~/tex/tina/Decay/fig/ppvar_power_comp.ps
;  mv idl.ps ~/tex/tina/Decay/fig/ppvar_power_comp2.ps
;
!p.charsize=1.7
!x.margin=[7.0,0.5]
!y.margin=[3.2,2.5]
;
dir1='short2304pm1_kf60b_noforce'
dir2='short288pm1_kf60a_noforce'
dir2='1024_4x32x32_Pm1'
file1='../'+dir1+'/data/VAR3_16_power.sav'
file2='../'+dir2+'/data/VAR1_2_power.sav'
file2='../'+dir2+'/data/VAR5_8_power.sav'
;!p.multi=[0,1,2]
yr=[5e-8,1e-4]
;
restore,file1
plot_oo,k(1:*),EM(1:*),/nodata,xtit='!8k!6/!8k!6!d1!n',yr=yr
oplot,k(1:*),EM(1:*),col=122
oplot,k(1:*),EK(1:*),col=122,li=2
xyouts,col=122,4,3e-6,file1
;
s=.08
s=.03
restore,file2
;plot_oo,k(1:*),EM(1:*),tit=dir2,xtit='!8k!6/!8k!6!d1!n'
oplot,k(1:*),s*EM(1:*),col=55
oplot,k(1:*),s*EK(1:*),col=55,li=2
xyouts,col=55,4,1e-7,file2
;
xx=[1.1,7] & oplot,xx,1e-7*xx^4
!p.multi=0
END