;$Id: pcomp_nlin.pro,v 1.2 2021/06/16 21:06:16 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
!p.charsize=1.7
!x.margin=[7.8,.5]
!y.margin=[3.2,.5]
print,"$mv idl.ps ~/tex/axel/vertical/fig/pbzmmax.ps"
;
;  sample to compare time series form different directories
;
dir1='M512_k1_kf100_nlin0'
dir2='M512_k1_kf100_nlin8'
dir3='M512_k1_kf100_nlin4'
;
dir1='K64irrok6c_StT_cont_tst0'
dir2='K64irrok6c_StT_cont_tst3'
;
pc_read_ts,o=ts1,datadir='../'+dir1+'/data'
pc_read_ts,o=ts2,datadir='../'+dir2+'/data'
;pc_read_ts,o=ts3,datadir='../'+dir3+'/data'
;
!p.multi=[0,1,3]
!p.multi=0
!x.range=[0,3]
;
yr=[0.,4e-6]
yr=[0.,4e-0]
plot,ts1.t,ts1.hrms,yr=yr
oplot,ts2.t,ts2.hrms,col=122,li=2
;oplot,ts3.t,ts3.hrms,col=55,li=2
;
END