Joined: 06 May 2014
|Posted: October 01 2016
I am trying to change cosmomc for adding two extra, new parameters.
In new model, the scalar power spectrum multiplied by an exponential cut-off:
I wanna to get the bestfit value of "kc" and "lc", following the procedure in CosmoMC readme I added the "kc" successfully.
for "lc", I use an easier and faster way, to modify the parameter "nrun" which not be use in my model and use it to instead of "lc".
However, For Planck TT,TE,EE+lowP, the best-fit model has nrun=0.63, and ln(kc/Mpc−1)=−7.96. (the correct bestfit is lc=0.53(nrun=0.53) and ln(kc/Mpc−1)=−7.98) Maybe there is something wrong with the modification about nrun.
I made some modifications for "nrun" in the following files:
in cosmomc/CAMB/power_tilt.f90, I modify the scalar power-law spectrum
as my new model.
in cosmomc/source/CosmologyParameterizations.f90, I modify the
|derived(ix) = CMB%InitPower(ns_index) +CMB%InitPower(nrun_index)*lograt +&
|derived(ix) = CMB%InitPower(ns_index) +&
in cosmomc/batch2/params_CMB_defaults.ini, I make the modification
|param[nrun] = 0.6, 0, 1, 0.01, 0.01 |
and in others files, I didn't make modification about nrun, and anyone can help me to point out the mistake in this process ? It puzzled me for a long time, thanks a lot!
PS:the version of CosmoMC is 2015 July, and ifort is 14+