function rl = Alb_Cox(mV,Tc,Td,c1,c2,c3)
%{
    Albrecht-Cox IR flux correction with three coefficients
    mV = thermopile voltage in mV
    Tc, Td = case & dome temperature in °K
%}

sig_sb = 5.67e-8;    % Stefan Boltzmann constant

rl = (mV/c1) + c2*sig_sb*Tc.^4 - c3*sig_sb*(Td.^4 - Tc.^4);

end