  The thermodynamic Bethe ansatz equations for the N = 6 Coqblin-Schrieffer
model with crystal fields have been solved numerically. The realistic case of
three Kramers doublets with arbitrary splittings has been studied for the first
time. The specific heat has been calculated for representative combinations of
the ionic energy splittings providing ample material for comparison with
experimental results for Cerium impurities and compounds.
