function [Estimates,Uncertainties] = ESPER_TSO(lat, lon, depth, salinity, tc, doxy, sdn);

% input can be in matricies; will vectorize them here
% 
% 


DesireVar = [1, 2, 3, 4, 5, 6]; %TA, DIC, pH, phosphate, nitrate, silicate
OutCoords = horzcat(lon(:), lat(:), depth(:)); %lon, lat, depth
PredictorTypes = [1 2 6]; % PSAL, TEMP, DOXY_ADJ, 
MeasEsper = [salinity(:), tc(:), doxy(:)];
refyear = year(sdn(:))+month(sdn(:))/12; % reference year for OA adjustment
Equations = 7; % S, T, O2




% ESPER-LIR for TA, DIC, pH, and Nitrate
disp('Starting ESPER_Mixed...');
[Estimates,Uncertainties] = ESPER_Mixed(DesireVar, OutCoords, MeasEsper,...
    PredictorTypes, 'Equations', Equations, 'EstDates', refyear); 
disp('Finished ESPER_LIR');

 return; 