c has been redefined such that the hole burning term is 1 / (1 + c^2 |psi|^2), i.e. all the rest of the garbage has been removed. Newton still works perfectly, because this amounts to a simple rescaling. For printing of data, we use the old convention. However, for outputting of c's to files, we're using the new convention.