; generate the angular file for random conical tilted images from mra
; writes all the angles psi,theta and phi into a document file.
do lb1 x19=1,4726       ;for each particle
ud ic x19,x21          ;get particle number and micrograph it comes from
particle_list
ud s,124,x50,x51,x52        ;load theta gamma and phi from picking into memory
dcb{***x21}     ;
ud ic x19,x99,x14          ;don't care about shift!!!  just load rotation
ali_50
                ;
x14=-x14            ;change of sign due to imagic spider convention difference!
                ;if the in-plane rotation alignment is done in spider, mark
                ;this line off!
x14=-x14            ;AGAIN! this is following spider batch. We are rotating the
                ;original volume (thus the reference to match the particles)
x14=x14-x51         ;put gamma and rotation together
if(x14.gt.0.0)goto lb20     ;make positive
x14=x14+360.0
lb20
if(x14.lt.360.0)goto lb21
x14=x14-360.0
lb21
sd x19,x52,x50,x14      ;save angles into angular file: psi, theta, phi
angular_file

lb1
ud ice
ali_50
ud ice
particle_list
en