if (Data_Initial == 1)
    Data = load('TS2.dat');

    Tc = 25;               % Temp in Celcius 

    pH = 8.886;            % us: 8.886   
    NaHCO3 = 0.015;        % [NaHCO3] in mol/L  us: 0.015   
    d18O_w = -3.0309;      % average d18O of H2O (VSMOW)
    d13C_SB = -2.7831;     % d13C of NaHCO3 used in experiments (VPDB)

    % 1st non-NaN d18O_BaCO3 (VPDB)
    n =  find(~isnan(Data(:,5)), 1);
    d18O_1st = Data(n, 5);
        % Converting to VSMOW
        d18O_dcp_t0 = 1.03091*(d18O_1st)+30.91;     

    % 1st non-NaN D47_BaCO3 data 
    n =  find(~isnan(Data(:,9)), 1);
    D47_DCP = Data(n, 9);   

    D47_co2 = 0.6;         % value essentially irrelevant 

    % Time Zero for the model
    T0 = 210; % sec

    % CA concentration (mol/L)
    CA = 0;
end

if(Data_Initial == 2)
    Data = load('TS2_CAmn.dat');
    Tc = 25;               % Temp in Celcius 

    pH = 8.886;            % us: 8.886   ss18: 9.2-9.4 (9.3)
    NaHCO3 = 0.015;        % [NaHCO3] in mol/L  us: 0.015   ss18: 0.0357
    d18O_w = -3.077;      % average d18O of H2O (VSMOW)
    d13C_SB = -2.7831;     % d13C of NaHCO3 used in experiments (VPDB)

    % 1st non-NaN d18O_BaCO3 data (VPDB)
    n =  find(~isnan(Data(:,5)), 1);
    d18O_1st = Data(n, 5);
        % Converting to VSMOW
        d18O_dcp_t0 = 1.03091*(d18O_1st)+30.91;     

    % 1st non-NaN D47_BaCO3 data 
    n =  find(~isnan(Data(:,9)), 1);
    D47_DCP = Data(n, 9);    

    D47_co2 = 0.6;          % value essentially irrelevant 

    % Time Zero for model
    T0 = 220;  

    % CA concentration (mol/L)
    CA = 3.7E-9;
end

if(Data_Initial == 3)
    Data = load('TS2_CA2.dat');
    Tc = 25;               % Temp in Celcius 

    pH = 8.886;            % us: 8.886   ss18: 9.2-9.4 (9.3)
    NaHCO3 = 0.015;        % [NaHCO3] in mol/L  us: 0.015   ss18: 0.0357
    d18O_w = -3.13;      % average d18O of H2O (VSMOW)
    d13C_SB = -2.7831;     % d13C of NaHCO3 used in experiments (VPDB)

    % 1st non-NaN d18O_BaCO3 data (VPDB)
    n =  find(~isnan(Data(:,5)), 1);
    d18O_1st = Data(n, 5);
        % Converting to VSMOW
        d18O_dcp_t0 = 1.03091*(d18O_1st)+30.91;     

    % 1st non-NaN D47_BaCO3 data 
    n =  find(~isnan(Data(:,9)), 1);
    D47_DCP = Data(n, 9);    

    D47_co2 = 0.6;          % value essentially irrelevant 

    % Time Zero for the model
    T0 = 160;  

    % CA concentration (mol/L)
    CA = 9.3E-9;
end

if(Data_Initial == 4)
    Data = load('TS2_CA3.dat');
    Tc = 25;               % Temp in Celcius 

    pH = 8.886;            % us: 8.886   ss18: 9.2-9.4 (9.3)
    NaHCO3 = 0.015;        % [NaHCO3] in mol/L  us: 0.015   ss18: 0.0357
    d18O_w = -3.10;        % average d18O of H2O (VSMOW)
    d13C_SB = -2.7831;     % d13C of NaHCO3 used in experiments (VPDB)


    % 1st non-NaN d18O_BaCO3 data (VPDB)
    n =  find(~isnan(Data(:,5)), 1);
    d18O_1st = Data(n, 5)-0.5;
        % Converting to VSMOW
        d18O_dcp_t0 = 1.03091*(d18O_1st)+30.91;     

    % 1st non-NaN D47_BaCO3 data 
    n =  find(~isnan(Data(:,9)), 1);
    D47_DCP = Data(n, 9)-0.018;    

    D47_co2 = 0.6;          % value essentially irrelevant 

    % Time Zero for model
    T0 = 190;  

    % CA concentration (mol/L)
    CA = 1.9E-8;
end

if(Data_Initial == 5)
    Data = load('SS18.dat');
    Tc = 25;               % Temp in Celcius 

    pH = 9.3;              % ss18: 9.2-9.4 (9.3)
    NaHCO3 = 0.0357;       % [NaHCO3] in mol/L  ss18: 0.0357
    d18O_w = -0.65;        % average d18O of H2O (VSMOW)

    n =  find(~isnan(Data(:,3)), 1);
    d13C_SB = Data(n, 3);
                           % d13C of NaHCO3 not reported, d13C of BaCO3
                           % changed over time. 
                           % Taking the first data-point in their time-series
                           % used as d13C_SB
    % or 
    % d13C_SB = mean(Data(:,3));   
                           % this selection won't make much difference...

    % 1st non-NaN d18O_BaCO3 data (VPDB)
    n =  find(~isnan(Data(:,5)), 1);
    d18O_1st = Data(n, 5);
        % Converting to VSMOW
        d18O_dcp_t0 = 1.03091*(d18O_1st)+30.91;     

    D47_co2 = 0.6;         % values almost irrelevant 

    % Finding 1st non-NaN D47_BaCO3 data (VPDB)
    n =  find(~isnan(Data(:,9)), 1);
    D47_DCP = Data(n, 9);

    % Time Zero for the model
    T0 = 50;  

    % CA concentration (mol/L)
    CA = 0;
end