Cd computation enabled.

The variable cl(2,3) that contains may contain the force coefficients in
each row. x in row 1, y in row 2, and z in row 3.
Save out if cl->infty