Background: Unbound single-particle states become important in determining the properties of a hot nucleus as its temperature increases. We present relativistic mean field (RMF) for hot nuclei considering not only the self-consistent temperature and density dependence of the self-consistent relativistic mean fields but also the vapor phase that takes into account the unbound nucleon states.Purpose: The temperature dependence of the pairing gaps, nuclear deformation, radii, binding energies, entropy, and caloric curves of spherical and deformed nuclei are obtained in self-consistent RMF calculations up to the limit of existence of the nucleus.Method: We perform Dirac-Hartree-Bogoliubov (DHB) calculations for hot nuclei using a zero-range approximation to the relativistic pairing interaction to calculate proton-proton and neutron-neutron pairing energies and gaps. A vapor subtraction procedure is used to account for unbound states and to remove long range Coulomb repulsion between the hot nucleus and the gas as well as the contribution of the external nucleon gas.Results: We show that $p\ensuremath{-}p$ and $n\ensuremath{-}n$ pairing gaps in the $^{1}\mathrm{S}_{0}$ channel vanish for low critical temperatures in the range ${T}_{{c}^{p}}\ensuremath{\approx}0.6--1.1\phantom{\rule{0.28em}{0ex}}\mathrm{MeV}$ for spherical nuclei such as $^{90}\mathrm{Zr}$, $^{124}\mathrm{Sn}$, and $^{140}\mathrm{Ce}$ and for both deformed nuclei $^{150}\mathrm{Sm}$ and $^{168}\mathrm{Er}$. We found that superconducting phase transition occurs at ${T}_{{c}^{p}}=1.03{\mathrm{\ensuremath{\Delta}}}_{pp}(0)$ for $^{90}\mathrm{Zr}$, ${T}_{{c}^{p}}=1.16{\mathrm{\ensuremath{\Delta}}}_{pp}(0)$ for $^{140}\mathrm{Ce}$, ${T}_{{c}^{p}}=0.92{\mathrm{\ensuremath{\Delta}}}_{pp}(0)$ for $^{150}\mathrm{Sm}$, and ${T}_{{c}^{p}}=0.97{\mathrm{\ensuremath{\Delta}}}_{pp}(0)$ for $^{168}\mathrm{Er}$. The superfluidity phase transition occurs at ${T}_{{c}^{p}}=0.72{\mathrm{\ensuremath{\Delta}}}_{nn}(0)$ for $^{124}\mathrm{Sn}$, ${T}_{{c}^{p}}=1.22{\mathrm{\ensuremath{\Delta}}}_{nn}(0)$ for $^{150}\mathrm{Sm}$, and ${T}_{{c}^{p}}=1.13{\mathrm{\ensuremath{\Delta}}}_{nn}(0)$ for $^{168}\mathrm{Er}$. Thus, the nuclear superfluidity phase---at least for this channel---can only survive at very low nuclear temperatures and this phase transition (when the neutron gap vanishes) always occurs before the superconducting one, where the proton gap is zero. For deformed nuclei the nuclear deformation disappear at temperatures of about ${T}_{{c}^{s}}=2.0--4.0\phantom{\rule{0.28em}{0ex}}\mathrm{MeV}$, well above the critical temperatures for pairing, ${T}_{{c}^{p}}$. If we associate the melting of hot nuclei into the surrounding vapor with the liquid-gas phase transition our results indicate that it occurs at temperatures around $T=8.0--10.0\phantom{\rule{0.28em}{0ex}}\mathrm{MeV}$, somewhat higher than observed in many experimental results.Conclusions: The change of the pairing fields with the temperature is important and must be taken into account in order to define the superfluidity and superconducting phase transitions. We obtain a Hamiltonian form of the pairing field calibrated by an overall constant ${c}_{\mathrm{pair}}$ to compensate for deficiencies of the interaction parameters and of the numerical calculation. When the pairing is not zero, the states close to the Fermi energy make the principal contribution to the anomalous density that appears in the pairing field. By including temperature through the use of the Matsubara formalism, the normal and anomalous densities are multiplied by a Fermi occupation factor. This leads to a reduction in the anomalous density and in the pairing as the temperature increases. When the temperature increases $(T\ensuremath{\ge}4\phantom{\rule{0.28em}{0ex}}\mathrm{MeV})$, the effects of the vapor phase that take into account the unbound nucleon states become important, allowing the study of nuclear properties of finite nuclei from zero to high temperatures.