/*oilevl get OI level population with Ly-beta pumping */
void oilevl(double abund, 
  double *coloi);