１．機　能
  ・magafd ルーチンは，磁化ソースの磁化の強度と伏角・偏角および地球磁場の伏角・
    偏角を設定する．
  ・mpoint／mvline／mhrect／mprism の各ルーチンは，計算するソースのモデルを，点
    ／鉛直な線／水平な矩形の面／直方体 に設定し，それぞれのモデルパラメータを設
    定する．
  ・calma ルーチン（関数）は，観測点の位置座標を指定して，設定済みのモデルの磁
    気異常を計算し，その結果を関数値とする．
 (注) 本ルーチン群の中では，北向きを X，東向きを Y，下向きを Z とする右手系のデ
  カルト座標を想定しており，すべての座標（長さ）について同一の単位を使用する限
  り，座標の単位は任意である．

２．呼出し手順と引数

      call magafd(am, the, phi, dip, dec)
  am        [float]  磁化強度を A/m単位で与える．
  the, phi  [float]  磁化方向の伏角と偏角を度単位で与える．
  dip, dec  [float]  地球磁場方向の伏角と偏角を度単位で与える．

      call mpoint(xc, yc, zc, vol)
      call mvline(xc, yc, zs, w, area)
      call mhrect(xs, xt, ys, yt, zc, h)
      call mprism(xs, xt, ys, yt, zs, w)
  xc, yc    [float]  ソース中心の水平位置座標を与える．
  zc        [float]  ソース中心の深度座標を与える．
  xs, xt    [float]  ソースの X座標方向のひろがりの下限と上限を与える．
  ys, yt    [float]  ソースの Y座標方向のひろがりの下限と上限を与える．
  zs, w     [float]  ソースの深度方向のひろがりの下限（上面深度）と厚さ（線の
                     長さ）を与える．但し，w＝0. のとき底なしモデル（w→∞）と
                     みなす．
  vol       [float]  点で代表させる立体の体積を与える．
                     （単位は，座標の単位の 3乗）
  area      [float]  線で代表させる立体の断面積を与える．
                     （単位は，座標の単位の 2乗）
  h         [float]  面で代表させる立体の厚さを与える．
                     （単位は，座標の単位と同じ）

      f = calma(xp, yp, zp)
  xp, yp    [float]  計算点の水平位置座標を与える．
  zp        [float]  計算点の深度座標を与える．
  f         [float]  計算結果．全磁力異常の値（nT単位）を得る．