/*ufunct one of K. Volk's quantum heating routines for grains */
double ufunct(double temp, 
  long int ndust1);