Let $(M,g)$ be a (bad) Riemannian $2$-orbifold with one cone point $x_{1}$ of order $p\geq2$ and underlying space a topological $2$-sphere, i.e., a teardrop with Euler characteristic $\chi(M)=1+\frac{1}{p}$. For $\varepsilon>0$ small, we have length $\operatorname{L}(\partial B_{\varepsilon}(x_{1}))\approx2\pi p^{-1}\varepsilon$ and geodesic curvature $\kappa(\partial B_{\varepsilon }(x_{1}))\approx\varepsilon^{-1}$. The Gauss-Bonnet formula is, using $\chi(M-B_{\varepsilon}(x_{1}))=1$, $$ \int_{M}Rd\mu=\lim_{\varepsilon\rightarrow0}\int_{M-B_{\varepsilon}(x_{1} )}Rd\mu=2\pi+\lim_{\varepsilon\rightarrow0}\int_{\partial B_{\varepsilon }(x_{1})}\kappa ds=2\pi\chi(M). $$
A vector field $X$ is smooth on $M$ if it is smooth on $M-x_{1}$ and its $p$-fold lift via the orbifold chart $\pi:B\rightarrow U$, where $U$ is a neighborhood of $x_{1}$, is smooth (the same applies to functions). This implies that $\lim_{x\rightarrow x_{1}}X\left( x\right) =X\left( x_{1}\right) =0$. Moreover, the divergence theorem is true for $X$: $\int _{M}\operatorname{div}X\,d\mu=0$. In particular, if $u$ is a smooth function on $M$, then $\int_{M}\Delta u\,d\mu=0$.
Conversely, we may solve $\Delta f=h$ if (and only if) $\int_{M}hd\mu=0$. E.g., use the calculus of variations to obtain a minimizer $f\in W^{1,2}(M)$ of the functional $J\left( f\right) =\int_{M}(\frac{1}{2}\left\vert \nabla f\right\vert ^{2}+hf)d\mu$ with the constraint $\int_{M}fd\mu=0$ (so we can apply the Poincar'{e} inequality; the constraint anchors the function). By regularity theory, the minimizer $f$ is smooth and satisfies $-\Delta f+h=\operatorname{const}$, which must be zero. This works on an orbifold with isolated singularites.
Let $r$ be the average scalar curvature. Since $\int_{M}(r-R)d\mu=0$, there exists $f$ (unique up to an additive constant) such that $\Delta f=r-R$, called the curvature potential. Note that $\nabla f\left( x_{1}\right) =0$.
The only closed bad $2$-orbifolds are the teardrop and the spindle (i.e., a bad orbifold with two cone points, of orders $p,q\geq2$, where $p\neq q$ and $\chi=\frac{1}{p}+\frac{1}{q}$). On these, Lang-Fang Wu proved that for any initial smooth metric with positive curvature, the normalized Ricci flow (NRF) asymptotically aproaches a (nonconstant curvature) shrinking gradient Ricci soliton. More specifically, consider the NRF modified by diffeomorphisms, defined by $\frac{\partial}{\partial t}g=-2\left( \operatorname{Ric} +\nabla^{2}f-\frac{r}{2}g\right) $, where the curvature potential $f$ satisfies $\frac{\partial f}{\partial t}=\Delta f-\left\vert \nabla f\right\vert ^{2}+rf$. (The right side of the first equation is actually trace-free, which implies that the area form is pointwise constant under the flow). Then, as $t\rightarrow\infty$, $g\left( t\right) $ and $f\left( t\right) $ converge to a $C^{\infty}$ metric $g_{\infty}$ and a $C^{\infty}$ function $f_{\infty}$, each exponentially fast in all $C^{k}$-norms. They satisfy $\operatorname{Ric}(g_{\infty})+\nabla_{g_{\infty}}^{2}f_{\infty }-\frac{r}{2}g_{\infty}=0$ and $\Delta_{g_{\infty}}f_{\infty}-\left\vert \nabla f_{\infty}\right\vert _{g_{\infty}}^{2}+rf_{\infty}=0$.
So, Wu proved that under the NRF on a bad closed $2$-orbifold there exists a gradient Ricci soliton structure $(\tilde{g}_{\infty},\tilde{f}_{\infty})$ equivalent to $(g_{\infty},f_{\infty})$ such that $g\left( t\right) $ asymptotically approaches $\varphi_{t}^{\ast}\tilde{g}_{\infty}$, where $\{\varphi_{t}\}$ is generated by $\operatorname{grad}_{\tilde{g}_{\infty} }\tilde{f}_{\infty}$. If the initial metric yields a shrinking gradient Ricci soliton structure, then we have equality, i.e., $g\left( t\right) =\varphi_{t}^{\ast}g\left( 0\right) $.
Disclaimer: I haven't been thinking about orbifolds recently, so some of the above is what I think is true from memory without checking the details; please look out for errors.