22.05.2014 Views

CERN Program Library Long Writeup W5013 - CERNLIB ...

CERN Program Library Long Writeup W5013 - CERNLIB ...

CERN Program Library Long Writeup W5013 - CERNLIB ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

DIMENSION ROTMAT(3,3), ROWMAT(3), PHI(3), THETA(3)<br />

LOGICAL ROTATE<br />

.<br />

.<br />

.<br />

DO 10 I=1,3<br />

ROWMAT(1) = ROTMAT(I,1)<br />

ROWMAT(2) = ROTMAT(I,2)<br />

ROWMAT(3) = ROTMAT(I,3)<br />

CALL GFANG(ROWMAT,COSTH,SINTH,COSPH,SINPH,ROTATE)<br />

THETA(I) = ATAN2(SINTH,COSTH)<br />

PHI(I) = ATAN2(SINPH,COSPH)<br />

10 CONTINUE<br />

.<br />

. {\sl Transform to degrees}<br />

.<br />

CALL GSROTM(IROT,THETA(1),PHI(1),THETA(2),PHI(2),THETA(3),PHI(3))<br />

CALL GPROTM<br />

(IROT)<br />

Prints the rotation matrix elements and angles.<br />

IROT<br />

(INTEGER) rotation matrix number: if IROT=0 all rotation matrixes will be printed, if IROT

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!