/*tav computes average of old and new optical depths for new scale at end of iter */
void tav(EmLine * t , 
  double f1);