We have constructed a non-spherical model for the hot oxygen corona of Mars by including the effects of planetary rotation and diurnal variation of the Martian ionosphere. Exospheric oxygen densities are calculated by integrating ensemble of ballistic and escaping oxygen atoms from the exobase over the entire planet. The hot oxygen atoms are produced by dissociative recombination of $O^+_2$, the major ion in the Martian ionosphere. The densities of hot oxygen atoms at the exobase are estimated from electron densities which have been measured to vary with solar zenith angle. Our model shows that the density difference of hot oxygen atoms between noon and terminator is about two orders of magnitude near the exobase, but reduces abruptly around altitudes of 2000 km due to lateral transport. The diurnal variation of hot oxygen densities remains significant up to the altitude of 10000 km. The diurnal variation of the hot oxygen corona should thus be considered when the upcoming Nozomi measurements are analyzed. The non-spherical model of the hot oxy-gen corona may contribute to building sophisticate solar wind interaction models and thus result in more accurate escaping rate of oxygens from Mars.