We calculate charmonium correlators on the lattice with 2+1-flavors of sea quarks and charm valence quark both described by the Möbius domain-wall fermion. Temporal moments of the correlators are calculated and matched to perturbative QCD formulae to extract the charm quark mass m c (µ) and strong coupling constant α s (µ). Lattice data at three lattice spacings, 0.044, 0.055, and 0.080 fm, are extrapolated to the continuum limit. The correlators in the vector channel are confirmed to be consistent with the experimental data for e + e − → cc, while the pseudo-scalar channel is used to extract m c (µ) and α s (µ). We obtain m c (3 GeV) = 1.003(10) GeV and α MS(4) s (3 GeV) = 0.253(13). Dominant source of the error is the truncation of perturbative expansion at α 3