c
c ------------------------------------------------
c
	subroutine rod2ang(r1,angle)
	real r1(3),angle
c  extracts a rotation angle from a Rodrigues vector
c
	angle=2.*atan((sqrt(r1(1)**2+r1(2)**2+r1(3)**2)))
	return
	end
c
c
c ------------------------------------------------
c