Phase diagram and conformal string excitations of square ice using gauge invariant matrix product states