$Id: p_antisym.pro,v 1.1 2018/03/20 13:30:04 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
;
;  mv idl.ps ~/tex/sayan/Helicity-Realizability/axel/fig/p_antisym1.ps
;  mv idl.ps ~/tex/sayan/Helicity-Realizability/axel/fig/p_antisym2.ps
;
!p.charsize=1.7
!x.margin=[7.8,.5]
!y.margin=[3.2,.5]
;
pc_read_var,obj=var
u2=dot2(var.uu)
power3d,u2,wav=k,spec=spec,/noplot
power3d_antisym,u2,wav=k,Sspec=S,Aspec=A,/noplot
;
!x.title='!8k!6'
!y.title='!6'
plot,k,spec,yr=[-1,1]*7e-7
oplot,k,A,col=55
oplot,k,S,col=122
;oplot,k,S+A,col=188,li=2
END