Compute steady state denudation using Eulerian description for the evolution of concentration

Eulerian formulation

solv_ss_eros_eul(Cobs, p, S, L, Cobs_e = 0, method = "single", n = 10000)

Arguments

Cobs

Uncertainty on measured concentration (at/g), optional but required ffor method MC and MCMC

p

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)

S

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

L

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

method

One of single (default), MC, MCMC