One of the cornerstone effects in spintronics is spin pumping by dynamical magnetization that is steadily precessing (around, e.g., the $z$-axis) with frequency $\omega_0$, due to absorption of low-power microwaves of frequency $\omega_0$ under the resonance conditions and in the {\em absence} of any applied bias voltage. The two-decades-old ``standard model'' of this effect, based on the scattering theory of adiabatic quantum pumping, predicts that component $I^{S_z}$ of spin current vector $\big( I^{S_x}(t),I^{S_y}(t),I^{S_z} \big) \propto \omega_0$ is time-independent while $I^{S_x}(t)$ and $I^{S_y}(t)$ oscillate harmonically in time with a {\em single} frequency $\omega_0$; whereas pumped charge current is zero $I \equiv 0$ in the same adiabatic $\propto \omega_0$ limit. Here we employ more general than ``standard model'' approaches, time-dependent nonequilibrium Green's function (NEGF) and Floquet-NEGF, to predict unforeseen features of spin pumping---precessing localized magnetic moments within ferromagnetic metal (FM) or antiferromagnetic metal (AFM), whose conduction electrons are exposed to spin-orbit coupling (SOC) of either intrinsic or proximity origin, will pump both spin $I^{S_\alpha}(t)$ and charge $I(t)$ currents. All four of these functions harmonically oscillate in time at {\em both even an odd integer multiples} $N\omega_0$ of the driving frequency $\omega_0$. The cutoff order of such high-harmonics increases with SOC strength, reaching $N_\mathrm{max} \simeq 11$ in the chosen-for-demonstration one-dimensional FM or AFM models. Higher cutoff $N_\mathrm{max} \simeq 25$ can be achieved in realistic two-dimensional (2D) FM models defined on the honeycomb lattice, where we provide prescription on how to realize them using 2D magnets and their heterostructures.