;$Id: pcomp_M512k2_large.pro,v 1.2 2021/09/10 05:22:54 brandenb Exp $
if !d.name eq 'PS' then begin
  ;device,xsize=18,ysize=12,yoffset=3
  device,xsize=18,ysize=5,yoffset=3
  !p.charthick=1.5 & !p.thick=1.5 & !x.thick=1.5 & !y.thick=1.5
end
;
!p.charsize=1.0
!x.margin=[8.3,1.5]
!y.margin=[3.1,0.4]
!p.multi=[0,2,1]
!x.range=[1,5.1]
;!x.tickv=[1.,1.5,2.,2.5,3.,3.5]
;!x.tickname=['1','1.5','2','2.5','3','3.5']
;
th1=2
th2=5
xout1=1.11
xout2=1.11
siz=.8
;
;-----------------------------------------------------------------------------
dir1='M512sig1_k6_ramp1a'
dir2='M512sig1_k6_ramp1a_f001'
dir4='M512sig1_k6_ramp1c_f002'
dir5='M512sig1_k6_ramp1b_f007'
;7.0e-03 5.0e-05 4.05e-01 3.00e-05 3.08e-05  2.07 600  M512sig1_k6_ramp1b_f007
;2.0e-03 2.0e-06 9.67e-02 5.62e-07 5.23e-06  2.38 600  M512sig1_k6_ramp1c_f002
;1.0e-03 2.0e-07 3.35e-02 3.53e-08 1.34e-06  2.09 600  M512sig1_k6_ramp1a_f001
;6.0e-04 2.0e-07 1.68e-02 5.31e-09 7.07e-07  1.98 600  M512sig1_k6_ramp1a
;
yrM=[1e-3,.60] & yrGW=[1e-10,5e-5]
;
default,iread,1
if iread then begin
  pc_read_ts,o=ts1,datadir='../'+dir1+'/data'
  pc_read_ts,o=ts2,datadir='../'+dir2+'/data'
  ;pc_read_ts,o=ts3,datadir='../'+dir3+'/data'
  pc_read_ts,o=ts4,datadir='../'+dir4+'/data'
  pc_read_ts,o=ts5,datadir='../'+dir5+'/data'
  iread=0
endif
;
!y.ticks=0
!x.title='!8t!6'
!y.title='!13E!6!dM!n(!8t!6)'
plot_io,ts1.t,ts1.EEM,yr=yrM
oplot,ts2.t,ts2.EEM,col= 55,th=th1
oplot,ts4.t,ts4.EEM,col=155,th=th1
oplot,ts5.t,ts5.EEM,col=122,th=th1
xyouts,xout1,.330,'!6(a)',siz=siz
xyouts,3.0,2.4e-1,'!6A',siz=siz,col=122
xyouts,2.2,1.2e-1,'!6B',siz=siz,col=155
xyouts,2.5,3.5e-2,'!6C',siz=siz,col=55
xyouts,2.7,1.4e-2,'!6D',siz=siz
;
!y.ticks=0
!y.title='!13E!6!dGW!n(!8t!6)'
plot_io,ts1.t,ts1.EEGW,yr=yrGW
oplot,ts2.t,ts2.EEGW,col= 55,th=th1
oplot,ts4.t,ts4.EEGW,col=155,th=th1
oplot,ts5.t,ts5.EEGW,col=122,th=th1
xyouts,xout2,1.6e-5,'!6(b)',siz=siz
;
print,"$mv idl.ps ~/tex/tina/BBN/fig/pcomp_M512k2_large.eps"
!p.multi=0
END