contErr is based on rhoPhi not on phi. The term contErr/rho*p in TEq makes a large correction which occasionally makes T unbounded
