u^{m-0}_h)^0 \, \rho ^{m-0}\, | \nabla \phi ^{m-0}|\nonumber \\&\qquad + \tau _m \int _\varOmega |\nabla u^m_h|^0 \, \rho ^m\, |\nabla \phi ^{m}|+ \gamma \tau _m^0 \int _\varOmega I^m_h \tilde{\rho }^m | \nabla u^m_h|^0 \nonumber \\&\q