correct EgySpec in chimie_apply_thermal_feedback()
add two new flags:
COMPUTE_PHYSICAL_UNITS WITH_ID_IN_DENSITY