CosmoCoffee Forum Index CosmoCoffee

 
 FAQFAQ   SearchSearch  MemberlistSmartFeed   MemberlistMemberlist    RegisterRegister 
   ProfileProfile   Log inLog in 
Arxiv New Filter | Bookmarks & clubs | Arxiv ref/author:

CAMB: cmbmain.f90: fixq and density perturbation evolution
 
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
Dahee Lee



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

PostPosted: April 16 2018  Reply with quote

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
Back to top
View user's profile [ Hidden ]
Antony Lewis



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

PostPosted: April 16 2018  Reply with quote

I suggest you use the python wrapper, examples in

http://camb.readthedocs.io/en/latest/CAMBdemo.html
Back to top
View user's profile [ Hidden ] Visit poster's website
Dahee Lee



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

PostPosted: April 17 2018  Reply with quote

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.
Back to top
View user's profile [ Hidden ]
Antony Lewis



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

PostPosted: April 18 2018  Reply with quote

q=k in a flat model.
Back to top
View user's profile [ Hidden ] Visit poster's website
Dahee Lee



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

PostPosted: April 18 2018  Reply with quote

Thank you. I got it.
That question was my misunderstanding.
Back to top
View user's profile [ Hidden ]
Display posts from previous:   
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software All times are GMT + 5 Hours
Page 1 of 1

 
Jump to:  
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


Powered by phpBB © 2001, 2005 phpBB Group. Sponsored by WordWeb online dictionary and dictionary software.