pro test
  read_fits_map,'wmap_140.h.fits',map
  readmatrix,aa,'wmap.140.cluster.dat',o=1
  nside=npix2nside(N_ELEMENTS(map))
  map2=fltarr(nside)

  FOR i=0L,10 DO BEGIN
     ang2vec,aa[i,2],aa[i,1],x0,/astro
     query_disc,nside,x0,0.05,ilist,/nested
     map2[ilist]=map[ilist]
  END
  mollview,map2,/nest,/log
end