added Tb as a constant 42+273 as per Pieter's note. Check that this is correct. Code does not run without defining a default Tb here, unless we add it as an input