The geometric Neumann problem for the Liouville equation

Let Ω denote the upper half-plane ${\mathbb{R}_+^2}$ or the upper half-disk ${D_{\varepsilon}^+\subset \mathbb{R}_+^2}$ of center 0 and radius ${\varepsilon}$ . In this paper we classify the solutions ${v\inC^2(\overline{\Omega}\setminus\{0\})}$ to the Neumann problem $$\left\{\begin{array}{lll}{\Delta v+2 Ke^v=0\quad {\rm in}\,\Omega\subseteq \mathbb{R}^2_+=\{(s, t)\in \mathbb{R}^2: t >0 \},}\\ {\frac{\partial v}{\partial t}=c_1e^{v/2}\quad\quad\quad{\rm on}\,\partial\Omega\cap\{s >0 \},}\\ {\frac{\partial v}{\partial t}=c_2e^{v/2}\quad\quad\quad{\rm on}\,\partial\Omega\cap\{s <0 \},}\end{array}\right.$$ where ${K, c_1, c_2 \in \mathbb{R}}$ , with the finite energy condition ${\int_{\Omega} e^v < \infty}$ As a result, we classify the conformal Riemannian metrics of constant curvature and finite area on a half-plane that have a finite number of boundary singularities, not assumed a priori to be conical, and constant geodesic curvature along each boundary arc.