Main parameter file

The main parameter file contains user's choice on tuning parameters. This filename must be passed to SEIS_FILO's programs as the first command-line argument. See parameter list for descriptions of each parameter.

EXAMPLE

recv_func_fwd [path to parameter file]
disper_fwd [path to parameter file]
mpirun -np 20 joint_inv [path to parameter file]

Format

• Specify one parameter per line
• Parameter's name should be on the left side of "=" and the value on the right side
• Comment out by "#"
• Order insensitive
• Empty lines are just ignored

# Snippet of the parameter file
n_iter = 100000
n_chain = 20
# The following lines give the same results.
n_burn = 50000
n_ bu rn= 500 0 0

Required parameters for disper_fwd

• freq_or_period: dispersion measurement type (freq: dispersion curve as a function of frequency; period: as a function of period)
• xmin: minimum frequecny or period of dispersion curve (Hz or s)
• xmax: maximum frequency or period of dispersion curve (Hz or s)
• dx: freqeuncy or period interval of dispersion curve (Hz or s)
• cmin: minimum phase velocity for root search (km/s)
• cmax: maximum phase velocity for root search (km/s)
• dc: inverval of phase velocity (km/s)
• disper_phase: phase type (R: Rayleigh, L: Love)
• n_mode: mode number (n_mode >= 0)
• vmod_in: filename for input velocity model
• disper_out: output filename

Required parameters for recv_func_fwd

• n_smp: number of samples in receiver function data
• delta: time interval of receiver function (s)
• rayp: ray parameter (s/km)
• a_gauss: Gaussian low-pass filter parameter
• t_pre: time before direct P arrival in receiver function (s)
• rf_phase: phase type (P: P-wave, S: S-wave)
• deconv_flag: whether deconvolved by Z component or not (T: with deconvolution, F: without deconvolution)
• correct_amp: amplitude correction to account for energy loss by filtering (T: with correction, F: without correction)
• damp: water-level damping factor
• vmod_in: filename for input velocity model
• recv_func_out: output filename

Required parameters for joint_inv

• recv_func_in : filename of receiver function data summary file
• disper_in: filename of dispersion curve data summary file
• solve_anomaly: solve for velocity anomaly (T) or aboslute velocity (F)
• ref_vmod_in: Filename of reference velocity model
• solve_vp: whether solve for Vp (T: solve Vp, F: not)
• solve_rf_sig: whether solve receiver function data error (T) or not (F)
• solve_disper_sig: whether solve dispersion curve data error (T) or not (F)
• is_sphere: Earth flattening is applied (T) or not (F).
• is_ocean: whether model includes ocean (T: with ocean layer, F: without ocean layer)
• ocean_thick: ocean layer thickness (only required when is_ocean=T)
• vp_bottom: Vp of bottom layer (km/s).
• vs_bottom: Vs of bottom layer (km/s).
• rho_bottom: Density of bottom layer (g/cm^3).
• n_iter: number of iterations
• n_burn: number of iterations for burn-in period
• n_corr: number of iterations each time before saving a model
• n_chain: number of MCMC chains per process
• n_cool: number of non-tempered MCMC chains per process
• temp_high: maximum temperature
• k_min: minimum number of layer interfaces

• k_max: maximum number of layer interfaces

• vp_min: minimum Vp (km/s)

• vp_max: maximum Vp (km/s)
• dvp_sig: Prior width for Vp anomaly (km/s)
• vs_min: minimum Vs (km/s).
• vs_max: maximum Vs (km/s).
• dvs_sig: Prior width for Vs anomaly (km/s)
• z_min: minimum interface depth (km)
• z_max: maximum interface depth (km)
• dev_vs: standard deviation of proposed Vs perturbation (km/s)
• dev_vp: standard deviation of proposed Vp perturbation (only required when solve_vp=T) (km/s).
• dev_z: standard deviation of proposed layer depth perturbation (km).
• n_bin_vs: number of Vs bin for display.
• n_bin_vp: number of Vp bin for display
• n_bin_z: number of depth bin for display.