I have been able to locate the source of the error. The weight files generated by lapw2 are supposed to have the actual value of EF in the first line. With the TETRA or GAUSS choices, they do. With the TEMP choice, they are displaced by exactly 30*T. Correspondingly, JOINT uses a wrong Fermi energy.