/*nphot one of Kevin Kolk's quantum heating routines for grains */
void nphot(float umin, 
	  float umax, 
	  float *xphot, 
	  long int ndust1);