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