/*qheat1 one of Kevin Volk's quantum heating routines for grains */
void qheat1(long int ndust1, 
  double *tmean, 
  double ted1);