Why is running cosmomc 2010 with cls = 3 (TT data only) apparently a problem ?

Use of Cobaya. camb, CLASS, cosmomc, compilers, etc.
Post Reply
Patrice Okouma
Posts: 19
Joined: November 05 2009
Affiliation: University of Cape Town

Why is running cosmomc 2010 with cls = 3 (TT data only) appa

Post by Patrice Okouma » August 04 2011

Dear all ,
I want to run cosmomc (ver. 2010) using WMAP TT data only.
From http://cosmocoffee.info/viewtopic.php?t ... f113780d3e ,
I need to comment out WMAP_LnLike_TE contribution at the bottom of cmbdata.f90 (in old cosmomc version).

I am using cosmomc 2010 and the section to comment out is not used anymore.
So I went into WMAP_7yr_options.F90 (inside likelihood code) where I set :
logical :: use_TE = .false.
logical :: use_lowl_pol = .false. ! include TE,EE,BB pixel likelihood for l<24
logical :: use_TE = .false. ! include MASTER TE in likelihood

Now, when ran, cosmomc gets frozen at :
Random seeds: 2982, 15661 rand_inst: 0
Computing tensors: F
Doing CMB lensing: T
lmax = 2100
Number of C_ls = 3

Could you plz suggest me a way out of this ?

Thanking you,
Patrice

Patrice Okouma
Posts: 19
Joined: November 05 2009
Affiliation: University of Cape Town

Why is running cosmomc 2010 with cls = 3 (TT data only) appa

Post by Patrice Okouma » August 07 2011

Modifying WMAP_7yr_options.F90 (inside likelihood) as follows does the trick (+ using a cosmomc ver May 2010 in my case) :

!---------------------------------------------------
! various likelihood options
! change these to include/ exclude various likelihood aspects
!---------------------------------------------------
logical :: use_lowl_TT = .false. ! include TT pixel likelihood, for l<=lowl_max
logical :: use_lowl_pol = .true. ! include TE,EE,BB pixel likelihood for l<24
logical :: use_TT = .true. ! include MASTER TT in likelihood
logical :: use_TT_beam_ptsrc = .true. ! include beam/ptsrc errors
logical :: use_TE = .false. ! include MASTER TE in likelihood

logical :: use_gibbs_pol_cleaning = .false.

Patrice

Post Reply