/*level3 compute three level atom, 10, 21, and 20 are line */
void level3(EmLine * t10, 
  EmLine * t21, 
  EmLine * t20);