docs: def lon_lat_trt_pmf(matrices):
      """
      Fold full disaggregation matrices to lon / lat / TRT PMF.
  
      :param matrices:
          a matrix with T submatrices
      :returns:
          3d array. First dimension represents longitude histogram bins,
          second one latitude histogram bins, third one trt histogram bins.
      """
      res = numpy.array([lon_lat_pmf(mat) for mat in matrices])
      return res.transpose(1, 2, 0)