| stg | 異なる方法で表した位置を対応付けるための構造体。
この構造体では
半格子点の通し番号ig(方法2)
を元に同じ地点における
特定の量のための通し番号(配列要素番号用; 方法3-2)
を求めるための対応関係を与える。 A structure to relate the locations expressed by different methods. This structure defines the consecutive index for a specific quantity (for array indices; method 3-2) for each location given by the consecutive index for a half-grid node (method 2). |
| parameters | コマンドライン引数で与えたパラメータ一式を格納した構造体。 A structure composed of the parameters given by command-line arguments. |
| Nigv | 下記のigv2ig[i]の配列要素数を並べた1次元配列(引数:i)。
igv2igは2次元配列であるが、その第2引数の個数は第1引数iによって異なるので、
個数をiを引数とする1次元配列で与える必要がある。 A 1-D array (argument: i) composed of the numbers of array components of igv2ig[i] defined below. This must be a 1-D array with argument i because the number of 2nd arguments of the 2-D array igv2ig depends on the 1st argument i. |
| igv2ig | 異なる方法で表した位置を対応付けるための2次元配列。
第1引数は速度\(V_i^k\)の成分番号\(i\)、
第2引数は位置を表す\(V_i^k\)用の
通し番号igv(ループ制御用; 方法3-1)
であり、igv2ig[i][igv]はその地点における
半格子点の通し番号ig(方法2)
を表す。 A 2-D array that relates the locations expressed by different methods. The 1st argument is the index i of a velocity component \(V_i^k\), and the 2nd argument is the consecutive index igv (for loop; method 3-1) that represents a location for \(V_i^k\). Each array component igv2ig[i][igv] indicates the consecutive index ig of the half-grid node (method 2) at that location. |
| it | 速度場を計算したい時刻\(t\)に対応する時刻番号。
計算の時間刻み(パラメータdtの値)を\(\Delta t\)として
\(t=\mbox{(it)}\Delta t\)である。 An index of time defined as \(t=\mbox{(it)}\Delta t\), where \(t\) is the target time of the computation of the velocity field and \(\Delta t\) is the time stepping of the computation given by parameter dt. |
| type | 半格子点の属性を並べた2次元配列。詳細は
こちらのページ
参照。 A 2-D array composed of the attributes of half-grid nodes. For detail, see this page. |
| c0 | 差分方程式の係数\(c_0^k(\posx)\)
(計算式
の(7)式)
を格納した2次元配列。
第1引数が\(k\)の値、
第2引数が位置\(\posx\)に対応する\(c_0^k\)用の
通し番号(位置の表現方法3-2)。 A 2-D array composed of the coefficients \(c_0^k(\posx)\) of difference equations (Eq. 7 of the formula). The 1st argument is \(k\), and the 2nd argument is the consecutive index (method 3-2 to represent the location) for \(c_0^k\) at the location \(\posx\). |
| c1 | 差分方程式の係数\(c_1^k(\posx)\)
(計算式
の(8)式)
を格納した2次元配列。
第1引数が\(k\)の値、
第2引数が位置\(\posx\)に対応する\(c_1^k\)用の
通し番号(位置の表現方法3-2)。 A 2-D array composed of the coefficients \(c_1^k(\posx)\) of difference equations (Eq. 8 of the formula). The 1st argument is \(k\), and the 2nd argument is the consecutive index (method 3-2 to represent the location) for \(c_1^k\) at the location \(\posx\). |
| T | 時刻\(t-\Delta t/2\)における応力成分\(\tau_{ij}^l\)を並べた3次元配列。
第1引数は\((i,j)\)の組を表す
インデックスiT、
第2引数は\(l\)であり、
第3引数は位置を表す\(\tau_{ij}^l\)用の
通し番号(位置の表現方法3-2)
である。 A 3-D array composed of the stress components \(\tau_{ij}^l\) at time \(t-\Delta t/2\). The 1st argument is the index iT corresponding to the components \((i,j)\), the 2nd argument is \(l\), and the 3rd is the consecutive index (method 3-2) that represents a location for \(\tau_{ij}^l\). |
| f | 地震波動ソースの等価体積力(位置の関数)を表す3次元配列。
地震波動ソースの与え方
で説明したように、waterPMLコマンドでは地震波動ソースを
位置の関数\(F_i^{(is)}(\posx)\)と時間の関数\(s^{(is)}(t)\)の積の
重ね合わせとして表現するが、
このうちの\(F_i^{(is)}(\posx)\)をこの引数で与える。
3次元配列の第1引数はソース番号is、第2引数は力の成分番号i、
第3引数は位置を表す地震波動ソース用の
通し番号(位置の表現方法3-2)
である。 A 3-D array to represent the equivalent body force (a function of position) of each seismic wave source. As is explained in Defining the seismic wave source, the waterPML command expresses the seismic wave source by summing the product of a function of position \(F_i^{(is)}(\posx)\) and a function of time \(s^{(is)}(t)\). This argument gives \(F_i^{(is)}(\posx)\). The 1st argument of the 3-D array is the index of source is, the 2nd is the direction i and the 3rd is the consecutive index (method 3-2) that represents a location for the seismic wave source. |
| stfun | 地震波動ソースの等価体積力(時間の関数)を並べた1次元配列。
stfun[is]は上記\(s^{(is)}(t)\)を表す。 A 1-D array composed of the equivalent body forces (functions of time) of the seismic wave source; stfun[is] represents \(s^{(is)}(t)\) introduced above. |
| V | 速度成分\(V_i^k\)を並べた3次元配列。
第1引数は\(i\)、第2引数は\(k\)であり、
第3引数は位置を表す\(V_i^k\)用の
通し番号(位置の表現方法3-2)
である。
関数呼び出し時には時刻\(t-\Delta t\)での値を与える。
関数内で時刻\(t\)での値に書き換えられる。 A 3-D array composed of the velocity components \(V_i^k\). The 1st argument is \(i\), the 2nd is \(k\), and the 3rd is the consecutive index (method 3-2) that represents a location for \(V_i^k\). Give the values at time \(t-\Delta t\) as the inputs to the function. The values are replaced by the values at time \(t\) within the function. |
| ファイル名 File name |
出力されるデータ Data written |
output_dir/debug/Tijl.ttime
|
時刻\(t-\Delta t/2\)での応力場。 The stress field at time \(t-\Delta t/2\). |
output_dir/debug/Vik.ttime
|
時刻\(t\)での速度場。 The velocity field at time \(t\). |
| (\ref{eq.Vik})式 Eq. (\ref{eq.Vik}) |
ソースコード Source code |
| \(V_i^k(\posx_{i-1/2},t_{-1})\) | 更新前(Before update):
|
| \(V_i^k(\posx_{i-1/2},t)\) | 更新後(After update):
|
| \(\tau_{ik}^0\left(\posx_{i-1/2}^{k+1/2},t_{-1/2}\right)\) |
|
| \(\tau_{ik}^0\left(\posx_{i-1/2}^{k-1/2},t_{-1/2}\right)\) |
|
| \(\tau_{ik}^1\left(\posx_{i-1/2}^{k+1/2},t_{-1/2}\right)\) |
|
| \(\tau_{ik}^1\left(\posx_{i-1/2}^{k-1/2},t_{-1/2}\right)\) |
|
| \(\tau_{ik}^2\left(\posx_{i-1/2}^{k+1/2},t_{-1/2}\right)\) |
|
| \(\tau_{ik}^2\left(\posx_{i-1/2}^{k-1/2},t_{-1/2}\right)\) |
|
| \(c_0^k(\posx_{i-1/2})\) |
|
| \(c_1^k(\posx_{i-1/2})\) |
|
| \(f_i(\posx_{i-1/2},t_{-1/2})\) |
|
| \(\posx_{i-1/2}\) |
|
| \(\posx_{i-1/2}^{k+1/2}\) |
|
| \(\posx_{i-1/2}^{k-1/2}\) |
|
| |||||||||||||||||||||