Adding Hdot to the background and perturbations in CAMB

[Maryam A. Abchouyeh]

Hi,

I am trying to modify CAMB for a cosmological model that includes Hdot in the first Friedman equation, with an additional constant coefficient.
As I understand, first of all I should change dtauda in equation.f90, but to modify the code properly, at least results.f90 and inidriver.f90 should also be modified.
Would you please let me know which part of these three files should be changed to create a reliable result?
And is there any other file that needs a modification accordingly?

Thanks

[Antony Lewis]

The background equations are also used in derivs() for the perturbations. However, if the background Einstein equation changes, usually the perturbed one does as well, so you'd need to modify the perturbation equations consistently anyway. Depending on your particular model, and what outputs you want to calculate, you'll have to go through the code checking relevant things that have changed.

[Maryam A. Abchouyeh]

The background equations are also used in derivs() for the perturbations. However, if the background Einstein equation changes, usually the perturbed one does as well, so you'd need to modify the perturbation equations consistently anyway. Depending on your particular model, and what outputs you want to calculate, you'll have to go through the code checking relevant things that have changed.

Thanks.
Equations for perturbations are only in drives(), or I should also look for them in other files?

[Antony Lewis]

derivs for evolution, there's also initial conditions in equations.f90. And the NonLinear model for late-time lensing etc.