Repeating this hundreds of thousands of times generates a sequence of configurations that are in thermodynamic equilibrium for that temperature. Then any kind of averages for that temperaure can be determined, such as the total magnetization.
Try out the xmc visualization program found in the zipped file at this link xmc. (You will need a c-compiler, unix system and X11 libraries.)
Try out the xmd simulation program found in the zipped file at this link xmd. (You will need a c-compiler, unix system and X11 libraries.)