/*popexc do level population for simple two level atom, no radiative transfer */
double popexc(double omega, 
  double g1, 
  double g2, 
  double a21, 
  double bltz, 
  double abund);