Update of axion CDM energy density

We improve the estimate of the axion CDM energy density by considering the new values of current quark masses, the QCD phase transition effect and a possible anharmonic effect.