depth_averaged_prod: Average nuclide production rate in soil

In VincentGodard/TCNtools: Terrestrial Cosmogenic Nuclides Calculations

Average nuclide production rate in soil

Description

Compute average nuclide production rates for vertically mixed soil column

Usage

depth_averaged_prod(h, p, L, S)

Arguments

h	soil thickness (in g/cm2)
p	production and decay parameters for the first nuclide (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 length (3 elements vector in g/cm2) 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

Examples

#' #' data("prm") # production and decay data
p = prm
data("Lambda") # attenuation length data
L = Lambda
altitude = 1000 # site elevation in m
latitude = 20 # site latitude in degrees
P = atm_pressure(alt=altitude,model="stone2000") # atmospheric pressure at site
S = scaling_st(P,latitude) # Stone 2000 scaling parameters
rhob = 2.65 # bedrock density (g/cm3)
rhos = rhob/2 # soil density (g/cm3)
h = seq(0,300,length.out = 100)
plot(NA,xlim=range(h)/100,ylim=c(2,6),
     xlab="Mixing depth (m)",ylab="Production rate (at/g/a)")
prod_soil = depth_averaged_prod(h*rhos,p[,1],L,S)
lines(h/100,prod_soil,lwd=3)

VincentGodard/TCNtools documentation built on Jan. 29, 2023, 9:23 a.m.