update the adiabatic time for both SNIa and SNII and make it work for comobile integr.

In order to make the adiabatic time for both SNIa and SNII, we introduce SNIaThermalTime
and SNIIThermalTime in the SphP structure.
These thwo times are firstly initalized to -1. In chimie.c, once a SN explode in the
neighborhood of the particles, it records the current time. In cooling.c,
we compute the time difference between the current time and SNIIThermalTime and
SNIaThermalTime to decide if the cooling is stopped or not.