Add function to check covariances for consistency

Still to be done: if covariances are specified, check M_.Sigma_e and M_.H by calling this function after setting those matrices using the preprocessor.