Maybe the definition of hTinv should be

hTinv = CellField(1 ./ collect(get_array(∫(1)dΩp) .^ (1/2)), Ωp)

i.e. with a power 1/2, since it is a 2D problem.