Hello everybody,
I have uncommented in cmbmain.f90,in particular in the CalcScalarSources subroutine,the part "Example code for plotting out variable evolution" but when I run the program the output is "Not a Number" for every y except the first one(y(1));could anyone tell me where I am wrong?
Thanks
Daniela
CAMB:plotting out the variable evolution
-
- Posts: 2
- Joined: September 21 2006
- Affiliation: bologna university
-
- Posts: 1942
- Joined: September 23 2004
- Affiliation: University of Sussex
- Contact:
Re: CAMB:plotting out the variable evolution
Don't uncomment the first line with
- so uncommented code would read
(also setting desired k a few lines up of course)
Code: Select all
call GaugeInterface_EvolveScal(EV,tau,y,CP%tau0,tol1,ind,c,w)
Code: Select all
tol1=tol/exp(AccuracyBoost-1)
do j=1,6000
tauend = taustart * exp(j/6000._dl*log(CP%tau0/taustart))
call GaugeInterface_EvolveScal(EV,tau,y,tauend,tol1,ind,c,w)
write (*,'(8E15.5)') y(1:7), grhoc*y(1)/(grhog)
end do
stop