c
c  ____________________________
c
c
	subroutine rod2eul(d1,d2,d3,rodr)
	real d1,d2,d3,rodr(3),sum,diff
c  converts Rodrigues vector to Bunge Euler angles 
c  d1== phi1,  d2== Phi,  d3== phi2
c
	sum=atan(rodr(3))
	diff=atan(rodr(2)/rodr(1))
	d1=sum+diff
	d2=2.*atan(rodr(1)*cos(sum)/cos(diff))
	d3=sum-diff
c
	return
	end
c
c
c  ____________________________
c