CosmoMC modifications to include the actlite_3yr_v1 likelihood. The parameters reported in Dunkley et al. 2013 and Sievers et al. 2013 are obtained with the CosmoMC January 2012 version. =============================================================================== ------------------------------------------------------------------------------- 1) source/cmbdata.F90 ------------------------------------------------------------------------------- < logical :: init_MAP = .true. > logical :: init_MAP = .true. > logical :: init_highell = .true. < if (aname == 'MAP' .or. aname == 'WMAP') then < aset%name = 'WMAP' < datasets(num_datasets) = aset < return > if (aname == 'MAP' .or. aname == 'WMAP') then > aset%name = 'WMAP' > datasets(num_datasets) = aset > return > elseif (aname == 'Highell' ) then > aset%name = 'Highell' > datasets(num_datasets) = aset > return < if (datasets(i)%name == 'WMAP') then < tot(i) = MAPLnLike(szcl) < else > if (datasets(i)%name == 'WMAP') then > tot(i) = MAPLnLike(szcl) > else if (datasets(i)%name == 'Highell') then > tot(i) = HighellLnLike(cl) > else < end module cmbdata > function HighellLnLike(cl) > #ifndef NOHIGHELL > use ACTlite_3yr_like > #endif > real cl(lmax,num_cls_tot),HighellLnLike > real(8), dimension(lmax) :: cl_tt > #ifndef NOHIGHELL > real(8) :: like_act > integer l > if (Init_highell) then > call act_likelihood_init > if (Feedback>0) write(*,*) 'reading ACT cmb data' > Init_highell = .false. > end if > cl_tt(2:lmax) = 0.d0 > do l =2, tt_lmax > cl_tt(l) = cl(l,1)*l*(l+1)/twopi > enddo > like_act = 0.d0 > call act_likelihood_compute(cl_tt,like_act) > HighellLnLike = like_act > write(*,*), ' ACT chi2 = ', 2*HighellLnLike > #else > HighellLnLike==cl(2,1) > #endif > end function > end module cmbdata ------------------------------------------------------------------------------- 2) source/Makefile.F90 ------------------------------------------------------------------------------- < WMAP = "your WMAP likelihood path" > WMAP = "your WMAP likelihood path" > HIGHELL = "your ACT likelihood path" < else < F90FLAGS += -DNOWMAP < endif > else > F90FLAGS += -DNOWMAP > endif > ifneq ($(HIGHELL),) > F90FLAGS += $(IFLAG)$(cfitsio)/include $(IFLAG)$(HIGHELL) > LINKFLAGS += -L$(cfitsio)/lib -L$(HIGHELL) -lcfitsio > OBJFILES += $(HIGHELL)/ACTlite_3yr_like.o \ > else > F90FLAGS += -DNOHighell > endif ------------------------------------------------------------------------------- 3) params.ini ------------------------------------------------------------------------------- > cmb_numdatasets = 2 > cmb_dataset1 = WMAP > cmb_dataset2 = Highell