c
c ________________________________________
c
	subroutine charlow(axis,aenerg,amobility,numtypes)
c
	integer numtypes
	real axis(20,3),aenerg(20),amobility(20),rnorm,remax,rmmax
c  AXIS describes the rotation axes, AENERG the energy
c  AMOBILITY the mobility
c
	include 'common.f'
c
	numtypes=13
c
	do 10, i=1,20
		aenerg(i)=0.
		amobility(i)=0.
		do 10, j=1,3
			axis(i,j)=0.
10			continue
c
	axis(1,1)=0.
	axis(1,2)=0.
	axis(1,3)=1.
	aenerg(1)=0.33
	amobility(1)=0.01
c
	axis(2,1)=0.
	axis(2,2)=1.
	axis(2,3)=5.
	aenerg(2)=0.3
	amobility(2)=0.02
c
	axis(3,1)=0.
	axis(3,2)=2.
	axis(3,3)=5.
	aenerg(3)=0.3
	amobility(3)=0.1
c
	axis(4,1)=0.
	axis(4,2)=2.
	axis(4,3)=3.
	aenerg(4)=0.26
	amobility(4)=0.02
c
	axis(5,1)=0.
	axis(5,2)=1.
	axis(5,3)=1.
	aenerg(5)=0.3
	amobility(5)=0.17
c
	axis(6,1)=2.
	axis(6,2)=7.
	axis(6,3)=7.
	aenerg(6)=0.26
	amobility(6)=0.1
c
	axis(7,1)=2.
	axis(7,2)=3.
	axis(7,3)=3.
	aenerg(7)=0.23
	amobility(7)=0.05
c
	axis(8,1)=1.
	axis(8,2)=1.
	axis(8,3)=1.
	aenerg(8)=0.23
	amobility(8)=0.95
c
	axis(9,1)=3.
	axis(9,2)=3.
	axis(9,3)=5.
	aenerg(9)=0.23
	amobility(9)=0.02
c
	axis(10,1)=1.
	axis(10,2)=1.
	axis(10,3)=3.
	aenerg(10)=0.23
	amobility(10)=0.02
c
	axis(11,1)=1.
	axis(11,2)=1.
	axis(11,3)=7.
	aenerg(11)=0.33
	amobility(11)=0.05
c
	axis(12,1)=1.
	axis(12,2)=2.
	axis(12,3)=5.
	aenerg(12)=0.3
	amobility(12)=0.1
c
	axis(13,1)=4.
	axis(13,2)=8.
	axis(13,3)=11.
	aenerg(13)=0.26
	amobility(13)=0.17
c
	do 100, i=1,numtypes
		rnorm=0.
		do 90, j=1,3
			rnorm=rnorm+axis(i,j)**2
90		continue
		rnorm=sqrt(rnorm)
		do 95, j=1,3
			axis(i,j)=axis(i,j)/rnorm
95		continue
100	continue
c  make unit vectors
c
		remax=aenerg(1)
		rmmax=amobility(1)
	do 200, i=2,numtypes
		if(aenerg(i).gt.remax) remax=aenerg(i)
		if(amobility(i).gt.rmmax) rmmax=amobility(i)
200		continue
	do 210, i=1,13
		aenerg(i)=aenerg(i)/remax
		amobility(i)=amobility(i)/rmmax
210		continue
c  normalize the energies and mobilities
c
	return
	end
c
c
c ________________________________________
c