c	TEXTUR - sets up required arrays
c	CHARLOW - creates properties for low angle bdies
c	QUATB  convert Bunge euler angles to a quaternion
c	MISQUAT  combines two quaternions and finds the minimum angle
c	MISQUATINV  combines two quaternions (reversed order) and finds the minimum angle
c	CSLQUAT  gets distance (angle) from a CSL boundary: misorient must be in 
c		the fundamental zone,  q1<q2<q3<q4
c	Q2EULB       converts quaternion to Bunge Euler angles
c	Q2MAT  converts quaternion to a 3x3 matrix
c	QRVEC  uses a quaternion to rotate a vector
c	COMQUAT  combines two quaternions (no disorientation)
c	PRESYMM  applies symmetry element from QUATSYMM as 1st q
c	POSTSYMM applies symm elem as 2nd q
c	DISQUAT  finds the disorientation the hard way but also determines
c		which symmetry operators put you in the fundamental zone
c	AXISANG2ROT(axis,ang,rot)
c  		AXIS  is the normalized axis vector
c  		ANG is the rotation angle in radians
c  		ROT is the 3x3 matrix
c	SCALAR    returns scalar product of vectors A and B
c	DONORM   normalizes a vector
c	AXVECB2VECC   transforms vector B by A to get vector C
c	VECNEG	takes negative of vector
c	RANEULER    picks a set of random Euler angles
c	EULER  converts Euler angles to/from matrix (active rotation)
c  	M2R   converts a MISOR matrix to Rodrigues vector
c	EUL2ROD  converts Euler angles to Rodrigues vector
c	ROD2EUL  converts a Rodrigues vector to Euler angles
c	SYMM_APPLY   postmultiplies by transpose of the symmetry operator
c	SYMM_PRE_APPLY  premultiplies by symmetry operator
c	DIS2MAT  
c  		analyzes EIN as a rotation matrix (assumed!) to find
c  		the disorientation and return the corresponding matrix, EOUT,
c  		the rotation ANGLE (degrees) and the AXIS of rotation (in the SST)
c  		INDEX1 is the index to the first symm oper, INDEX2, the second
c  		INVERT tells whether the rotations were taken in the opposite order
c  		ROUT is the Rodrigues vector of the result
c	READSIGMAS  reads in special boundary parameters
c	RODCOM  combines two Rodrigues vectors
c	RODCOMABT.f   
c  		combines two rotations as Rodrigues vectors
c  		ROUT= R1+R2- R1 cross R2 / 1 - (R1dotR2)
c  		except use the transpose of second rotation
c	RODCOMATB
c  		combines two rotations as Rodrigues vectors
c  		ROUT= R1+R2- R1 cross R2 / 1 - (R1dotR2)
c  		except use the transpose of first rotation
c	ROD2ANG  extracts a rotation angle from a Rodrigues vector
c	CSL   tests a misorientation for proximity to a CSL value
c	AXB2C   multiplies two matrices (3x3)
c	ATXB2C  multiplies two matrices (3x3) using transpose of 1st
c	AXBT2C  multiplies two matrices (3x3) using transpose of 2nd
c	VECPRO2  vector product of two vectors
c	VECPRO  vector product/ special version from LApp