changed real equality statements to ULP*spacing(X)-type statements, based upon warnings from gfortran 4.8.1