Monte Carlo exploration for a depth profile

Eulerian approach

depth_profile_mc(C, C_e, depth, prm, L, S, age, ero, inh, rho, n1 = 0, n2 = 0)

Arguments

C

concentrations vector (at/g)

C_e

uncertainties on the concentrations vector (at/g)

prm

production and decay parameters (4 elements vector) p[1] -> unscaled spallation production rate (at/g/a) p[2] -> unscaled stopped muons production rate (at/g/a) p[3] -> unscaled fast muons production rate (at/g/a) p[4] -> decay constant (1/a)

L

attenuation lengths (3 elements vector) L[1] -> neutrons L[2] -> stopped muons L[3] -> fast muons

S

scaling factors (2 elements vector) S[1] -> scaling factor for spallation S[2] -> scaling factor for muons

age

age range to explore (a). 2 elements vector. If 1 element this parameter is fixed

ero

erosion rate range to explore (m/Ma). 2 elements vector. If 1 element this parameter is fixed

inh

inherited concentration range to explore (at/g). 2 elements vector. If 1 element this parameter is fixed

rho

density range to explore (g/cm3). 2 elements vector. If 1 element this parameter is fixed

n1

number of runs with random

n2

z

depth vector (g/cm2)