Inclusive photoproduction of $\dspm$ in ep collisions at HERA has been measured with the ZEUS detector for photon-proton centre of mass energies in the range \linebreak \wrang and photon virtuality Q~2 < 4 \g2. The cross section $\sigma_{ep \to \ds X} $ integrated over the kinematic region \ptrangand \etarang is {\xsecs}. Differential cross sections as functions of $p_{\perp}~{\ds}$, $\eta~{\ds}$ and W are given. The data are compared with two next-to-leading order perturbative QCD predictions. For a calculation using a massive charm scheme the predicted cross sections are smaller than the measured ones. A recent calculation using a massless charm scheme is in agreement with the data.