input_Paris700_pass_2

# --------
# Input file 
# --------
# Data_file Paris700.txt
Data_file archeomagnetic_data/Paris700_within_credible_interval.txt
#
# Column format for datafile, in the order age, delta age, intensity, delta intensity, stratification
# in Python ordering - where column 0 is the first, column 1 is the second etc.
# Note: stratification = -1 means no stratification column is loaded (and so the data is assumed unstratified)
File_format 0 1 2 3 -1
# RJMCMC parameters
Burn_in 50000
Nsamples  1050000
# Discretistion size of the model
model_discretisation 200
# Chain-forming parameters:  show, thinning
Chain_parameters 10000 100 
# Age distribution: U = uniform, N = normal
Age_distribution U
# Running mode: 1 - normal, 0 - set all likelihoods to 1 to recover the prior distributions.
running_mode 1 
# Age Min and Max for model parametrisation
Age_bounds -1160 2000
# Parameters that describe the various model perturbations:
# sigma_move, sigma_change, sigma_birth
Sigmas 200 20 8
# Fraction of ages to change in a proposal (i.e. 20 means that 1/20 of them will change)
Age_frac 20.0
#
# --------
# Prior information
# --------
#  Intensities: I_min, I_max in micro Tesla
Intensity_prior 30 100 
# Number of change points: K_min, K_max
Num_change_points 1 50 
#
# --------
# Outputs
# --------
# Directory for all outputs:
Outputs_directory Outputs_Paris700_pass_2
# Credible interval or any non-positive number for none.
Credible 95
# Number of bins for posterior marginals:
Nbins 200
# Name and write frequency of model file. The line below causes the code to write every 100th model to "models.dat". Enter a frequency of -1 for no output. 
output_model models.dat 10
# Joint distribution output: if frequency is -1 then nothing is output; otherwise this defines the frequency
output_joint_distribution_freq 1