Minor fix in calculation s_crit in ViscoelasticPM.