Reactivated writing out the temperature if stream(3) == 1
and the code is running non-isothermal