$Id: pcomp_ko.pro,v 1.7 2022/10/19 02:47:02 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=19,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ~/tex/roper/GW/fig/pcomp.ps
;
!p.charsize=1.7
!x.margin=[11.0,.5]
!y.margin=[3.2,.5]
;
;  sample to compare time series form different directories
;
diro='K512irrok6c_rep1'
dir1='K512irrok6c_rep1_cont'
dir2='K512sig0_k6_ramp1a_cont'
dir3='K512irro_k6_ramp1c_rep'
dir4='K512irro_k6_dtf0p05e' & s4=10.
dir4='K512irro_k6_dtf0p05h2' & s4=1.
;
pc_read_ts,o=tso,datadir='../'+diro+'/data'
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'
;
!p.multi=0
!p.multi=[0,1,2]
xr=[1.,3.6]
yr=[.0,.017]
yr2=[.0,2e-8]
yr2=[.0,3e-8]
;
!x.title="!6"
!y.title="!6<!8u!6!u2!n>"
plot,ts1.t,ts1.urms^2,/nodata,xr=xr,yr=yr
oplot,tso.t,tso.urms^2,col=55
oplot,ts1.t,ts1.urms^2,col=55
oplot,ts2.t,ts2.urms^2,col=122
loadct,6
oplot,ts3.t,ts3.urms^2,col=122
loadct,5
oplot,ts4.t,s4*ts4.urms^2,col=155
oplot,[2,2],yr,li=1
;
!x.title="!8t!6"
!y.title="!13E!6!dGW!n"
plot,ts1.t,ts1.EEGW,/nodata,xr=xr,yr=yr2
oplot,tso.t,tso.EEGW,col=55
oplot,ts1.t,ts1.EEGW,col=55
oplot,ts2.t,ts2.EEGW,col=122
loadct,6
oplot,ts3.t,ts3.EEGW,col=122
loadct,5
oplot,ts4.t,s4^2*ts4.EEGW,col=155
oplot,[2,2],yr2,li=1
;
;xx=[6,30] & oplot,xx,3*xx^2,col=155
;xx=[3,30] & oplot,xx,.1*xx,col=155
;xyouts,20,1.,siz=siz,'slope=1',col=155
;xyouts,9,2e3,siz=siz,'slope=2',col=155
;
print,'$mv idl.ps ~/GitHub/Axel/VortexGen/fig/pcomp_ko.eps'
END