In this paper, an adaptive numerical method is proposed to solve the Gierer-Meinhardt (GM) system on irregular domain. The method works for domains defined by level sets of implicit functions and the generated mesh is of high quality. The method is shown to be effective by comparing with asymptotic result. Boundary spike solutions of the GM system are obtained and studied numerically, including stability of boundary spike and spike motion along the boundary.