CosmoCoffee Forum Index CosmoCoffee

 
 FAQFAQ   SearchSearch  MemberlistSmartFeed   MemberlistMemberlist    RegisterRegister 
   ProfileProfile   Log inLog in 
Arxiv New Filter | Bookmarks & clubs | Arxiv ref/author:

Some troubles about adding new parameters in CosmoMC
 
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
Yutong Wang



Joined: 06 May 2014
Posts: 5
Affiliation: university of Chinese Academy of Sciences

PostPosted: October 01 2016  Reply with quote

Hi everyone!

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:
Ps(k)=P0(k)*(1-exp(-(k/k_c)lc))),

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
Code:
derived(ix) = CMB%InitPower(ns_index) +CMB%InitPower(nrun_index)*lograt +&
            CMB%InitPower(nrunrun_index)*lograt**2/2

as
Code:
derived(ix) = CMB%InitPower(ns_index) +&
            CMB%InitPower(nrunrun_index)*lograt**2/2

in cosmomc/batch2/params_CMB_defaults.ini, I make the modification
Code:
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+
Back to top
View user's profile  
Display posts from previous:   
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software All times are GMT + 5 Hours
Page 1 of 1

 
Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum


Powered by phpBB © 2001, 2005 phpBB Group. Sponsored by WordWeb online dictionary and dictionary software.