;$Id: pppol2_log.pro,v 1.2 2020/07/31 04:00:56 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=[8.6,.5]
!y.margin=[3.2,.5]
bar='!20!s!A$!n!r!6'
bar2='!20!s!u$!n!r!6'
tilde='!9!s!aA!n!r!6'
;
;  sample to compare time series form different directories
;
file='pol.txt'
a=rtable(file,head=1,6)
sig=reform(a[0,*])
pol1=reform(a[2,*])
pol2=reform(a[3,*])
hel1=reform(a[4,*])
hel2=reform(a[5,*])
;
xr=[-1,1]
yr=[-1,1]
xr2=[.0005,1.2]
yr2=[.0005,1.2]
xr3=[-1,1]*.4
yr3=[-1,1]*.4
xr3=[-1,1]
yr3=[-1,1]
circ_sym,1.3,1
;
!x.title=tilde+'!13P!6!dM!n'
!y.title=tilde+'!13P!6!dGW!n'
plot_oo,hel2,pol2,ps=1,xr=xr2,yr=yr2,/nodata
oplot,hel2,pol2,ps=8,col=122
oplot,-hel2,-pol2,ps=8,col=55
xx=[.002,.4] & oplot,xx,1.7*xx
;
xxx=grange(-1.,1.,200) & oplot,xxx,4.*xxx*(1.+xxx^2)/((1.+xxx^2)^2+4.*xxx^2)
print
print,"$mv idl.ps ~/tex/roper/helical/fig/pppol2_log.ps"
print
!p.multi=0
END