Added source code from crystallization to main trunk of plumed

The source code from crystallization allows you to calculate
Q4 and Q6 parameters.  You can also calculate local Q4 and local
Q6 parameters and you create compound collective variables based
on the relative orientations of internal vectors in adjacent molecules.

Regtests for this functionality are reasonably thorough.  Manual is
not so great as yet