cbk.lame.delivery: Estimate delivery rate of chem

Description Usage Arguments Details Value Examples

View source: R/cbk.lame.delivery.R

Description

Estimate delivery rate of chem. Calculation is made in dimension of CGS.

Usage

1
2
cbk.lame.delivery(pmlame, intlame, sputter.rate = 1.1e-09,
  verbose = FALSE)

Arguments

pmlame

A pmlame with row of stone and column of chem [g/g].

intlame

A pseudo-pmlame of ion intensity with row of acq and column of isomeas [cps]. Typical rowname is ‘trc_cpx_klb1@1’. This is like array of I(Li7), ..., I(Sr88).

sputter.rate

Sputter rate [g/s] (default: 1.1e-9 that corresponds to abbration of 40 s to produce a pit with diameter of 25 micron, depth of 30 micron, and density of 3 g/cc). This also can be vector with name of acqlist, that should be consistent with rownames of intlame.

verbose

Output debug info (default: FALSE).

Details

With ion Si29 intensity 3.85e+05 [cps], element SiO2 abundance 59 wt estimated to be 1.25 ppm.

Value

Delivery rate of element from solid target to detector.

Examples

1
2
3
4
5
6
7
8
9
pmlame0      <- data.frame(row.names=c("ref-gl-tahiti","ref-cpx-klb1"),SiO2=c(0.590,0.520),Li=c(20.8e-06,1.37e-06),Sr=c(107e-06,36.1e-06))
pmlame1      <- cbk.lame.reduce(pmlame0)
intlame0     <- data.frame(row.names=c("ref_gl_tahiti@1","ref_cpx_klb1@1"),Si29=c(3.85e+05,3.95e+05),Li7=c(1.87e+03,1.73e+02),Sr88=c(3.24e+04,1.50e+04))
radius       <- 12.5e-6 # Radius of a pit [m]
depth        <- 30e-6   # Depth of a pit [m]
t            <- 40      # Period of irradiation [s]
density      <- 3       # Density of solid target [g/cc]
sputter.rate <- 3.14 * (radius * 100)^2 * (depth * 100) * density / t # [g/s]
cbk.lame.delivery(pmlame1,intlame0,sputter.rate)

misasa/chelyabinsk documentation built on Nov. 24, 2020, 5:47 a.m.