/*PutHetCol save heating, cooling, and temperature in stack for numerical derivatives */
void PutHetCol(double te, 
  double htot, 
  double ctot);