Global existence and convergence of a flow to Kazdan-Warner equation with non-negative prescribed function

We consider an evolution problem associated to the Kazdan-Warner equation on a closed Riemann surface $(Sigma,g)$ begin{align*} -Delta_{g}u=8pileft(frac{he^{u}}{int_{Sigma}he^{u}{rm d}mu_{g}}-frac{1}{int_{Sigma}{rm d}mu_{g}}right) end{align*} where the prescribed function $hgeq0$ and $max_{Sigma}h>0$. We prove the global existence and convergence under additional assumptions such as begin{align*} Delta_{g}ln h(p_0)+8pi-2K(p_0)>0 end{align*} for any maximum point $p_0$ of the sum of $2ln h$ and the regular part of the Green function, where $K$ is the Gaussian curvature of $Sigma$. In particular, this gives a new proof of the existence result by Yang and Zhu [Proc. Amer. Math. Soc. 145 (2017), no. 9, 3953-3959] which generalizes existence result of Ding, Jost, Li and Wang [Asian J. Math. 1 (1997), no. 2, 230-248] to the non-negative prescribed function case.