cosmomc: adding an extra parameter
-
- Posts: 1944
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
-
- Posts: 14
- Joined: May 29 2015
- Affiliation: BITS-Pilani, India
cosmomc: adding an extra parameter
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
I added
(2) In CosmologyParameterizations.f90, I modified the code line
by
(3) In Calculator_CAMB.f90, below the code line
I added
In the same file, below the code line
I added
(4) In params_CMB.paramnames, below the code line
I added
(5) In params_CMB_defaults.ini, I added
After the above changes, when I run the code, it stops with the following error:
Please reply in detail so that I can run the code correctly with the new parameter.
(1) In CosmologyTypes.f90, below the code line
Code: Select all
real(mcp) w, wa
Code: Select all
real(mcp) delta_coupled
Code: Select all
call this%SetTheoryParameterNumbers(16,last_power_index)
Code: Select all
call this%SetTheoryParameterNumbers(17,last_power_index)
Code: Select all
real(dl) neff_massive_standard
Code: Select all
real(dl) delta_coupled
Code: Select all
wa_ppf = CMB%wa
Code: Select all
delta_coupled = CMB%delta_coupled
Code: Select all
wa w_a #w_a variation
Code: Select all
delta_coupled delta_coupled
Code: Select all
param[delta_coupled] = 0 -1 1 0.01 0.01
Code: Select all
Fast divided into 1 blocks
25 parameters (11 slow ( 0 semi-slow), 14 fast ( 0 semi-fast))
starting Monte-Carlo
Compile with CosmoRec to use fdm
MpiStop: 0
-
- Posts: 1944
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: cosmomc: adding an extra parameter
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.
-
- Posts: 14
- Joined: May 29 2015
- Affiliation: BITS-Pilani, India
cosmomc: adding an extra parameter
It is working now. Thanks.
-
- Posts: 14
- Joined: May 29 2015
- Affiliation: BITS-Pilani, India
cosmomc: adding an extra parameter
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?
-
- Posts: 1944
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: cosmomc: adding an extra parameter
I don't think it's interfaced to CAMB.
-
- Posts: 19
- Joined: March 03 2017
- Affiliation: IIT KANPUR
- Contact:
cosmomc: adding an extra parameter
I am also getting the same error .Please explain how to resolve it?