I'm working on a project where I need to modify the parameterisation w(a) = w0 + wa * (1-a) from the equations_PPF.f90 file used in CAMB and CosmoMC.

I can tell that I definitely need to modify line 98:

Code: Select all

`w_de=w_lam+wa_ppf*(1._dl-a)`

Code: Select all

`grho_de=grhov*a**(1._dl-3.*w_lam-3.*wa_ppf)*exp(-3.*wa_ppf*(1._dl-a))`

Code: Select all

`f(a) = a**(-3.*(1+w_lam+wa_ppf)) * exp(-3*wa_ppf*(1.-a))`

So, in summary - what role does line 143 of equations_PPF.f90 play in CAMB, and if I want to modify the w(a) = w0 + wa * (1-a) parameterisation do I also need to modify line 143 and not just line 98?

Thank you!

- Per