This is LK with S=2 and correct boundary current conservation & This is GS with gamma & kappa = 0, and T eqn in (m_b+m_f) form.

This is with gamma=0.5 and T eqn in (m_b+m_f) form
All S=2, iteration = 10000, deltaT=0.001
gamma(0.5),m=1,
gamma(2),m=1,
gamma(3),m=1,
gamma(5),m=1,

gamma(1),m=0.5,
gamma(2),m=0.5,
gamma(3),m=0.5,
gamma(5),m=0.5,