CAMB:Get overdensities as a function of scale factor

Use of Healpix, camb, CLASS, cosmomc, compilers, etc.
Post Reply
Tomohiro Nozawa
Posts: 1
Joined: September 09 2016
Affiliation: graduate student

CAMB:Get overdensities as a function of scale factor

Post by Tomohiro Nozawa » October 25 2016

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

Antony Lewis
Posts: 1369
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

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

Post by Antony Lewis » October 25 2016

There are examples of how to make time evolution plots in the sample python Jupyter notebook:

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

Erick Almaraz
Posts: 25
Joined: June 10 2012
Affiliation: University of Mexico

CAMB:Get overdensities as a function of scale factor

Post by Erick Almaraz » October 28 2016

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

Post Reply