Using the phenomenological relativistic harmonic model (RHM) for quarks, we have obtained the masses of S wave charmonium states. The full Hamiltonian used in the investigation has Lorentz scalar plus vector confinement potential, along with the confined one gluon exchange potential (COGEP). A good agreement with the experimental masses for the ground state and the radially excited states is obtained for both the triplet and singlet S wave mesons. The calculated charge radii, meson decay constants, leptonic decay width, two photon decay width, and the radiative M1 decay width are in good agreement with the experimental results.