How can make dropolet be exothermic
 

In lagrangian multiphase problem, how can make droplet be exothermic when it evaporates and burns in the core flow.

I try to use drheat.f (defined heat by EMENTH ), however,no matter what expression to be used , the droplet always be endothermic as long as it evaporates.

Should I use sorent.f ?