;$Id: pcomp.pro,v 1.8 2021/02/19 21:02:03 brandenb Exp $ if !d.name eq 'PS' then begin device,xsize=18,ysize=18,yoffset=3 !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3 end ; ; mv idl.ps ~/tex/axel/vertical/fig/pbzmmax.ps ; !p.charsize=1.7 !x.margin=[9.2,.5] !y.margin=[3.2,.5] ; ; sample to compare time series form different directories ; dir1='M64sig0_k6_wav1d' dir2='M64sig0_k6_wav10_2p34em2a' dir3='M64sig0_k6_wav3_3p2em2' dir4='M64sig0_k6_wav03_2p0em1' ; 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' ; xr=[1,6.3] yr=[0,.3] !p.multi=[0,1,2] !x.range=xr !x.title='!6' !y.title='!7X!6!dM!n' plot_oi,xr,yr,/nodata ; oplot,ts1.t,ts1.brms oplot,ts2.t,ts2.brms,col=122 oplot,ts3.t,ts3.brms,col=55 oplot,ts4.t,ts4.brms,col=155 ; !x.title='!8t!6' !y.title='!7X!6!dGW!n' plot_io,ts1.t,ts1.EEGW,yr=[1e-6,5e-4] oplot,ts2.t,ts2.EEGW,col=122 oplot,ts3.t,ts3.EEGW,col=55 oplot,ts4.t,ts4.EEGW,col=155 ; print,"$mv idl.ps ~/tex/emma/k1GW/fig/pcomp.ps" print,"$mv idl.ps ~/tex/emma/QCD-GWs/fig/pcomp.ps" END