/*rrfit calculates rates of radiative recombination for all ions */
void rrfit(long int iz, 
	  long int in, 
	  double t, 
	  float *r);