Improvement of Transformer Saturation Modeling for Electromagnetic Transient Programs