c
c ________________________________________
c
	subroutine vecnorm(vec)
	real vec(3),rnorm
	rnorm=0.
	do 10, i=1,3
	rnorm=rnorm+vec(i)**2
10	continue
	if(rnorm.le.0.0) return
	rnorm=sqrt(rnorm)
	do 20, i=1,3
	vec(i)=vec(i)/rnorm
20	continue
	return
	end
c
c ________________________________________
c