In Introduction,  reasons why we choose FEM are added. Surely, 
FVM requires less memory than FEM, but any types of unstructured 
grid can be set using FEM.

   In Section 2, some papers are referred for the mesh generation 
on a sphere surface, and the process to obtain the mesh data for 
the local domains is described. This process is performed by the 
petitioner provides by GeoFEM. The local mesh data with 
communication table are generated by the partitioner.

   Discussion for performances of the present simulation and the 
spectral method is added.  Significantly high efficiency is utilized 
in the case of the spectral method because of the architecture of 
an SX-4 and methods of parallelization. However, the number of 
processor is controlled  by the number of radial grid points in this 
scheme. Furthermore, one of methods for increasing the 
performance of GeoFEM on an SR8000 is described.

   As a future plan, we describe that a MHD simulation will be 
performed on the Earth Simulator, which has 512 nodes of SMP 
type 8 vector processor.

   Finally, a lot of grammatical errors are corrected throughout this paper.