We present an updated global fit of neutrino oscillation data in the simplest three-neutrino framework. In the present study we include up-to-date analyses from a number of experiments. Concerning the atmospheric and solar sectors, besides the data considered previously, we give updated analyses of IceCube DeepCore and Sudbury Neutrino Observatory data, respectively. We have also included the latest electron antineutrino data collected by the Daya Bay and RENO reactor experiments, and the long-baseline T2K and NOνA measurements, as reported in the Neutrino 2020 conference. All in all, these new analyses result in more accurate measurements of θ13, θ12, Delta {m}_{21}^2 and left|Delta {m}_{31}^2right| . The best fit value for the atmospheric angle θ23 lies in the second octant, but first octant solutions remain allowed at ∼ 2.4σ. Regarding CP violation measurements, the preferred value of δ we obtain is 1.08π (1.58π) for normal (inverted) neutrino mass ordering. The global analysis still prefers normal neutrino mass ordering with 2.5σ statistical significance. This preference is milder than the one found in previous global analyses. These new results should be regarded as robust due to the agreement found between our Bayesian and frequentist approaches. Taking into account only oscillation data, there is a weak/moderate preference for the normal neutrino mass ordering of 2.00σ. While adding neutrinoless double beta decay from the latest Gerda, CUORE and KamLAND-Zen results barely modifies this picture, cosmological measurements raise the preference to 2.68σ within a conservative approach. A more aggressive data set combination of cosmological observations leads to a similar preference for normal with respect to inverted mass ordering, namely 2.70σ. This very same cosmological data set provides 2σ upper limits on the total neutrino mass corresponding to Σmν< 0.12 (0.15) eV in the normal (inverted) neutrino mass ordering scenario. The bounds on the neutrino mixing parameters and masses presented in this up-to-date global fit analysis include all currently available neutrino physics inputs.