If <tt>constrainedMD = 1</tt>, then during MD simulation, the scaled coordinates of the
atoms specified in array <tt>constrainatoms</tt> will be constrained in the hyperplane 
perpendicular to <tt>SR2 - SR1</tt>.
