CosmoCoffee

Posted: **August 18 2016**

It is in the readme:

http://cosmologist.info/cosmomc/readme.html#FAQ

Posted: **August 25 2016**

Thanks. But still it is not working. Following the read me, I did the following to add a new parameter namely "delta_coupled".

(1) In CosmologyTypes.f90, below the code line

Code: Select all

 real&#40;mcp&#41; w, wa

I added

Code: Select all

real&#40;mcp&#41; delta_coupled

(2) In CosmologyParameterizations.f90, I modified the code line

Code: Select all

call this%SetTheoryParameterNumbers&#40;16,last_power_index&#41;

by

Code: Select all

call this%SetTheoryParameterNumbers&#40;17,last_power_index&#41;

(3) In Calculator_CAMB.f90, below the code line

Code: Select all

 real&#40;dl&#41; neff_massive_standard

I added

Code: Select all

real&#40;dl&#41; delta_coupled

In the same file, below the code line

Code: Select all

 wa_ppf = CMB%wa

I added

Code: Select all

delta_coupled = CMB%delta_coupled

(4) In params_CMB.paramnames, below the code line

Code: Select all

wa            w_a       #w_a variation

I added

Code: Select all

delta_coupled delta_coupled

(5) In params_CMB_defaults.ini, I added

Code: Select all

param&#91;delta_coupled&#93; = 0 -1 1 0.01 0.01

After the above changes, when I run the code, it stops with the following error:

Code: Select all

Fast divided into            1  blocks
 25 parameters &#40;11 slow &#40; 0 semi-slow&#41;, 14 fast &#40; 0 semi-fast&#41;&#41;
 starting Monte-Carlo
 Compile with CosmoRec to use fdm
 MpiStop&#58;            0

Please reply in detail so that I can run the code correctly with the new parameter.

Posted: **August 25 2016**

I think you need to get/set the value of your new CMB%delta_coupled parameter in CosmologyParameterizations.f90, updating indexing of other params if necessary.

Posted: **August 25 2016**

It is working now. Thanks.

Posted: **August 26 2016**

I see that the parameter cs2_lam is used in equations.f90, but it is not listed in params_CMB_defaults.ini and params_CMB.paramnames just like w_lam. Is it hard coded or listed somewhere else?

Posted: **August 27 2016**

I don't think it's interfaced to CAMB.

Posted: **April 05 2017**

I am also getting the same error .Please explain how to resolve it?

CosmoCoffee

cosmomc: adding an extra parameter

Re: cosmomc: adding an extra parameter

cosmomc: adding an extra parameter

Re: cosmomc: adding an extra parameter

cosmomc: adding an extra parameter

cosmomc: adding an extra parameter

Re: cosmomc: adding an extra parameter

cosmomc: adding an extra parameter