function [B] = loadBOG()
%LOADBOG Summary of this function goes here
%   Detailed explanation goes here
filepath = 'C:\Users\bwerb\Documents\Spray_crossovers\BOG\O2_TCO2_ALK_pH_loadcheck.xlsm';
opts = detectImportOptions(filepath);
B = readtable(filepath,opts);
B.DATE_TIME.TimeZone = "UTC";

% var names to match spray data
B = renamevars(B, 'pH', 'pH20atm'); 
B = renamevars(B, 'PRESSURE', 'pres'); 
B = renamevars(B, 'TMP', 'tc'); 
B = renamevars(B, 'SAL', 'psal');
B = renamevars(B, 'O2', 'doxy');

trex = CO2SYSv3(B.TCO2, B.pH20atm, 2, 3, B.psal, 20, 25, B.pres, 0, B.SIO4, B.PO4, 0, 0, 1, 10, 1, 2, 2); % convert from pH20atm to pH25atm


B.pH25atm = reshape(trex(:,20), size(B.pres));
B.pH25atm(B.pH25atm == -999) = NaN;

trex = CO2SYSv3(B.TCO2, B.pH20atm, 2, 3, B.psal, 20, B.tc, B.pres, 0, B.SIO4, B.PO4, 0, 0, 1, 10, 1, 2, 2); % convert from pH20atm to pH25atm
B.pHin = reshape(trex(:,20), size(B.pres));
B.pHin(B.pHin == -999) = NaN;

end