/*HydroEinstA calculates Einstein A's from  osillator strengths*/
double HydroEinstA(long int lower, 
	  long int iupper);