#  表題  回転球殻モデルにおける MHD ダイナモシミュレーション(松井宏晃)
#
#  履歴  1999/08/25  小高正嗣

■モデル概要

  球殻対流モデル
  ・球面上の MHD 方程式
    - 渦度方程式
    - 温度の式
    - 磁場誘導方程式
  ・内核と外核の半径比は 2:5.

  境界条件
  ・ICB, CMB 上で rigid-rid (回転系)
  ・磁場境界; 球殻内部, 外部は絶縁体
  ・温度境界; 球殻内部熱源なし, 境界で一様固定温度.

  パラメータ([]内は実際の地球)
  ・Ra(レイリー数)	1.5E+4〜1.8Rac		[1.0E+30]
  ・Pr(プランドル数)	1			[0.1]
  ・Ta(テイラー数)	2.5E+5			[1.0E+30]
  ・Pm(磁気プランドル)	15			[1.0E-6]

  解法
  ・角度方向は球面調和関数展開, P_n^m の n=18
    (非線形項も波数空間で解く)
  ・動径方向は 2 次中心差分, 境界を含め 64点

■結果

  ・Pm を小さくすると(磁気拡散を大きくすると), 磁場の強さに違いが現れる.
    磁場が大きく成長する解と, 周期的に成長・減衰する.