/*widlin determine half width of any line with known optical depths */
double widlin(double tauin, 
	  double tauout, 
	  double a, 
	  double vth);