CosmoCoffee

[Starla Talbot]

Hello,

I am trying to add a new parameter w_a. I followed the information of the ReadMe file of cosmomc and the relevant posts here, but when I run cosmomc I get the message:

> ./cosmomc params.ini
Random seeds: 15052, 14868 rand_inst: 0
Computing tensors: F
Doing CMB lensing: F
lmax = 2100
Number of C_ls = 3
parameter ranges not found, param 14 - omegal

======
Here are the modification to cosmomc
======

In settings.f90

integer, parameter :: num_hard = 8

In cmbtypes.f90

Type CMBParams

real w_a

In params_CMB.f90

subroutine SetForH
CMB%w_a = Params(8)

subroutine CMBParamsToParams
Params(8) = CMB%w_a

In CMB_Cls_simple.f90

subroutine CMBToCAMB
P%w_a = CMB%w_a

—–
In CAMB

In modules.f90

type CAMBparams

real(dl) :: w_a

In inidrivers.f90

w_a = Ini_Read_Double('w_a')

—–

In params.ini

#parameter start center, min, max, start width, st. dev. estimate
param[omegabh2] = 0.0223 0.0223 0.0223 0 0
param[omegadmh2] = 0.150 0.150 0.150 0 0
param[theta] = 1.04 1.04 1.04 0 0

param[tau] = 0 0 0 0 0

param[omegak] = 0 0 0 0 0
param[fnu] = 0 0 0 0 0
param[w] = −1 −1.5 −0.5 0.01 0.01
param[w_a] = 0 −0.5 0.5 0.01 0.01


param[ns] = 0.95 0.95 0.95 0 0

param[nt] = 0 0 0 0 0
param[nrun] = 0 0 0 0 0

#log[10¹⁰ A_s]
param[logA] = 3 3 3 0 0
param[r] = 0 0 0 0 0


—–

I am only running cosmomc with supernovae data. Any help is appreciated.

Thanks,

Starla

[Lukas Hollenstein]

Hi Starla,

I think you have to add a line in "params_CMB.paramnames" after the line describing "w". Something like

[Starla Talbot]

Hello Lukas,

That indeed was the silly error. Thanks.

Cheers,

Starla.