/
CRN_global_variables.hpp
117 lines (87 loc) · 4.52 KB
/
CRN_global_variables.hpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
/*-----------------------------------------------------------------------------------------
Global parameters for Rocky Coast Cosmogenic Radionuclide Model
Martin Hurst
British Geological Survey
December 2014
-----------------------------------------------------------------------------------------*/
/**
@file CRN_global_variables.hpp
@author Martin D. Hurst (mhurst@bgs.ac.uk)
@date 14/09/2015
@brief Global variables for modelling cosmogenic nuclide accumulation.
*/
#include <cmath>
#include <vector>
using namespace std;
#ifndef CRN_GLOBALS_HPP
#define CRN_GLOBALS_HPP
//Sea level high latitude production rates
///@brief 10Be spallation reference production rate.
///@details Spallation (a/g/yr) calibrated 10Be production rate of Lifton et al. (2014), reported in Borchers et al. (2016.
float Po_10Be_Spal = 4.009;
///@brief 14C spallation reference production rate.
///@details Spallation (a/g/yr) calibrated 14C production rate of Lifton et al. (2014), reported in Borchers et al. (2016)..
float Po_14C_Spal = 12.72;
///@brief 26Al spallation reference production rate.
///@details Spallation (a/g/yr) calibrated 26Al production rate of Lifton et al. (2014), reported in Borchers et al. (2016)..
float Po_26Al_Spal = 28.61;
///@brief 36Cl spallation of Ca reference production rate.
///@details Spallation of calcium (a/g(Ca)/yr) calibrated 36Cl production rate of Lifton et al. (2014), reported in Borchers et al. (2016)..
float Po_36Cl_Ca_Spal = 56.61;
///@brief 36Cl spallation of K reference production rate.
///@details Spallation of pottasium (a/g(K)/yr) calibrated 36Cl production rate for Sf scaling scheme of Lifton et al. (2014), reported in Borchers et al. (2016).
float Po_36Cl_K_Spal = 153.95;
///@brief 10Be muogneic reference produciton rate.
///@details Total muogenic production rate (a/g/yr) following Balco (2017).
float Po_10Be_Muon_Fast = 0.084;
float Po_10Be_Muon_Slow = 0;
///@brief 14C muogneic reference produciton rate.
///@details Total muogenic production rate (a/g/yr) following Heisinger et al. (2002)
float Po_14C_Muon_Fast = 3.34;
float Po_14C_Muon_Slow = 0.44;
///@brief 26Al muogneic reference produciton rate.
///@details Total muogenic production rate (a/g/yr) following Balco (2017) approximation
float Po_26Al_Muon_Fast = 0.761;
float Po_26Al_Muon_Slow = 0.;
///@brief 36Cl muogneic reference produciton rate.
///@details Total muogenic production rate (a/g/yr)
float Po_36Cl_Muon_Fast = 3.34;
float Po_36Cl_Muon_Slow = 0.44;
// Attenuation Lengths
///@brief Spallogenic attenuation length (kg/m^2) following Balco et al (2008).
float Lamb_Spal = 1600.0;
// might need different lambda values for different productions?
///@brief Muogenic attenuation length (kg/m^2) (for 10Be) following Balco (2017).
float Lamb_Muon = 25010;
///@brief Muogenic attenuation length (kg/m^2) for 26Al following Balco (2017)
///float Lamb_26Al_Muon = 24170.0;
// density of rock and sea water respectively
float rho_r = 2600.;
float rho_w = 1035.;
//Decay length scales
float z_rs = Lamb_Spal/rho_r; //Decay length scale chalk spallation
float z_ws = Lamb_Spal/rho_w; //Decay length scale sea water spallation
float z_rm = Lamb_Muon/rho_r; //Decay length scale chalk muons
float z_wm = Lamb_Muon/rho_w; //Decay length scale sea water muons
///@brief Half life of 10Be (Korschineck et al. 2010).
float HalfLife_10Be = 1.387*pow(10.,6.);
float Lambda_10Be = log(2.)/HalfLife_10Be;
///@brief Half life of 14C (ref).
float HalfLife_14C = 5730.;
float Lambda_14C = log(2.)/HalfLife_14C;
///@brief Half life of 26Al (ref).
float HalfLife_26Al = 7.05*pow(10.,5.);
float Lambda_26Al = log(2.)/HalfLife_26Al;
///@brief Half life of 36Cl (ref).
float HalfLife_36Cl = 3.01*pow(10.,5.);
float Lambda_36Cl = log(2.)/HalfLife_36Cl;
//float Po_Spal = 4.34; //Spallation (a/g/yr) to match Regard et al (2012)
//float Po_Fast = 0.039; //Fast Muons (a/g/yr) to match Regard et al (2012)
//float Po_Slow = 0.012; //Slow Muons (a/g/yr) to match Regard et al (2012)
//float Lamb_Fast = 15000.0; //Slow Muons (kg/m^2) to match Regard et al (2012)
//float Lamb_Slow = 43200.0; //Fast Muons (kg/m^2) to match Regard et al (2012)
//float z_rfast = Lamb_Fast/rho_r; //Decay length scale chalk fast muons , to match Regard et al (2012)
//float z_wfast = Lamb_Fast/rho_w; //Decay length scale sea water fast muons, to match Regard et al (2012)
//float z_rslow = Lamb_Slow/rho_r; //Decay length scale chalk slow muons, to match Regard et al (2012)
//float z_wslow = Lamb_Slow/rho_w; //Decay length scale sea water slow muons, to match Regard et al (2012)
#endif