bgc_normalise.m

function nCq_nAq = bgc_normalise(Aq_Cq)
%bgc_normalise Normalise Aq and Cq for any process to unit vector case
% === INPUT ===============================================================
%   Aq_Cq: change in [AT CT] due to biogeochemical process / micromol/kg
% === OUTPUT ==============================================================
% nAq_nCq: input [Aq Cq] normalised such that sqrt(Aq^2 + Cq^2) = 1
% =========================================================================
% --- Written by Matthew P. Humphreys -------------------------------------
%       v1: last updated 2017-06-23
% Citation: Humphreys MP, Daniels CJ, Wolf-Gladrow DA, Tyrrell T, & 
%           Achterberg EP (2017): "On the influence of marine
%           biogeochemical processes over CO2 exchange between the
%           atmosphere and ocean", Marine Chemistry, in press,
%           doi:10.1016/j.marchem.2017.12.006
% =========================================================================

% Unpack input
Aq = Aq_Cq(1);
Cq = Aq_Cq(2);

% Normalise PC
if Cq == 0, nCq = 0;
else, nCq = (abs(Cq)/Cq) * sqrt(Cq^2/(Cq^2+Aq^2));
end %if

% Normalise PA
if Aq == 0, nAq = 0;
else, nAq = (abs(Aq)/Aq) * sqrt(Aq^2/(Cq^2+Aq^2));
end %if

% Generate output
nCq_nAq = [nAq nCq];

end %function bgc_normalise