/*hmole determine populations of hydrogen molecules */
void hmole(float *hmovh1, 
  double h2ovh1);