We propose a new approach, based on the puncture method, to construct black hole initial data in the so-called trumpet geometry, i.e. on slices that asymptote to a limiting surface of non-zero areal radius. Our approach is easy to implement numerically and, at least for non-spinning black holes, does not require any internal boundary conditions. We present numerical results, obtained with a uniform-grid finite-difference code, for boosted black holes and binary black holes. We also comment on generalizations of this method for spinning black holes.PACS numbers: 04.20. Ex, 04.25.dg, 04.70.Bw Numerical simulations of black hole spacetimes have recently experienced a dramatic breakthrough (see [1,2,3] as well as numerous later publications). Most of these simulations now adopt some variation of the BSSN formulation [4,5] together with the moving puncture [2,3] method to handle the black hole singularities.The moving-puncture method is based on a set of empirically found coordinate conditions, namely the "1+log" slicing condition for the lapse [6] and a "Γ-freezing" gauge condition for the shift [7]. As demonstrated by [8,9,10,11], dynamical simulations of a Schwarzschild spacetime using these coordinate conditions settle down to a spatial slice that terminates at a non-zero areal radius, and hence does not encounter the spacetime singularity at the center of the black hole. An embedding diagram of such a slice, which suggests the name trumpet data, is shown in Fig. 2 of [11].Typically, moving-puncture simulations adopt initial data that are constructed using the puncture method [12,13,14]. As we explain in more detail below, the central idea of the puncture method is to write the conformal factor as a sum of an analytically known, singular background term, and a correction term that is unknown but regular. The equations for the correction term can then be solved everywhere, without any need for excision or any other means of dealing with the black hole singularity. To date, all applications of this method have adopted Schwarzschild data on a slice of constant Schwarzschild time as the background solution. These data connect spatial infinity in one universe with spatial infinity in another universe; the resulting initial data therefore represent wormhole data.Clearly, it would be desirable to produce initial data that represent black holes as trumpets rather than wormholes, since otherwise moving-puncture evolutions will drive the individual black holes to a trumpet geometry. Problems with one possible approach, based on a stationary 1+log slicing in the context of the conformal thin-sandwich decomposition, were described in [15]. Recently, trumpet initial data on hyperboloidal slices using * Also at Department of Physics, University of Illinois at UrbanaChampaign, Urbana, IL 61801 an excision method were constructed in [16]. In this paper we demonstrate how such data can be produced by generalizing the puncture method, which does not require any excision or internal boundary conditions, and is very easy to imple...