ThDDGOXupdatelattice_temp [SECS2D]
Tl = ThDDGOXupdatelattice_temp(mesh,Dnodes,Tl,Tn,Tp,n,p,...
kappa,Egap,tn,tp,...
twn0,twp0,twn1,twp1,n0,p0)