/*opinPL add opacity of individual species, with stimulated emission, for power law opacity */
void opinPL(long int n, 
  long int ipZ, /* ipZ=0 for h, 1 for he, etc */
  float abundance, 
  double b);