Assuming that your mesh is fine and solution has converged sufficiently, your Cd will depend on those reference values that you have defined. I would imagine that the experimental results that you mention would have their corresponding values documented.