Input Files

Input files

Parameter file (params.in)

A parameter file, in which other input files and tuning parameters are specified, must exist in the current directory with the name 'params.in'. Otherwise, users need to give a path to the parameter file (with any file name) as the first arguments in the command line: mpirun -np 20 rf_inv somewhere/params.in, for example.

In the parameter file, all parameters must appear in the exact order as shown in the example below. Note that lines beginning with '#' are ignored.

Format

Line #	Notation	Description
1	Dir	Directory name where outputs are placed
2	N_burn	Number of MCMC steps in burn-in period
3	N_iter	Number of total MCMC steps
4	N_skip	MCMC sampling interval (sample/step)
5	N_chains	Number of MCMC chains per process
6	N_cool	Number of non-tempered MCMC chains per process
7	T_high	Maximum temperature used in parallel tempering
8	I_seed	Seed number to generate random numbers
9	N_trc	Number of input waveforms
9+(1, ..., N_trc)	Rayp	Ray parameters (s/km)
9+N_trc+(1, ..., N_trc)	A_gauss	Parameters for Gaussian low-pass filter
9+2*N_trc+(1, ..., N_trc)	I_pha	Incident phase type (P: 1, S: -1)
10+3*N_trc	N_fft	Number of elements in time series used for forward computing
10+3*N_trc+(1, ..., N_trc)	Obs_file	SAC format file names for input waveforms
11+4*N_trc	T_start, T_end	Start and end times of time window to be analyzed
12+4*N_trc	deconv_mode	Deconvoluition mode (0: w/o deconvolution, 1: w/ deconvolution)
13+4*N_trc	Sta_dep	Station depth below the sea surface (km)
14+4*N_trc	Vel_file	File names for reference velocities
15+4*N_trc	Vp_mode	Treatment of P wave velocity (0: Fixed at reference model, 1: Solved as unknown)
16+4*N_trc	k_min, k_max	Minimum and maximum numbers of layter interfaces
17+4*N_trc	z_min, z_max	Minimum and maximum of interface depths
18+4*N_trc	h_min	Minimum thickness of layer (km)
19+4*N_trc	dVs_min, dVs_max	Minimum and maximum of S-wave velocity perturbation (km/s)
20+4*N_trc	dVp_min, dVp_max	Minimum and maximum of P-wave velocity perturbation (km/s) (Not used when Vp_mode=0)
21+4*N_trc	prior_mode	Form of prior distribution for velocity perturbations (1: Laplacian, 2: Gaussian)
22+4*N_trc	sig_dVs_prior	Standard deviation of prior probability of S wave velocity perturbation (km/s)
23+4*N_trc	sig_dVp_prior	Standard deviation of prior probability of P wave velocity perturbation (km/s) (Not used when Vp_mode=0)
23+4*N_trc+(1, ..., N_trc)	sig_min, sig_max	Minimum and maximum of noise standard deviation
24+5*N_trc	del_z	Standard deviation for changing interface depth (km)
25+5*N_trc	del_dVs	Standard deviation for changing S wave velocity perturbation (km/s)
26+5*N_trc	del_dVp	Standard deviation for changing P wave velocity perturbation (km/s)
27+5*N_trc	del_sig	Standard deviation for changing standard devation of noise
28+5*N_trc	N_bin_z	Number of bins in terms of depths
29+5*N_trc	N_bin_dVs	Number of bins in terms of S wave velocity perturbation
30+5*N_trc	N_bin_dVp	Number of bins in terms of P wave velocity perturbation
31+5*N_trc	N_bin_vpvs	Number of bins in terms of Vp/Vs ratio
32+5*N_trc	N_bin_sig	Number of bins in terms of noise standard deviation
33+5*N_trc	amp_min, amp_max	Minimum maximum amplitudes to be saved
34+5*N_trc	Vp_min, Vp_max	Minimum maximum P wave velocity to be saved (km/s)
35+5*N_trc	Vs_min, Vs_max	Minimum maximum S wave velocity to be saved (km/s)
36+5*N_trc	VpVs_min, VpVs_min	Minimum maximum Vp/Vs ratio to be saved

