パラメタファイルの記法
本コードの挙動は原則として一つの入力パラメタファイルで制御される．
本章では主として3次元コードの入力を説明するが，2次元P-SV,
SHコードについても， 一部変数が不要なだけで同じパラメータで動作する．
1行に1変数を
の形式で記述する．変数値の表現にはFortranの記法を用いる．たとえば論理型は
.true. や .false.で記述される．
原則として自由書式であり，上記変数行以外には何を書いてあってもよいが，特に!と#で始まる行は
明示的にコメント行として読みとばされる．コメントは変数と同じ行にあっても良い．たとえば
| nx = 1024 ! number of grids
|
のような記述が可能である．
変数の順番は問わず，どの変数がどの行に書かれていても良い．
ただし，指定されていない変数がある場合には，プログラム動作の支障の具合に応じて，デフォルト値が使われる場合と動作を停止する場合がある．
もし同じ変数の記述が複数回ある場合には，最初に読み込まれたものが採用される．
変数指定の等号の両側には空白があって良い．ただし，整数もしくは実数のパラメタが負の値の場合，マイナス記号と数字の間に空白を置くことは許されない．また，ディレクトリの区切り記号を含む文字列変数は，シングルもしくはダブルクォーテーションで括って与える必要がある．
Caution
動作に必須ではないパラメタがパラメタファイルに記載されていない場合，利用されたデフォルト値が以下のように標準エラー出力に表示される．
| [info] key XXXX is not found.
[info] Use default value YYYY instead.
|
この記載を見逃すと予期せぬ動作につながることがある．
パラメタファイルの例
以下にサンプルパラメタを示す．以降の各節で，機能ごとに挙動とパラメタの解説を行う．
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215 | !! ----------------------------------------------------------------------- !!
!!
!! SWPC input file
!!
!! ----------------------------------------------------------------------- !!
!! ----------------------------------------------------------------------- !!
!! Control
!!
title = 'swpc' !! exe title: used for output filenames
odir = './out' !! output directory
ntdec_r = 50 !! screen report timing (1/cycle)
strict_mode = .false. !! all parameters to be explicitly definied
!! ----------------------------------------------------------------------- !!
!! Model/Grid Size and Area
!!
nproc_x = 2 !! parallelization in x-dir
nproc_y = 2 !! parallelization in x-dir
nx = 384 !! total grid number in x-dir
ny = 384 !! total grid number in y-dir
nz = 384 !! total grid number in z-dir
nt = 1000 !! time step number
dx = 0.5 !! grid width in x-dir
dy = 0.5 !! grid width in y-dir
dz = 0.5 !! grid width in z-dir
dt = 0.02 !! time step width
vcut = 1.5 !! minimum velocity
!- smaller velocities will be increased
xbeg = -96.0 !! minimum in x-dir
ybeg = -96.0 !! minimum in y-dir
zbeg = -10.0 !! minimum in z-dir
tbeg = 0.0 !! start time
clon = 139.7604 !! center longitude
clat = 35.7182 !! center latitude
phi = 0.0 !! horizontal coordinate rotation
!- measured clockwise from the north
fq_min = 0.02 !! minimum freq. for Q-constant model
fq_max = 2.00 !! maximum freq. for Q-constant model
fq_ref = 1.0 !! ref. freq. for physical dispersion
!! ----------------------------------------------------------------------- !!
!! Snapshot Output
!!
snp_format = 'netcdf' !! snapshot format (netcdf)
xy_ps%sw = .false. !! P&S amp. for xy section
xz_ps%sw = .true. !! P&S amp. for xz section
yz_ps%sw = .false. !! P&S amp. for yz section
fs_ps%sw = .false. !! P&S amp. for free surface
ob_ps%sw = .true. !! P&S amp. for ocean bottom
xy_v%sw = .false. !! 3-comp. velocity for xy section
xz_v%sw = .true. !! 3-comp. velocity for xz section
yz_v%sw = .false. !! 3-comp. velocity for yz section
fs_v%sw = .false. !! 3-comp. velocity for free surface
ob_v%sw = .true. !! 3-comp. velocity for ocean bottom
xy_u%sw = .false. !! 3-comp. disp. for xy section
xz_u%sw = .true. !! 3-comp. disp. for xz section
yz_u%sw = .false. !! 3-comp. disp. for yz section
fs_u%sw = .false. !! 3-comp. disp. for free surface
ob_u%sw = .true. !! 3-comp. disp. for ocean bottom
z0_xy = 7.0 !! depth for xy cross section
x0_yz = 0.0 !! x-value for yz cross section
y0_xz = 0.0 !! y-value for xz cross section
ntdec_s = 5 !! time decimation of snapshot
!- (specify 1 for no decimation)
idec = 2 !! x-decimation for snapshot
jdec = 2 !! y-decimation for snapshot
kdec = 2 !! z-decimation for snapshot
!! ----------------------------------------------------------------------- !!
!! Waveform Output
!!
sw_wav_v = .true. !! velocity trace output at stations
sw_wav_u = .false. !! displacement trace output at stations
sw_wav_stress = .false. !! stress tensor trace
sw_wav_strain = .false. !! strain tansor trace
ntdec_w = 5 !! time decimation of waveform output
st_format = 'xy' !! station format: 'xy' or 'll'
fn_stloc = './example/stloc.xy' !! station location file
wav_format = 'sac' !! 'sac' or 'csf' ('sac' recommended)
wav_calc_dist = .false. !! Calculate epicentral distance
!! ----------------------------------------------------------------------- !!
!! Earthquake Source
!!
!! Moment tensor source format:
!! xymoij / xym0dc / llm0ij / llm0dc / xymwij / xymwdc / llmwij / llmwdc
!! Body force source format:
!! xy or ll
stf_format = 'xym0ij'
!! Basis source time function
!! 'boxcar' / 'triangle' / 'herrmann' / 'kupper' / 'cosine' / 'texp'
stftype = 'kupper'
fn_stf = "./example/source.dat" !! Source grid file name
!! source depth correction
!! 'asis':use z value, 'bd{i}': i-th boundary (i=0...9)
sdep_fit = 'asis'
!! --------------------------------------------------------------------- !!
!! Body force source mode
!!
bf_mode = .false.
!! --------------------------------------------------------------------- !!
!! Plane wave source mode
!!
pw_mode = .false. !! plane wave input. neglects fn_stf
pw_ztop = 100. !! top z-coordinate of the initial plane wave
pw_zlen = 30. !! wavelength of the initial plane wave
pw_ps = 'p' !! 'p' P-wave 's' S-wave
pw_strike = 0.0 !! strike direction of plane wave (deg.)
pw_dip = 0.0 !! dip of plane wave (deg.)
pw_rake = 0.0 !! rake of plane S-wave polarization (deg.)
!! ----------------------------------------------------------------------- !!
!! Absorbing Boundary Condition
!!
abc_type = 'pml' !! 'pml' or 'cerjan'
na = 20 !! absorbing layer thickness
stabilize_pml = .false. !! avoid low-v layer in PML region
!! ----------------------------------------------------------------------- !!
!! Velocity model
!!
vmodel_type = 'lhm' !! velocity model type 'uni'/'grd'/'lhm'
is_ocean = .true. !! topography z<0 is covered by ocean
topo_flatten = .false. !! Force topography variation to zero (formerly is_flatten)
munk_profile = .true. !! velocity gradient inside the seawater column
earth_flattening = .false. !! Earth-flattening tranformation
!! --------------------------------------------------------------------- !!
!! For uniform velocity model 'uni'
!!
vp0 = 5.0 !! P-wave velocity [km/s]
vs0 = 3.0 !! S-wave velocity [km/s]
rho0 = 2.7 !! mass density [g/cm^3]
qp0 = 200 !! Qp
qs0 = 200 !! Qs
topo0 = 0 !! topography location
!! --------------------------------------------------------------------- !!
!! For GMT grid file input 'grd' ( requires netcdf library )
!!
dir_grd = '${DATASET}/vmodel/ejivsm/' !! directory for grd file
fn_grdlst = './example/grd.lst' !! grd file list
node_grd = 0 !! input MPI node
!! --------------------------------------------------------------------- !!
!! For layered homogeneous medium model ('lhm')
!!
fn_lhm = 'example/lhm.dat' !! 1D velocity structure
!! --------------------------------------------------------------------- !!
!! For random medium models
!!
dir_rmed = './in/' !! location of random medium file
fn_grdlst_rmed = './example/grd.lst' !! grd file list
rhomin = 1.0 !! minimum density threshold
fn_rmed0 = 'dummy.ns' !! vel. purturb. on a uniform media
!! ----------------------------------------------------------------------- !!
!! Checkpoint/Restart
!!
is_ckp = .false. !! perform checkpoint/restart
ckpdir = './out/ckp' !! output directory
ckp_interval = 1000000 !! interval for checkpoint check（1/cycle）
ckp_time = 1000000. !! checkpoint time
ckp_seq = .true. !! sequential output mode
!! ----------------------------------------------------------------------- !!
!! Reciprocity Green's Function Mode
!!
green_mode = .false. !! reciprocity Green's function mode
green_stnm = 'st01' !! virtual station name from fn_stlst
green_cmp = 'z' !! virtual source direction 'x', 'y', 'z'
green_trise = 1.0 !! rise time
green_bforce = .false. !! also calc. body force Green's function
green_maxdist = 550. !! horizontal limit of source grid
green_fmt = 'llz' !! list file format: 'xyz' or 'llz'
fn_glst = 'example/green.lst' !! Green's function grid point list
!! ----------------------------------------------------------------------- !!
!! MISC
!!
stopwatch_mode = .false. !! measure computation time at routines
benchmark_mode = .false. !! benchmark mode
ipad = 0 !! memory padding size for tuning
jpad = 0 !! memory padding size for tuning
kpad = 0 !! memory padding size for tuning
|
以下の各節では，機能ごとにパラメタの詳細について述べる．