Systems that have more than one conserved quantity (i.e. energy plus momentum, density etc.), can exhibit quite interesting temperature profiles in non-equilibrium stationary states. I will present some numerical experiment and mathematical result. I will also expose some other connected problems, always concerning thermal boundary conditions in hydrodynamic limits.