理論磁気異常計算(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単位)を得る.