/*hdexct compute collisional deexcitation for hydrogen atom, from Vriens and Smeets */
double hdexct(long int ihigh, 
	  long int ilow);