function [ kd ] = Ctemp_m( Temp )
%
%This routine computes the temperature correction factor
%
%  t = temperature in Centigrade

%  First convert to farenheit

f = (9 * (Temp + 40) / 5) - 40;
if ((f >= 70) && (f <= 1000))
    kd = 0.975 + 0.432e-3*f - 0.115e-5*f^2 + 0.104e-8*f^3 - 0.595e-12*f^4;
else
    error('MATLAB:Ctemp_m', ' Temperature %6.1f is out of range', Temp);
end

