A finite temperature hydrodynamic model is derived for the spin-1 ultracold bosons by the many-particle quantum hydrodynamic method. It is presented as the two fluid model of the Bose-Einstein condensate (BEC) and normal fluid. The continuity, Euler, spin evolution, and nematic tensor evolution equations are derived for each fluid. The linear and quadratic Zeeman effects are included. Scalar and spin-spin like short-range interactions are considered in the first order by the interaction radius. Obtained hydrodynamic equations are also represented as the set of two nonlinear Pauli equations. The spectrum of the bulk collective excitations is considered for the ferromagnetic phase in the small temperature limit. The spin wave is not affected by the presence of the small temperature in the described minimal coupling model, where the thermal part of the spin-current of the normal fluid is neglected. The two sound waves are affected by the spin evolution in the same way as the change of spectrum of the single sound wave in BEC, where speed of sound is proportional to g1+g2 with gi as the interaction constants. © 2021 Author(s).