/*abscf convert gf into absorption coefficient */
double abscf(double gf, 
	  double enercm, 
	  double gl);