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'. 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 computation, 
#                            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' 
# 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