function tu81 = tu2tu81(year,trit)

% Samar Khatiwala (spk@ldeo.columbia.edu)

l = 0.05575;
y81 = 1981*ones(size(year));
tu81 = trit.*exp(l*(year-y81));
