In this paper, we consider the existence and limiting behaviour of solutions to a semilinear elliptic equation arising from confined plasma problem in dimension two\[ \begin{cases} -\Delta u=\lambda k(x)f(u) & \text{in}\ D,\\ u= c & \displaystyle\text{on}\ \partial D,\\ \displaystyle - \int_{\partial D} \frac{\partial u}{\partial \nu}\,{\rm d}s=I, \end{cases} \] where $D\subseteq \mathbb {R}^2$ is a smooth bounded domain, $\nu$ is the outward unit normal to the boundary $\partial D$, $\lambda$ and $I$ are given constants and $c$ is an unknown constant. Under some assumptions on $f$ and $k$, we prove that there exists a family of solutions concentrating near strict local minimum points of $\Gamma (x)=({1}/{2})h(x,\,x)- ({1}/{8\pi })\ln k(x)$ as $\lambda \to +\infty$. Here $h(x,\,x)$ is the Robin function of $-\Delta$ in $D$. The prescribed functions $f$ and $k$ can be very general. The result is proved by regarding $k$ as a $measure$ and using the vorticity method, that is, solving a maximization problem for vorticity and analysing the asymptotic behaviour of maximizers. Existence of solutions concentrating near several points is also obtained.