  We present full-sky maps of the Integrated Sachs-Wolfe effect (ISW) for the
MICE Grand Challenge lightcone simulation up to redshift 1.4. The maps are
constructed in the linear regime using spherical Bessel transforms. We compare
and contrast this procedure against analytical approximations found in the
literature. By computing the ISW in the linear regime we remove the substantial
computing and storage resources required to calculate the non-linear
Rees-Sciama effect. Since the linear ISW is $~10^{2}\times$ larger in $\Delta
T/T$ this has a neglegible impact on the maps produced and only becomes
relevant on scales which are dominated by cosmic microwave background (CMB)
anisotropies. The MICE simulation products have been extensively used for
studies involving current and future galaxy surveys. The availability of these
maps will allow MICE to be used for future galaxy and CMB cross-correlation
studies, ISW reconstruction studies and ISW void-stacking studies probed by
galaxy surveys such as DES, DESI, Euclid and Rubin LSST. The pipeline developed
in this study is provided as a public Python package pyGenISW. This could be
used in future studies for constructing the ISW from existing and future
simulation suites probing vast sets of cosmological parameters and models.
