* Plots data from Atkinson, et al, Z Phys C27(1985)233

opt *
set *
igset *
opt nbox

vec/delete *
hist/delete *

fort/file 66 phipol.eps
metafile 66 -113

set xsiz 20
set ysiz 10
set xwin  2

vec/create cth(9) R -0.85 -0.625 -0.40 -0.25 0 0.25 0.4 0.625 0.85
vec/create dth(9) R 175 350 650 770 640 720 500 310 90
vec/create eth(9) R 50 70 80 80 80 80 70 60 30
vec/create psi(9) R 20 60 100 140 180 220 260 300 340
vec/create dps(9) R 550 320 315 510 610 490 375 430 570
vec/create eps(9) R 80 60 60 60 60 60 60 60 75
vec/create zer(9) R 0 0 0 0 0 0 0 0 0

zone 2 1 1
hplot/null -1 1 0 900
hplot/error cth dth zer eth 9 20
set ltyp 0 ; func/plot 675*(1-x**2) -1 1 CS

hplot/null 0 360 0 720
hplot/error psi dps zer eps 9 20
set ltyp 0 ; func/plot 350*(1+0.75*cos(x*3.14/180)**2) 0 360 CS

close 66