real(kind=dp)::qbound_spin=-1!! spin mixing qbound in degrees (if rms of spin directions is < qbound all angles are fixed, defaut negative value means no automatic fixing)
real(kind=dp)::angles_cutoff=1.0e-4_dp!! MdSD: threshold for the magnitude of the spin moment below which the angles stop being updated in linear mixing
real(kind=dp)::angles_cutoff=-1._dp!! MdSD: threshold for the magnitude of the spin moment below which the angles stop being updated in linear mixing