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.

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
# Min./Max. of dVs (km/s)
-1.5 1.5
# Min./Max. of dVp (km/s)
-2. 2.0
# 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