def write_model_map(self, model_name, name=None):
        """Save the counts model map to a FITS file.

        Parameters
        ----------
        model_name : str
            String that will be append to the name of the output file.
        name : str
            Name of the component.

        Returns
        -------

        """
        maps = [c.write_model_map(model_name, name) for c in self.components]

        outfile = os.path.join(self.workdir,
                               'mcube_%s.fits' % (model_name))

        mmap = Map.from_geom(self.geom)
        for m in maps:
            mmap.coadd(m)
        mmap.write(outfile, overwrite=True, conv='fgst-ccube')
        return [mmap] + maps