code: def write_model_map(self, model_name, name=None):
          
          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
docs: 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
          -------