I have started using cosmomc very recently, so I do not know exactly how to add a parameter in cosmomc. Reading some posts I have manged (seems to be so) to add a new parameter. My parameter is included in the
initial power spectrum ("power_tilt.f90").
Here are my steps: Please check and reply whether I have done correct changes or not.
1. in /camb in "power_tilt.f90" I defined my model parameter
by: real(dl) p(nnmax)
then in the subroutine SetDefPowerParams(AP)
AP%p = 1.0
then in subroutine InitialPower_ReadParams(InitPower, Ini, WantTensors)
InitPower%p = Ini_Read_Double_Array_File(Ini,'parameter_1',i,1._dl)
And in the "params.ini" file I have added: parameter_1(1)=0.8
Thats all I have done in the CAMB (did not added it to inidriver.f90 or modules.f90).
2. in /source I have first defined a new parameter in "cmbtypes.f90"
by: real pp
in "settings.f90" I have changed: num_hard=7 to num_hard=8
in "params_CMB.f90" in subroutine SetCAMBInitPower(P,CMB,in) I have added P%InitPower%p(in) = CMB%pp
in subroutine SetforH I have added CMB%pp = Params(8)
in subroutine CMBParamsToParams I have added Params(8)= CMB%pp
Did not change "cmb_Cls_Simple.f90"
I also modified "params_CMB.paramsnames" and added param[pp]= 0.6 0.5 0.7 0.01 0.01 in params.ini file.
I did not change anything else. Have I done the right modifications??
Please help.
Sorry for poor english.