More and more studies indicate that the effects of quantum size and energy level statistics play a crucial role in the thermodynamic properties of ultrasmall metallic grains. This paper aims to investigate how they affect the specific heat of ultrasmall metallic grains in magnetic field. As the particle size decreases, fluctuation effects and the impact of energy level separation are becoming more and more important. The method of static path approximation (SPA) is adopted to handle the fluctuation effect. Random matrix theory (RMT) is adopted due to its successful description of the energy level of metal nanoparticles. The normalized specific heat of several typical temperatures and electron spins were taken in the calculation, and the results were analyzed. It was found that spin and the spin-orbit coupling affect the specific heat very obviously, and the suppressed high spin weakens the contribution of electrons to the heat capacity.