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).
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: Select all
real(dl) :: fixq = 0._dl !Debug output of one q
Code: Select all
!!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