A perturbation of spacetime Laplacian equation

We study a perturbation \begin{equation} \Delta u + P | \nabla u| = h | \nabla u|, \end{equation} of spacetime Laplacian equation in an initial data set $(M, g, p)$ where $P$ is the trace of the symmetric 2-tensor $p$ and $h$ is a smooth function.