Adding a new model parameter in CosmoMC
Posted: December 19 2019
Hi All,
I've unsuccessfully tried to add a new model parameter, alpha, to the latest CosmoMC code (sumer 2019) version. I did follow the instructions appearing in one of the readme's FAQ and was able to add new model parameters in an older, November 2016, version. However, all my efforts to modify the new version so as to accommodate a new model parameter have failed. Indeed, I have noticed significant structural changes between the old a new versions, so I guess this might be the reason.
Could someone instruct me as to the additional changes that should be applied to the new version to allow inclusion of a new cosmological parameter ?
Changes carried out in my most recent attempts to modify Cosmomc 2019:
(1) I defined the new parameter in CosmologyTypes.f90
(2) In CosmologyParameterizations.f90 I replaced:
call this%SetTheoryParameterNumbers(16,last_power_index)
with
call this%SetTheoryParameterNumbers(17,last_power_index)
in addition to defining (in subroutine SetForH):
CMB%alpha = Params(17)
(3) I defined P%alpha = CMB%alpha in subroutine CAMBCalc_CMBToCAMB in Calculator_CAMB.f90
(4) I modified the results.f90 file to effect my cosmological model changes.
(5) I've added the new parameter alpha with appropriate range to the parameter .ini file in batch3.
Thanks,
Meir
I've unsuccessfully tried to add a new model parameter, alpha, to the latest CosmoMC code (sumer 2019) version. I did follow the instructions appearing in one of the readme's FAQ and was able to add new model parameters in an older, November 2016, version. However, all my efforts to modify the new version so as to accommodate a new model parameter have failed. Indeed, I have noticed significant structural changes between the old a new versions, so I guess this might be the reason.
Could someone instruct me as to the additional changes that should be applied to the new version to allow inclusion of a new cosmological parameter ?
Changes carried out in my most recent attempts to modify Cosmomc 2019:
(1) I defined the new parameter in CosmologyTypes.f90
(2) In CosmologyParameterizations.f90 I replaced:
call this%SetTheoryParameterNumbers(16,last_power_index)
with
call this%SetTheoryParameterNumbers(17,last_power_index)
in addition to defining (in subroutine SetForH):
CMB%alpha = Params(17)
(3) I defined P%alpha = CMB%alpha in subroutine CAMBCalc_CMBToCAMB in Calculator_CAMB.f90
(4) I modified the results.f90 file to effect my cosmological model changes.
(5) I've added the new parameter alpha with appropriate range to the parameter .ini file in batch3.
Thanks,
Meir