function [Vct, Vcp, VrsHClCorrection]  = vrsHClCorr(TdegC, PdBar)

R = 8.3145;
F = 96485;

tDegK = TdegC + 273.15;

meanGamma = -0.000113 .* TdegC + 0.9067;
deltaV0 = 16.38 + (0.0896 .* TdegC) - (0.001264 .* TdegC.^2);
Vct = ((-R .* tDegK) ./ F) .* log((meanGamma .* 0.01) .^ 2);
Vcp = (deltaV0 ./ 23061) .* 0.0242 .* ((PdBar - 10) ./ 10) .* (1/1.01);

VrsHClCorrection = Vct + Vcp;