CAMB:plotting out the variable evolution

Use of Cobaya. camb, CLASS, cosmomc, compilers, etc.
Post Reply
daniela paoletti
Posts: 2
Joined: September 21 2006
Affiliation: bologna university

CAMB:plotting out the variable evolution

Post by daniela paoletti » September 21 2006

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

Antony Lewis
Posts: 1941
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

Re: CAMB:plotting out the variable evolution

Post by Antony Lewis » September 21 2006

Don't uncomment the first line with

Code: Select all

call GaugeInterface_EvolveScal(EV,tau,y,CP%tau0,tol1,ind,c,w)
- so uncommented code would read

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
(also setting desired k a few lines up of course)

Post Reply