COSMOMC April 2013: error with 'clik_try_lensing'

[Francesco Montesano]

I managed to compile WMAP 9 and plc-1.0.

I have linked the the directory with Planck data into data/click, defined the correct paths to WMAP and Planck codes in .bashrc and Makefile and managed to compile COSMOMC . Everything has been compiled with ifort, mpiifort and icc V13.1.1.

Then I have tried to run "mpirun -np 2 ./cosmomc params.ini" and got this error:

Number of MPI processes: 2
Random seeds: 14656, 2800 rand_inst: 1
WMAP options (beam TE TT) T T T
Using clik with likelihood file ./data/clik/CAMspec_v6.2TN_2013_02_26_dist.clik
clik_try_lensing(src/clik.c:428)::ForwardError
fopen_err(src/minipmc/io.c:17)::Error -101 (Cannot open file './data/clik/CAMspec_v6.2TN_2013_02_26_dist.clik' (mode "r"))

From other posts like http://cosmocoffee.info/viewtopic.php?t=2083 it's clear that at least in COSMOMC March 2013 the error was not there.
The only modifications that I have done are in the Makefiles to use the correct version of intel compilers and get and install libraries from/to the right places.

Does anyone had similar problem or suggestions?

Thanks,
Francesco

[Antony Lewis]

Does
./data/clik/CAMspec_v6.2TN_2013_02_26_dist.clik

exist?

(test.ini is the default Planck test file)

[Francesco Montesano]

Dear Antony,

stupid me: I did a typo when I did the link (click instead of clik).

Now it runs for a bit, checks the likelihood, add parameters for Planck, throw a

WARNING: zero padding ext cls in LoadFiducialHighLTemplate

does something with a non-linear Pk and than throw the following

[...]
Varying 20 parameters ( 6 slow ( 2 semi-slow), 14 fast ( 0 semi-fast))
key not found : test_likelihood

and terminates without error code. Any idea?

Francesco

[Antony Lewis]

Just set test_likelihood =F (this should be fixed by default in latest download)

[Francesco Montesano]

(this should be fixed by default in latest download)

Do you mean April 2013 version or later?

I have changed line 411 of paramdef.F90

from

Code: Select all

test_likelihood = Ini_read_Logical_file(Ini,'test_likelihood')

to

Code: Select all

test_likelihood = Ini_read_Logical_file(Ini,'test_likelihood', .false.)