  We prove (sub)mean value formulas at the point $0\in\Sigma$ for (sub)harmonic
functions a on a hypersurface $\Sigma\subset\mathbb{R}^{n+1}$ where the
differentiable structure and the surface measure depend on the ambient Grushin
structure.
