[x50,x51]

fr
?Doc file from AP MQ?
;<1> docapmq
fr
? name of input particles?
;<2> original particles
fr
? name of rotated particles?
;<3> rotated particles 

;  x51 - number of reference images

do lb1 i=1,x50
x55=x0
ud ic,x55,x21,x22,x23,x24,x25,x26
<1>
if(x21.gt.x51) then
rt sq
<2>x26
_1
x23
x24,x25
mr
_1
<3>x26
Y
else
rt sq
<2>x26
<3>x26
x23
x24,x25
endif
lb1

ud ice
<1>
re