The documentation of “waterPML

waterPML_sub/source.h マニュアル

(The documentation of waterPML_sub/source.h)

Last Update: 2024/2/19

waterPML_sub/source.hでは地震波動ソースの設定に基づいて等価体積力の空間分布を設定する関数が定義されている。
Functions to determine the spatial distribution of an equivalent body force based on the configuration of the seismic wave source are defined in waterPML_sub/source.h.

waterPMLコマンドでは地震波動ソースを等価体積力\(f_i(\posx,t)\)で与える。等価体積力は位置の関数と時間の関数の積の重ね合わせとして \[\begin{equation} f_i(\posx,t)=\sum_{is}F_i^{(is)}(\posx)s^{is}(t) \label{eq.force} \end{equation}\] の形で表す (地震波動ソースの与え方参照)。このヘッダファイルではこのうち位置の関数である\(F_i^{(is)}(\posx)\)を与える。なお時間の関数\(s^{is}(t)\)の設定にはymaeda_opentools共通の関数である create_timefunc が用いられる。
The waterPML command gives the seismic wave source as an equivalent body force \(f_i(\posx,t)\), which is expressed as a superposition of the products of spatial and temporal functions (Eq. \ref{eq.force}; see Defining the seismic wave source). This header file determines \(F_i^{(is)}(\posx)\), which gives the spatial distribution of an equivalent body force. The time function of the force is determined by function create_timefunc, which is not specific to the waterPML command but commonly used by ymaeda_opentools.

地震波動ソースの設定はファイルで与えるが、ファイルから設定を読み込む処理はwaterPMLコマンドに特有のものではなく、 ymaeda_opentools全体の共通関数である tws_set_source が用いられる。したがって、その戻り値である struct tws_source型構造体がこのヘッダファイル内の関数における入力となる。
Although the configuration of the seismic wave source is given by a file, it is read and processed by function tws_set_source, which is not specific to this command but is shared by ymaeda_opentools. Its return value, given as a struct tws_source-type structure, is the input to the functions in this header file.

このヘッダファイル内で定義されている関数を以下に示す。各関数の詳細は関数名をクリックしてリンク先を参照のこと。
Functions defined in this header file are listed below. For details of individual functions, click the links.

関数名 Function name	機能・用途 Purpose
TF_set_f	等価体積力の空間分布を設定する。 Determine the spatial distribution of an equivalent body force.
TF_calculate_source_intensity	各格子セルに与える地震波動ソースの強度を設定する。 Determine the intensity of a seismic wave source assigned to each grid cell.
TF_set_MandFandT	大きさと向きで定義された地震波動ソースを直交座標に沿った行列成分表記に変換し、その空間分布を設定する。 Convert the mechanism of a seismic wave source, given by the magnitude and orientation, to matrix representations based on a cartesian coordinate and determine their spatial distributions.
TF_set_f_index	等価体積力の代入先の配列要素番号(方法3-2)が半格子点の通し番号(方法2)から得られるように struct grid型構造体stgのメンバindexの値を設定する。また、等価体積力を与える半格子点の個数をカウントする。 Define the index of an array to represent the equivalent body force (method 3-2) for each consecutive index of a half-grid node (method 2) by setting the values of member index of a struct grid-type structure stg, and count the number of half-grid nodes where the equivalent body force is to be given.
TF_set_f_value	1つの地震波動ソース要素(is)について等価体積力の空間分布を設定する。 Determine the spatial distribution of an equivalent body force for a seismic wave source element (is).