Martens_ea_NatComm_SourceApp_MCMC.txt

Martens_ea_NatComm_SourceApp_MCMC.txt

@@ -11,7 +11,7 @@
 % Sea; 10. Central Ocean. % This script is based on the MCMC method published 
 % by Andersson et al. (ES&T, DOI: 10.1021/es503855e).
 
-% CITATION: Martens et al. 2022, in review in Nat. Comm.
+% CITATION: Martens et al. 2022, Nat. Comm.
 clear;clc;clf;close all;
 %% Load the data
 % Here an example data file is loaded from an Excel file
@@ -90,28 +90,28 @@ else
 end
     
 % Here the age-corrected mean D14C values of the SurfSoil and/or Peat end
-% members are calculated for each sample according the equation provided by Bröder et al. 2018
+% members are calculated for each sample based on the equation provided by Bröder et al. 2018
 fMC_F1(j,1)=(Shelfdata(j,6)+1000)/(1000*exp(1/8033));
 fMC_F1_transp(j,1)=fMC_F1(j,1)*exp(-(CSTT_F1(j))/8033);
 D14C_F1(j,1)=(fMC_F1_transp(j,1)*exp(1/8033)-1)*1000;
 M_A_av=D14C_F1(j,1);
 
 % Here the age-corrected standard deviation of the SurfSoil and/or Peat end
-% member mean value are calculated for each sample according the equation provided by Bröder et al. 2018
+% member mean value are calculated for each sample based on the equation provided by Bröder et al. 2018
 fMC_F1(j,2)=((Shelfdata(j,6)-Shelfdata(j,7))+1000)/(1000*exp(1/8033));
 fMC_F1_transp(j,2)=fMC_F1(j,2)*exp(-(CSTT_F1(j))/8033);
 D14C_F1(j,2)=(fMC_F1_transp(j,2)*exp(1/8033)-1)*1000;
 M_A_std=M_A_av-D14C_F1(j,2); 
 
 % Here the age-corrected mean D14C value of the ICD end member is
-% calculated according the equation provided by Bröder et al. 2018
+% calculated based on the equation provided by Bröder et al. 2018
 fMC_F2(j,1)=(Shelfdata(j,8)+1000)/(1000*exp(1/8033));
 fMC_F2_transp(j,1)=fMC_F2(j,1)*exp(-(CSTT_F2(j))/8033);
 D14C_F2(j,1)=(fMC_F2_transp(j,1)*exp(1/8033)-1)*1000;
 M_B_av=D14C_F2(j,1);
 
 % Here the age-corrected standard deviation of the ICD mean value is
-% calculated according the equation provided by Bröder et al. 2018
+% calculated based on the equation provided by Bröder et al. 2018
 fMC_F2(j,2)=((Shelfdata(j,8)-Shelfdata(j,9))+1000)/(1000*exp(1/8033));
 fMC_F2_transp(j,2)=fMC_F2(j,2)*exp(-(CSTT_F2(j))/8033);
 D14C_F2(j,2)=(fMC_F2_transp(j,2)*exp(1/8033)-1)*1000;