A perturbation of spacetime Laplacian equation

We study a perturbation begin{equation} Delta u + P | abla u| = h | abla 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.