Small cleaning of the code including a better organisation of scattering and thermalisation parts within the interaction function; upgrade of the keyword ndiffmax, replaced by nscattmax