macro third
exec init1
flag=0

set XLAB 2.5
set XMGL 3.5

if [flag].eq.1 then

nt/cre 10 ' ' 4 ' ' 10000 t p dt dp
nt/rea 10 scat1

set MTYP 7

null 0 180 -.9 .3
2d 100 ' ' 110 0 180 110 -.9 .3
nt/pl 10.dt%t ! ! ! ! n 100

set ncol 255

h/pl 100 's col'

endif

zone 1 2

nt/cre 10 ' ' 4 ' ' 10000 t p dt dp
nt/rea 10 scat_small4

null 0 180 -1.5 1.5
nt/pl 10.dp%t ! ! ! ! s

atit ' ' 'deviation in [f] from direction of primary'


opt linx
opt liny
null 0 1 0 1 sa

set TXCI 4
set chhe 0.5
itx .3 .9 with magnetic field


nt/cre 10 ' ' 4 ' ' 10000 t p dt dp
nt/rea 10 scat_small1

null 0 180 -1.5 1.5
nt/pl 10.dp%t ! ! ! ! s

set TXCI 1
atit '[q] "M# degrees "N#'


opt linx
opt liny
null 0 1 0 1 sa

set TXCI 4
set chhe 0.5
itx .3 .9 w/o magnetic field


exec plot1
return










macro second
exec init1
flag=2

nt/cre 10 ' ' 4 ' ' 10000 t p dt dp
nt/rea 10 scat_small4

null -1 1 -1 1

set MTYP 7

set ncol 58
palette 1

if [flag].eq.1 then

appl comis quit
        real function fput
        real t, p, dt, dp
        logical chain
        character*128 cfile
        common /pawchn/ chain,nchevt,ichevt
        common /pawchc/ cfile
        common /pawidn/ idnevt,obs(13), t, p, dt, dp
	dimension xx(1), yy(1)
	if ( p.gt.180 ) then
		p=360-p
	endif
	g=float(nint(p/4)+10)
	call igset('PMCI', g)
	xx(1)=dp
	yy(1)=dt
	call ipm(1, xx, yy)
        end
quit

nt/lo 10 fput

elseif [flag].eq.2 then

do i=360,1,-1

if [i].le.180 then
	n=[i]
else
	n=360-[i]
endif

set PMCI [n]/4+10

min=([i]-1)
max=[i]

nt/pl 10.dt%dp p.gt.[min].and.p.lt.[max] ! ! ! s

enddo

endif

atit 'deviation in [f] from direction of primary' 'deviation in [q] from direction of primary'


opt linx
opt liny
null 0 1 0 1 sa

set TXCI 4
set chhe 0.5
itx .3 .2 with magnetic field


itx .2 .8 [f] color map:
appl comis quit
	subroutine phi
	data x/0.65/, y/0.8/

	call isfais(1)
	do i=1,360
	if ( i.le.180 ) then
		n=i
	else
		n=360-i
	endif
	fmin=float(i-1)
	fmax=float(i)
	call isfaci(n/4+10)
	call igarc(x, y, 0.0, 0.1, fmin, fmax)
	enddo

	call isfaci(10)
	call igarc(x, y, 0.0, 0.2, -1., 1.)

	call isfaci(17)
	call igarc(x, y, 0.0, 0.2, 31., 33.)

	call isfaci(1)
	call igarc(x, y, 0.172, 0.175, 4., 28.)

	end
	call phi
	end
quit

x=0.8032
y=0.8815

x1=[x]+.02
y1=[y]-.02

x2=[x]
y2=[y]-.0283

dline [x] [x1] [y] [y1]
dline [x] [x2] [y] [y2]

itx 0.83 0.84 [f]

exec plot1
return










macro first
exec init1
pi=3.14159/180

opt logy
null -1 1 2.e3 3e4

vec/rea a fullcos4
set dmod 1
vec/pl a ! s

vec/del a
vec/rea a fullcos1
set dmod 2
vec/pl a ! s



opt linx
opt liny
null 0 1 0 1 sa

off=-0.63
set TXCI 4
set chhe 0.35
set dmod 1
dline .65+[off] .78+[off] .907 .907
itx .8+[off] .90 ver. 6.016 curved surface with [q] "G# 90 with magnetic field
set dmod 2
dline .65+[off] .78+[off] .857 .857
itx .8+[off] .85 ver. 6.016 curved surface with [q] "G# 90 w/o magnetic field

* set chhe 0.5
* itx .3 .2 with magnetic field

atit 'cos [q]' 'entries'

exec plot1
return










macro main
* exec read

pi=3.14159/180
fl=3


if [fl].eq.1 then

opt liny
null -1 1 -1 1
nt/pl 10.b%c a.lt.80.and.a.gt.78 ! ! ! s

elseif [fl].eq.2 then

opt logy
null -1 1 100 20000
set dmod 1
nt/pl 10.cos(a*[pi]) ! ! ! ! s

elseif [fl].eq.3 then

opt logy
null 0 1 10 20000
set dmod 1
vec/pl a ! s

endif


* exec plot1
return


macro readnt
exec init1
nt/cre 10 ' ' 3 ' ' 1000 a b c
nt/rea 10 scat1
return


macro readvec
exec init1

pi=3.14159/180

vec/rea a,b,c scat8
sigma a=cos(a*[pi])
return


macro init1
set *
opt *
pic/del *
vec/del *
hi/del *
opt zfl
opt NBOX
set XMGL 3.0
set YMGL 2.4
set XVAL 0.8
set YVAL 0.4
set XLAB 2.0
set YLAB 1.4
set KSIZ .50
set GSIZ .50
set TSIZ .50
set ASIZ .50
set CSIZ .50
set PSIZ .50
set VSIZ .50
set SSIZ .50
set 2SIZ .50
igset LASI .50
igset TMSI .50
set HWID 5
set BWID 1
set PWID 5
set FWID 5
set XWID 1
set YWID 1
set CWID 5
igset LWID 5
igset TXFP -130
set CHHE .50
set cfon -40
set gfon -60
set lfon -40
set tfon -40
set vfon -40
return


macro plot1
fort/file 66 a.ps
meta 66 -111
pict/plot *
close 66
return