steadyEplot.m

function [Ec,conc,h] = steadyEplot(p,E)

% Supplemental material for:
% Article Title: Quantifying erosion rates and weathering pathways that maximize soil organic carbon storage
% Journal: Biogeochemistry
% Authors: Joshua J. Roering*, Brooke D. Hunter, Ken L. Ferrier, Oliver A. Chadwick, Kyungsoo Yoo, Adrian A. Wackett, Peter C. Almond, Lucas Silva, and A. Mark Jelline
% Corresponding author: Josh Roering, Dept of Earth Science, University of Oregon, jroering@uoregon.edu

%this script summarizes and plots how mineral concentrations vary with 
%erosion rate.  all terms are non-dimensional according to Ferrier & West
%2017. the script takes in the p structure created by lag_parameters. 

%Example call
%E=[0.01:0.01:0.99];    %creates vector with range of E values
%p=lag_parameters;[Ec,conc,h]=steadyEplot(p,E); 

%how many non-dim erosion rates will be calculated
Enum=length(E);

%create a loop to cycle through E and create conc data matrix
for i=1:Enum
    B=steadyq_nd(E(i),p);
    Ec(i)=B;
    conc(i,:)=Ec(i).qs_nd;
    h(i)=Ec(i).H_nd;
end

plot(E,conc); legend(p.mineralname);
figure; plot(E,conc(:,6));

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%MIT License

% Copyright (c) 2023 Joshua J. Roering
% 
% Permission is hereby granted, free of charge, to any person obtaining a copy
% of this software and associated documentation files (the "Software"), to deal
% in the Software without restriction, including without limitation the rights
% to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
% copies of the Software, and to permit persons to whom the Software is
% furnished to do so, subject to the following conditions:
% 
% The above copyright notice and this permission notice shall be included in all
% copies or substantial portions of the Software.
% 
% THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
% IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
% FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
% AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
% LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
% OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
% SOFTWARE.