Harnack inequality for harmonic functions relative to a nonlinear $p$ homogeneous Riemannian Dirichlet form