c
c  ____________________________
c
c
	subroutine eul2rod(d1,d2,d3,rodr)
	real d1,d2,d3,rodr(3),sum,diff,csum,cdiff,sdiff
c  converts Bunge EUler angles to Rodrigues vector
c  d1== phi1,  d2== Phi,  d3== phi2
c
	sum=(d1+d3)/2.
	diff=(d1-d3)/2.
	csum=cos(sum)
	cdiff=cos(sum)
	sdiff=sin(diff)
	t2=tan(d2/2.)
c
	rodr(1)=t2*cdiff/csum
	rodr(2)=t2*sdiff/csum
	rodr(3)=tan(sum)
c
	return
	end
c
c
c  ____________________________
c