;$Id: pppol1.pro,v 1.7 2020/11/01 15:13:44 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=10,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
!p.charsize=1.5
!x.margin=[8.0,.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'
;
;  read file
;
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]*1.05
yr=[-1,1]*1.05
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,hel2,pol2,ps=1,xr=xr,yr=yr,/nodata
oplot,hel1,pol1,ps=8,col=122
;oplot,hel2,pol2,ps=8,col=55
;xx=[.002,.4] & oplot,xx,1.7*xx
;
xxx=2*grange(-1.,1.,200)
oplot,xxx,.6*tan(xxx)+.03
xxx=2*grange(0.,1.,200)
;oplot,xxx,.8*tan(xxx)-.05,col=55,li=1
oplot,xxx,.7*tan(xxx),col=55
;
xxx=2*grange(-1.,1.,200)
hel3=2.*xxx/(1.+xxx^2)
pol3=4.*xxx*(1.+xxx^2)/((1.+xxx^2)^2+4.*xxx^2)
oplot,hel3,pol3,li=1
;
;  read file
;
file='pol_ini.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,*])
oplot,hel1,pol1,ps=8,col=55
;
print
print,"$mv idl.ps ~/tex/roper/helical/fig/pppol1.ps"
print
!p.multi=0
END