;$Id: pppol2.pro,v 1.1 2019/11/19 19:03:36 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.2,.5]
!y.margin=[3.2,.5]
bar='!20!s!A$!n!r!6'
bar2='!20!s!u$!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='!13P!6!dM!n'
!y.title='!13P!6!dGW!n'
plot_oo,hel1,pol1,ps=1,xr=xr2,yr=yr2,/nodata
oplot,hel1,pol1,ps=8,col=122
oplot,-hel1,-pol1,ps=8,col=55
;xx=[.002,.4] & oplot,xx,1.7*xx
;
print
print,"$mv idl.ps ~/tex/roper/helical/fig/pppol2.ps"
print
!p.multi=0
END