/*expion VERY simple two level cooling, do de-excitation or transfer */
double expion(double bltz, 
  double abund);