The Gauss-Bonnet formula for harmonic surfaces

We consider harmonic immersions in R-d of compact Riemann surfaces with finitely many punctures where the harmonic coordinate functions are given as real parts of meromorphic functions. We prove that such surfaces have finite total Gauss curvature. The contribution of each end is a multiple of 2 pi, determined by the maximal pole order of the meromorphic functions. This generalizes the well known Gackstatter-Jorge-Meeks formula for minimal surfaces. The situation is complicated as the ends with their induced metrics are generally not conformally equivalent to punctured disks, nor do the surfaces generally have limit tangent planes at the ends.