理論磁気異常計算(CALMA)


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

2.呼出し手順と引数

      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単位)を得る.