関数Mogi マニュアル

(The documentation of function Mogi)
Last Update: 2021/12/8


◆機能・用途(Purpose)

茂木モデルに基づいて変位または傾斜変動の特定の一成分を計算する。
Compute a specified displacement or a tilt component based on the Mogi model.

◆形式(Format)

#include <tws/Mogi.h>
inline double Mogi
(const char ∗source,const double sx,const double sy,const double sz,
  const double source_intensity,
  const char ∗station,const double x,const double y,
  const char ∗surface,const double z0,
  const double theta_x,const double theta_y, const char ∗unit,
  const char ∗component,const char ∗component_name)


◆引数(Arguments)

この関数では実数の引数を多数用いるため、 指定する位置や順番を間違えるミスが発生しやすい。 それを避けるため、後に続くパラメータが何を表すのかを確認するための 文字列引数を設けている。 例えば"source"という文字列の後にソースに関するパラメータを並べ、 "station"という文字列の後に観測点に関するパラメータを並べる方が パラメータの値だけを書き並べるよりも分かりやすくミスが起きにくい。 これらの文字列はチェックのみに用いられる。
As this function uses many arguments of real numbers, mistakes may easily occur for the location and order of the values. To avoid it, there are several string arguments to confirm what the subsequent parameters represent. For example, it is more comprehensive, and thus mistakes are less frequent, by writing a string "source" before source parameters, and a string "station” before station parameters, than listing only the parameter values. These strings are used only for checks.
source ユーザがパラメータの順番を間違えないための引数であり、 文字列"source"のみ可(他の文字列を指定するとエラーになる)。 この後の引数がソースに関するパラメータであることを示す。
An argument to avoid a mistake of the order of parameters. Only the value "source" is allowed; if another string is given to here, the program finishes as an error. This argument indicates that the subsequent arguments are parameters of the source.
sx 震源位置の\(x\)座標。
The \(x\)-component of the source location.
sy 震源位置の\(y\)座標。
The \(y\)-component of the source location.
sz 震源位置の\(z\)座標。
The \(z\)-component of the source location.
source_intensity 変位や傾斜変動の計算式 に登場する係数\(\frac{3a^3 P}{4\mu}\)。 地震モーメントとの関係は\(M_0=3\pi a^3 P\)である。
The coefficient \(\frac{3a^3 P}{4\mu}\) in the formula of the displacement and tilt. This quantity is related to the seismic moment as \(M_0=3\pi a^3 P\).
station ユーザがパラメータの順番を間違えないための引数であり、 文字列"station"のみ可(他の文字列を指定するとエラーになる)。 この後の引数が観測点に関するパラメータであることを示す。
An argument to avoid a mistake of the order of parameters. Only the value "station" is allowed; if another string is given to here, the program finishes as an error. This argument indicates that the subsequent arguments are parameters of the station.
x 観測点の\(x\)座標。
The \(x\)-component of the station.
y 観測点の\(y\)座標。
The \(y\)-component of the station.
surface ユーザがパラメータの順番を間違えないための引数であり、 文字列"surface"のみ可(他の文字列を指定するとエラーになる)。 この後の引数が地表面に関するパラメータであることを示す。
An argument to avoid a mistake of the order of parameters. Only the value "surface" is allowed; if another string is given to here, the program finishes as an error. This argument indicates that the subsequent arguments are parameters of the ground surface.
z0 水平位置\((0,0)\)における地表面の\(z\)座標。
The \(z\)-component of the ground surface at a horizontal location \((0,0)\).
theta_x 地表面の\(x\)方向の傾き\(\theta_x\) (計算式参照)。
The slope, \(\theta_x\), of the ground surface along the \(x\)-direction; see the formula.
theta_y 地表面の\(y\)方向の傾き\(\theta_y\) (計算式参照)。
The slope, \(\theta_y\), of the ground surface along the \(y\)-direction; see the formula.
unit 引数theta_x, theta_yの単位で、 radian単位なら"rad"、 °単位なら"deg"を指定する。
The unit of the values of arguments theta_x and theta_y, which is either "rad" if the values are given in radian, or "deg" if the values are given in degrees.
component ユーザがパラメータの順番を間違えないための引数であり、 文字列"component"のみ可(他の文字列を指定するとエラーになる)。 この後の引数が成分に関するパラメータであることを示す。
An argument to avoid a mistake of the order of parameters. Only the value "component" is allowed; if another string is given to here, the program finishes as an error. This argument indicates that the subsequent arguments are parameters of the component.
component_name 計算する成分。以下のいずれかが指定可能である。
The component to calculate, which can be one of the followings.
  • "Ux" — 変位の\(x\)成分。
    The \(x\)-component of the displacement.
  • "Uy" — 変位の\(y\)成分。
    The \(y\)-component of the displacement.
  • "Uz" — 変位の\(z\)成分。
    The \(z\)-component of the displacement.
  • "Tx" — 傾斜変動の\(x\)成分。
    The \(x\)-component of the tilt.
  • "Ty" — 傾斜変動の\(y\)成分。
    The \(y\)-component of the tilt.

なお、観測点の\(z\)座標は指定した\((x,y)\)における地表面となる。
The \(z\)-component of the station is automatically set at the ground surface of \((x,y)\) specified.

◆戻り値(Return value)

指定したソースによって指定した\((x,y)\)座標の地表面に生じる地殻変動の 指定した成分の値。
The value of the specified component of the crustal deformation at the ground surface of the specified \((x,y)\) caused by the specified source.

◆使用例(Example)

double Tx=Mogi("source",0.0,0.0,-500.0,1.0,
    "station",200.0,100.0,
    "surface",0.0,10.0,12.0,"deg",
    "component","Tx");

この例では\((0,0,-500)\)の位置に \(\frac{3a^3 P}{4\mu}=1\)の等方圧力源を置いたときに 水平位置\((200,100)\)の地表面に生じる傾斜変動の\(x\)成分を計算する。 ここで地表面は水平位置\((0,0)\)で\(z=0\)、\(x\)方向の勾配が10°、 \(y\)方向の勾配が12°とする。
In this example, the \(x\)-component of tilt at the ground surface of a horizontal location \((200,100)\) caused by an isotropic source of \(\frac{3a^3 P}{4\mu}=1\) at a location \((0,0,-500)\). Here, the ground surface is set at \(z=0\) at a horizontal location \((0,0)\) and inclined by 10° and 12° in the \(x\)- and \(y\)-directions, respectively.