There is no code that accounts for temperature dependence of heat capacity, at least not any that's been implemented in a general way using the physconst module which is responsible for the values of properties which have traditionally been considered to be constant. However, this module now contains arrays for several constants, the heat capacity of air among them, to allow them to vary spatially and temporally. This was implemented to allow the properties of air to vary with composition in the WACCM-X model which has a top far above the homosphere. You could use a similar mechanism to allow for spatial/temporal variation of any physical property.