Page 1 of 1

CAMB:Get overdensities as a function of scale factor

Posted: October 25 2016
by Tomohiro Nozawa
Hello all,

I would like to make a graph of "overdensities of CDM, baryon, photon and neutrino vs. scale factor " for a particular wavenumber e.g. k = 1, 0.1, 0.01 /Mpc.
I have the overdensities as a function of wave number k in transfer.dat output file, but how could I get those as a function of scale factor?

I really appreciate any help you can provide.

Tomohiro

Re: CAMB:Get overdensities as a function of scale factor

Posted: October 25 2016
by Antony Lewis
There are examples of how to make time evolution plots in the sample python Jupyter notebook:

http://camb.readthedocs.io/en/latest/CAMBdemo.html

CAMB:Get overdensities as a function of scale factor

Posted: October 28 2016
by Erick Almaraz
Dear Tomohiro.

In cmbmain.f90 there is a subroutine called CalcScalarSources. There is an option to create a txt file for plotting the evolution of some quantity (eg, the overdensities of CDM, baryons, etc). You need to specify the mode you want to study in

Code: Select all

    real(dl) :: fixq = 0._dl !Debug output of one q
This sentence is located at the top of cmbmain.f90.

Regards