def inferred_heat_transport(self):
        """Calculates the inferred heat transport by integrating the TOA
        energy imbalance from pole to pole.

        The method is calculating

        .. math::

            H(\\varphi) = 2 \pi R^2 \int_{-\pi/2}^{\\varphi} cos\phi \ R_{TOA} d\phi

        where :math:`R_{TOA}` is the net radiation at top of atmosphere.


        :return: total heat transport on the latitude grid in unit :math:`\\textrm{PW}`
        :rtype: array of size ``np.size(self.lat_lat)``

        :Example:

            .. plot:: code_input_manual/example_EBM_inferred_heat_transport.py
                :include-source:

        """
        phi = np.deg2rad(self.lat)
        energy_in = np.squeeze(self.net_radiation)
        return (1E-15 * 2 * np.math.pi * const.a**2 *
                integrate.cumtrapz(np.cos(phi)*energy_in, x=phi, initial=0.))