  A formulation for thermal noise in the stochastic form of the Landau Lifshitz
Bloch equation used for modeling the magnetization dynamics at elevated
temperatures is presented. The diffusion coefficients for thermal fluctuations
are obtained via the Fokker-Plank equation using the mean field approximation
of the field in defining the free energy. The presented model leads to a mean
magnetization consistent with the equilibrium magnetization for small and large
particles. The distribution of the magnetization magnitude is of a Poisson-like
type rather than the classical Boltzmann distribution. The presented model was
tested by studying the equilibrium magnetization in macrospin particles at high
temperatures. The model is appealing for multi-scale modeling, such as modeling
heat assisted magnetic recording systems and all-optical magnetization reversal
