An important aspect of the proposed approach is a direct relation between low-level calculations and the finite element Monte Carlo magnetic simulations. The quantities describing the re-scaled system, e.g., spin vector, exchange integral parameter and magnetic anisotropy attributed to the super-cell, are derived from the unscaled one. In a situation when the nodes in the unscaled system represent atoms in a crystal structure, the mentioned quantities can be determined using, for example, density functional theory (DFT) based calculations.