tuse ratios from same file - cosmo - front and backend for Markov-Chain Monte Carlo inversion of cosmogenic nuclide concentrations
git clone git://src.adamsgaard.dk/cosmo
Log
Files
Refs
README
LICENSE
---
commit 593310ddd57dbd7fdad41c884702054acaf2c7ae
parent a8f727dac7d0609c888c9e8e562ea93e21758f75
Author: Anders Damsgaard 
Date:   Mon,  7 Dec 2015 20:44:50 +0100

use ratios from same file

Diffstat:
  M matlab/m_pakke2014maj11/gCosmoLong… |      13 +++++++------

1 file changed, 7 insertions(+), 6 deletions(-)
---
diff --git a/matlab/m_pakke2014maj11/gCosmoLongsteps.m b/matlab/m_pakke2014maj11/gCosmoLongsteps.m
t@@ -99,11 +99,12 @@ fm_rat10 = 0.0940/(0.1070+0.0940);
 nm_rat26 = 0.7/(0.7+0.6);
 fm_rat26 = (0.6/(0.7+0.6));
 
-% use ratios from 10Be since I (Anders) don't know better...
-nm_rat14 = nm_rat10;
-fm_rat14 = fm_rat10;
-nm_rat21 = nm_rat10;
-fm_rat21 = fm_rat10;
+% unconfirmed ratios
+nm_rat14 = 0.4/(0.4+0.35);
+fm_rat14 = 0.35/(0.4+0.35);
+nm_rat21 = 2.3/(2.3+2.1);
+fm_rat21 = 2.1/(2.3+2.1);
+
 
 %>>>BHJ: To be used in analytical expressions
 % L_spal = Tau_spal/rho; %Decay depth, exp(-z/L)
t@@ -115,7 +116,7 @@ fm_rat21 = fm_rat10;
 % Input fra Kasper
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% - Start
 P10_top_spal=5.33e3; %atoms/kg/yr
-P10_top_nm=0.106e3; %atoms/kg/yr
+1P10_top_nm=0.106e3; %atoms/kg/yr
 P10_top_fm=0.093e3; %atoms/kg/yr
 
 if exist('fixed_stuff.be_prod_spall', 'var') == 1 && ...