code: def inferred_heat_transport(self):
          
          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.))
docs: 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: