/* sum up the charge transfer heating */
double SumCTHeat(void);

/*HCTion H charge transfer ionization, H+ + A => H + A+ */
double HCTIon(long int ion, 
  long int nelem);

/*HCTRecom H charge transfer recombination, H + A+ => H+ + A */
double HCTRecom(long int ion, 
  long int nelem);