CosmoCoffee

Hi,

I am new to CAMB. I never did any serious calculations using Fortran as I use Mathematica. This question may have already been asked, but I cannot find it. Please pardon me. I will be grateful if you help me out.

I have interacting dark matter dark energy perturbation equations in conformal Newtonian gauge. The EoS parameter for dark energy w_{de} is not a constant but depends on the redshift, z. Other background quantities are also evolving with redshift.

My questions are:
(a) Do I have to change my equations into the synchronous gauge?
(b) Do I need to change the GuageInterface module of equations.f90 for my equations?
(c) How do I incorporate the varying background in the code?
(d) Where do I put the initial conditions?
(e) How to do these if I want to use pyCAMB?

Thank you a lot in advance.


Srijita

a,b: yes. c and d are also implemented in equations.f90.

Once you have modified the fortran, using the python version should be straightforward, you'll just need to add any new parameters to CAMBparams in both fortran and the python wrapper.