/*widla determine Lya line width for radiation pressure */
double widla(double tauin, 
	  double tauout, 
	  double a, 
	  double vth);