def diffusive_heat_transport(self):
phi = np.deg2rad(self.lat)<EOL>phi_stag = np.deg2rad(self.lat_bounds)<EOL>D = self.param['<STR_LIT:D>']<EOL>T = np.squeeze(self.Ts)<EOL>dTdphi = np.diff(T) / np.diff(phi)<EOL>dTdphi = np.append(dTdphi, <NUM_LIT:0.>)<EOL>dTdphi = np.insert(dTdphi, <NUM_LIT:0>, <NUM_LIT:0.>)<EOL>return (<NUM_LIT>*-<NUM_LIT:2>*np.math.pi*np.cos(phi_stag)*const.a**<NUM_LIT:2>*D*dTdphi)<EOL>
Compute instantaneous diffusive heat transport in unit :math:`\\textrm{PW}`
        on the staggered grid (bounds) through calculating:

        .. math::

            H(\\varphi) = - 2 \pi R^2 cos(\\varphi) D \\frac{dT}{d\\varphi}
                        \\approx - 2 \pi R^2 cos(\\varphi) D \\frac{\Delta T}{\Delta \\varphi}

        :rtype: array of size ``np.size(self.lat_bounds)``

        THIS IS DEPRECATED AND WILL BE REMOVED IN THE FUTURE. Use the diagnostic
        ``heat_transport`` instead, which implements the same calculation.