22 How do I calculate the bulk modulus as a function of temperature?

Use what='mur_lc_t' and the option lmurn=.TRUE.. This approach is rigorously valid only for cubic solids, for anisotropic solid it is only an approximation in which only celldm(1) is changed while the other crystal parameters are kept constant.