tfix 26 muon prod values, fix isotope numbers - 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 3df2de3130aac5b7b292e874004f38e845584362
 (DIR) parent 74c25ccd8c430a3d4b9c31412ca98ccb66da01cc
 (HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
       Date:   Thu,  3 Dec 2015 16:09:58 +0100
       
       fix 26 muon prod values, fix isotope numbers
       
       Diffstat:
         M matlab/generate_plots.m             |      24 ++++++++++++------------
         M matlab/m_pakke2014maj11/gCosmoLong… |       6 ++++--
       
       2 files changed, 16 insertions(+), 14 deletions(-)
       ---
 (DIR) diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m
       t@@ -874,15 +874,15 @@ html = ['\n' ...
            '            <td>' num2str(be_conc/1000.) ' atoms/g</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
       -    '            <td><sup>10</sup>Al concentration</td>\n'...
       +    '            <td><sup>26</sup>Al concentration</td>\n'...
            '            <td>' num2str(al_conc/1000.) ' atoms/g</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
       -    '            <td><sup>10</sup>C concentration</td>\n'...
       +    '            <td><sup>14</sup>C concentration</td>\n'...
            '            <td>' num2str(c_conc/1000.) ' atoms/g</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
       -    '            <td><sup>10</sup>Ne concentration</td>\n'...
       +    '            <td><sup>21</sup>Ne concentration</td>\n'...
            '            <td>' num2str(ne_conc/1000.) ' atoms/g</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
       t@@ -890,15 +890,15 @@ html = ['\n' ...
            '            <td>' num2str(be_uncer) ' %%</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
       -    '            <td><sup>10</sup>Al conc. uncertainty</td>\n'...
       +    '            <td><sup>26</sup>Al conc. uncertainty</td>\n'...
            '            <td>' num2str(al_uncer) ' %%</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
       -    '            <td><sup>10</sup>C conc. uncertainty</td>\n'...
       +    '            <td><sup>14</sup>C conc. uncertainty</td>\n'...
            '            <td>' num2str(c_uncer) ' %%</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
       -    '            <td><sup>10</sup>Ne conc. uncertainty</td>\n'...
       +    '            <td><sup>21</sup>Ne conc. uncertainty</td>\n'...
            '            <td>' num2str(ne_uncer) ' %%</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
       t@@ -910,15 +910,15 @@ html = ['\n' ...
            '            <td>' num2str(be_prod_spall/1000.) ' atoms/g/yr</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
       -    '            <td><sup>10</sup>Al production (spallation)</td>\n'...
       +    '            <td><sup>26</sup>Al production (spallation)</td>\n'...
            '            <td>' num2str(al_prod_spall/1000.) ' atoms/g/yr</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
       -    '            <td><sup>10</sup>C production (spallation)</td>\n'...
       +    '            <td><sup>14</sup>C production (spallation)</td>\n'...
            '            <td>' num2str(c_prod_spall/1000.) ' atoms/g/yr</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
       -    '            <td><sup>10</sup>Ne production (spallation)</td>\n'...
       +    '            <td><sup>21</sup>Ne production (spallation)</td>\n'...
            '            <td>' num2str(ne_prod_spall/1000.) ' atoms/g/yr</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
       t@@ -926,15 +926,15 @@ html = ['\n' ...
            '            <td>' num2str(be_prod_muons/1000.) ' atoms/g/yr</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
       -    '            <td><sup>10</sup>Al production (muons)</td>\n'...
       +    '            <td><sup>26</sup>Al production (muons)</td>\n'...
            '            <td>' num2str(al_prod_muons/1000.) ' atoms/g/yr</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
       -    '            <td><sup>10</sup>C production (muons)</td>\n'...
       +    '            <td><sup>14</sup>C production (muons)</td>\n'...
            '            <td>' num2str(c_prod_muons/1000.) ' atoms/g/yr</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
       -    '            <td><sup>10</sup>Ne production (muons)</td>\n'...
       +    '            <td><sup>21</sup>Ne production (muons)</td>\n'...
            '            <td>' num2str(ne_prod_muons/1000.) ' atoms/g/yr</td>\n'...
            '        </tr>\n'...
            '        <tr>\n'...
 (DIR) diff --git a/matlab/m_pakke2014maj11/gCosmoLongsteps.m b/matlab/m_pakke2014maj11/gCosmoLongsteps.m
       t@@ -96,6 +96,8 @@ c3 = 0.375;
        % ratios between fast muon capture and negative muons
        nm_rat10 = 0.1070/(0.1070+0.0940);
        fm_rat10 = 0.0940/(0.1070+0.0940);
       +nm_rat26 = 0.7/(0.7+0.6);
       +fm_rat26 = (0.6/(0.7+0.6));
        
        %>>>BHJ: To be used in analytical expressions
        % L_spal = Tau_spal/rho; %Decay depth, exp(-z/L)
       t@@ -159,8 +161,8 @@ end
        if ~isempty(fixed_stuff.al_prod_spall) && ...
                ~isempty(fixed_stuff.al_prod_muons)
            P26_top_spal = fixed_stuff.al_prod_spall;
       -    P26_top_nm = nm_rat10*fixed_stuff.al_prod_muons;
       -    P26_top_fm = fm_rat10*fixed_stuff.al_prod_muons;
       +    P26_top_nm = nm_rat26*fixed_stuff.al_prod_muons;
       +    P26_top_fm = fm_rat26*fixed_stuff.al_prod_muons;
        else
            P26_top_spal=31.1e3; %atoms/kg/yr
            P26_top_nm=0.7e3; %atoms/kg/yr