Compute concentration evolution along a time-depth history

Lagrangian formulation

solv_conc_lag(t, z, C0, Psp0, Pmu0, lambda, S, final = FALSE)

Arguments

t

time vector (a)

z

depth vector (g/cm2), same length as t

C0

initial concentration (at/g)

Psp0

SLHL spallation production profile (at/g/a) at depth corresponding to z

Pmu0

muonic production profil at depth corresponding to z

lambda

radioactive decay constant (1/a)

S

scaling parameters for PsP0 and Pmu0 (2 columns), at times corresponding to t (same number of rows)

final

if TRUE only compute the final concentration (default=FALSE)