1      subroutine rota(alfa)
2      implicit double precision (a-h,o-z)
3      common/rotxyz/amat(3,3),bmat(3,3),cmat(3,3),rmat(3,3),todeg
4
5c--   alfa expected in degrees
6      rada = alfa / todeg
7      amat(1,1) =  dcos(rada)
8      amat(1,2) =  dsin(rada)
9      amat(2,1) = -dsin(rada)
10      amat(2,2) =  dcos(rada)
11
12      return
13      end
14