/*hmrate compile molecular rates using Hollenback and Mckee fits */
double hmrate(double a, 
	  double b, 
	  double c);