tuse production rates from web interface, muons may be wrong - cosmo - front and backend for Markov-Chain Monte Carlo inversion of cosmogenic nuclide concentrations
(HTM) git clone git://src.adamsgaard.dk/cosmo
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) README
(DIR) LICENSE
---
(DIR) commit b1e24d4743bd7fb10ec52d6bc2b8ba5aabc61436
(DIR) parent 8bfbd1b6e11b787ba4356204901c7e4e891027a0
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Wed, 25 Nov 2015 09:02:32 +0100
use production rates from web interface, muons may be wrong
Diffstat:
M matlab/m_pakke2014maj11/gCosmoLong… | 27 +++++++++++++++++++--------
1 file changed, 19 insertions(+), 8 deletions(-)
---
(DIR) diff --git a/matlab/m_pakke2014maj11/gCosmoLongsteps.m b/matlab/m_pakke2014maj11/gCosmoLongsteps.m
t@@ -102,9 +102,11 @@ c3 = 0.375;
%10Be production
% Input fra Kasper
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% - Start
-P10_top_spal=5.33e3; %atoms/kg/yr
-P10_top_nm=0.106e3; %atoms/kg/yr
+%P10_top_spal=5.33e3; %atoms/kg/yr
+%P10_top_nm=0.106e3; %atoms/kg/yr
P10_top_fm=0.093e3; %atoms/kg/yr
+P10_top_spal = fixed_stuff.be_prod_spall;
+P10_top_nm = fixed_stuff.be_prod_muons;
%Reference values for Kasper
%P10_top_spal=5.33e3; %atoms/kg/yr
t@@ -137,10 +139,13 @@ P10_top_fm=0.093e3; %atoms/kg/yr
% Input fra Kasper
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% - Start
-P26_top_spal=31.1e3; %atoms/kg/yr
-P26_top_nm=0.7e3; %atoms/kg/yr
+%P26_top_spal=31.1e3; %atoms/kg/yr
+%P26_top_nm=0.7e3; %atoms/kg/yr
P26_top_fm=0.6e3; %atoms/kg/yr
+P26_top_spal = fixed_stuff.al_prod_spall;
+P26_top_nm = fixed_stuff.al_prod_muons;
+
%Reference values for Kasper
%P26_top_spal=31.1e3; %atoms/kg/yr
%P26_top_nm=0.7e3; %atoms/kg/yr
t@@ -162,10 +167,13 @@ P26_top_fm=0.6e3; %atoms/kg/yr
% P26_total = (P26_spal + P26_nm + P26_fm);
%21Ne production
-P21_top_spal=20.8e3; %atoms/kg/yr
-P21_top_nm=0.4e3; %atoms/kg/yr
+%P21_top_spal=20.8e3; %atoms/kg/yr
+%P21_top_nm=0.4e3; %atoms/kg/yr
P21_top_fm=0.35e3; %atoms/kg/yr
+P21_top_spal = fixed_stuff.ne_prod_spall;
+P21_top_nm = fixed_stuff.ne_prod_muons;
+
% P21_spal = P21_top_spal*exp(-z*rho/Tau_spal);
% P21_nm = P21_top_nm*exp(-z*rho/Tau_nm);
% P21_fm = P21_top_fm*exp(-z*rho/Tau_fm);
t@@ -173,10 +181,13 @@ P21_top_fm=0.35e3; %atoms/kg/yr
% P21_total = (P21_spal + P21_nm + P21_fm);
%14C production
-P14_top_spal=14.6e3; %atoms/kg/yr
-P14_top_nm=2.3e3; %atoms/kg/yr
+%P14_top_spal=14.6e3; %atoms/kg/yr
+%P14_top_nm=2.3e3; %atoms/kg/yr
P14_top_fm=2.1e3; %atoms/kg/yr
+P14_top_spal = fixed_stuff.c_prod_spall;
+P14_top_nm = fixed_stuff.c_prod_muons;
+
% P14_spal = P14_top_spal*exp(-z*rho/Tau_spal);
% P14_nm = P14_top_nm*exp(-z*rho/Tau_nm);
% P14_fm = P14_top_fm*exp(-z*rho/Tau_fm);