Let $f(x)$ be a harmonic function in $\mathbb{R}^n$, and let $Z$ be the zero level set $$ Z = \{(x)\in\mathbb{R}^n \mid f(x) = 0\}. $$ Arguing using the Harnack inequality, I believe one can show that $$ \frac{|\nabla f(x)|}{|f(x)|} \lesssim r^{-1} $$ when $d(x,Z) = \inf_{z\in Z} |x-z| \geq r$. I wonder if the opposite inequality holds.
Do there exist constants $c,C>0$ depending only on $n$ such that $$ c\,d(x,Z)^{-1} \leq \frac{|\nabla f(x)|}{|f(x)|} \leq C\, d(x,Z)^{-1}? $$
This sounds like it is too good to be true, so I would be interested to see a counterexample. I am also curious about the special case that $n=2$ and $f(x)$ is the real part of a polynomial.