Example

# Example of a parameter file for RF_INV (params.in)
#
# Output directory (must exist)
'Out'
#
#---------------------------------------------------------
# Iteration
#---------------------------------------------------------
# Iteration number in burn-in period
10000
# Iteration number after burn-in
50000
# Iteration number per saving one sample
10
# Number of chains per processor
5
# Number of non-tempered chains per processor (>= 1)
1
# Highest temperature (>= 1.0)
50.0
# Random number seed
56763298
#
#---------------------------------------------------------
# Input Data
#---------------------------------------------------------
# Number of input waveforms
3
# Ray parameter (s/km) (* Need to repeat as many time as the number of input waveforms) 
#     Note: If ray parameters are common for all input waveforms, then forward computation occurs only once.
#           Otherwise, the computation is repeated as many times as the number of input. 
0.06
0.06
0.06
# Gaussian parameter (i.e., a in G(omega) = exp(-(omega^2/(4a^2))) 
#                    (* Need to repeat as many time as the number of input waveforms)
8.0
4.0
1.0
# Incident phase type (P: 1, S: -1)
#                     (* Need to repeat as many time as the number of input waveforms)
1
1
1
# Sample number used in FFT (higher number takes longer time for forward computing, 
#                            but may offer better resolution in the frequency domain)
2048
# Observed RF file (* Need to repeat as many time as the number of input waveforms)
#    Note: Observed RF file must be SAC format. 
#          A header 'delta' is required to be common for all input waveforms. 
'KMB06_120_150_a8.0_MC3.sac' 
'KMB06_120_150_a4.0_MC3.sac' 
'KMB06_120_150_a1.0_MC3.sac' 
# Start/End time
#    Note: Direct P is assumed to arrive at t=0  
-1.0 10.0
# Forward computation mode (0: w/o deconvolution, 1: w/ deconvolution)
1
# Station depth (km below the sea surface)
#    Note: A station is assumed to be located on the seafloor. 
2.499
# Reference velocity file
#    Format: z(km)  Vp(km/s)  Vs(km/s)
#    Note: Need evenly spaced grids in depth
'KMB06.vel.formatted'
#
#---------------------------------------------------------
# Inversion mode
#---------------------------------------------------------
# Vp mode (0: Fixed at reference model, 1: Solved)
0
#
#---------------------------------------------------------
# Prior probability 
#---------------------------------------------------------
# Min./Max. # of interfaces [min., max.)
1 31
# Min./Max. of interface depth (km; below sea surface)
2.499 20.0
# Minimum thickness of layer (km)
0.05
# Min./Max. of dVs (km/s)
-1.5 1.5
# Min./Max. of dVp (km/s)
-2. 2.0
# Form of prior probability distribution for velocity perturbations (1: Laplacian, 2: Gaussian) 
# Standard deviation for dVs prior
0.2
# Standard deviation for dVp Prior
1.0 
# Min./Max. of noise sigma
0.005 0.08
#
#---------------------------------------------------------
# Proposal
#---------------------------------------------------------
# standard deviation for depth proposal
0.03
# standard deviation for dVs proposal
0.02
# standard deviation for dVp proposal
0.02
# standard deviation for noise sigma proposal 
0.005
#
#---------------------------------------------------------
# Figure
#---------------------------------------------------------
# Number of bins for depth
100
# Number of bins for Vs
50
# Number of bins for Vp
50
# Number of bins for noise sigma
50
# Number of bins for amplitudes
100
# Min./Max. amplitudes to be displayed
-0.6 0.6
# Min./Max. Vp to be displayed
0.1 8.6
# Min./Max. Vs to be displayed
0.001 5.0

Velocity file

Free-format ASCII file for a reference velocity model.

Format

1st Column	2nd Column	3rd Column
Depth of layer top	P wave velocity (km/s)	S wave velocity

Note

The depth increment must be constant.