CosmoCoffee

 FAQ   Search  SmartFeed   Memberlist    Register Profile   Log in Arxiv New Filter | Bookmarks & clubs | Arxiv ref/author:

 CAMB: cmbmain.f90: fixq and density perturbation evolution
Author Message
Dahee Lee

Joined: 16 Apr 2018
Posts: 8
Affiliation: Kyungpook National University

Posted: April 16 2018

Hello.
I'm beginner of CAMB.

I would like to plot density perturbation evolution using cmbmain.f90.
In cmbmain.f90, I use fixq = 0 for computing power spectrum(scalCls.dat).

 Quote: q is used instead of k in parts of the code (for non-flat models). You can write out whatever you like.

At this sentence, I understood that fixq is k for non-flat model.
So, how can I use it for flat model?
I'm so confused.

For plotting density perturbation evolution, I changed value of fixq(like 1, 0.1, 0.01)
but I think that is non-flat because of changing fixq.

I want to calculate density perturbation evolution and scalCls.dat at the same time.
How can I do for this?

Any help would be appreciated.

 Code: real(dl) :: fixq = 0._dl !Debug output of one q

 Code: !!Example code for plotting out variable evolution if (fixq/=0._dl) then tol1=tol/exp(AccuracyBoost-1) call CreateTxtFile('evolve.txt',1) do j=1,1000 tauend = taustart+(j-1)*(CP%tau0-taustart)/1000 call GaugeInterface_EvolveScal(EV,tau,y,tauend,tol1,ind,c,w) yprime = 0 call derivs(EV,EV%ScalEqsToPropagate,tau,y,yprime) adotoa = 1/(y(1)*dtauda(y(1))) ddelta= (yprime(3)*grhoc+yprime(4)*grhob)/(grhob+grhoc) delta=(grhoc*y(3)+grhob*y(4))/(grhob+grhoc) growth= ddelta/delta/adotoa write (1,'(7E15.5)') tau, delta, growth, y(3), y(4), y(EV%g_ix), y(1) end do close(1) stop end if
Antony Lewis

Joined: 23 Sep 2004
Posts: 1350
Affiliation: University of Sussex

 Posted: April 16 2018 I suggest you use the python wrapper, examples in http://camb.readthedocs.io/en/latest/CAMBdemo.html
Dahee Lee

Joined: 16 Apr 2018
Posts: 8
Affiliation: Kyungpook National University

 Posted: April 17 2018 Thank you very much for your answer. But I still don't know about fixq for flat model. Is it right that changing fixq(like 1, 0.1, 0.01) is non-flat? Any help would be appreciated.
Antony Lewis

Joined: 23 Sep 2004
Posts: 1350
Affiliation: University of Sussex

 Posted: April 18 2018 q=k in a flat model.
Dahee Lee

Joined: 16 Apr 2018
Posts: 8
Affiliation: Kyungpook National University

 Posted: April 18 2018 Thank you. I got it. That question was my misunderstanding.
 Display posts from previous: All Posts1 Day7 Days2 Weeks1 Month3 Months6 Months1 Year Oldest FirstNewest First
 All times are GMT + 5 Hours Page 1 of 1

 Jump to: Select a forum Arxiv paper discussion----------------arXiv papers Topic discussion----------------Early UniverseCosmological ModelCosmological Observations  Projects and Resources----------------Computers and softwareTeaching, Papers and PresentationsResearch projectsiCosmo Coming up----------------Job vacanciesConferences and meetings Management----------------CosmoCoffee
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum