The structure of geomagnetism is very complex, and there are still some problems left in magnetohydrodynamics (MHD) simulations. Recently, the macro spin model has been suggested. This is the idea that geomagnetism is described by interaction with many local dynamo elements (macro-spins). This model can reproduce some features of geomagnetism by solving equations of motion with only one variable. In this paper, we study this model to make more general by adding one variable. In this result, our model becomes possible for several things that are not treated in the previous study, for example, migration of the North (or South) Magnetic Pole, comparison with observed data of magnetic field distributions expressed in two directions, etc. Moreover, as a result of an application to the sun, we could reproduce periodic polarity reversals and the power index of the power spectrum, etc. In addition, we investigate the statistical properties of the pole migration.