A self-consistent model of a spherical microwave discharge is presented. The model takes into account thermal and ionization non-equilibrium of the discharge plasma. We adopt a partial local thermodynamic equilibrium model for the plasma in two-temperature approximation. Numerical experiments are carried out for the discharge in argon at atmospheric pressure. Results are presented for the characteristics of the discharge plasma against the external parameters (the power and frequency of the applied electromagnetic field and the size of the discharge chamber). Presented model results are compared with the results for the ionization equilibrium model of the spherical microwave discharge.