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
CAMB:Get overdensities as a function of scale factor
-
- Posts: 1
- Joined: September 09 2016
- Affiliation: graduate student
-
- Posts: 1951
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: CAMB:Get overdensities as a function of scale factor
There are examples of how to make time evolution plots in the sample python Jupyter notebook:
http://camb.readthedocs.io/en/latest/CAMBdemo.html
http://camb.readthedocs.io/en/latest/CAMBdemo.html
-
- Posts: 29
- Joined: June 10 2012
- Affiliation: University of Mexico
CAMB:Get overdensities as a function of scale factor
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
This sentence is located at the top of cmbmain.f90.
Regards
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
Regards