CAMB: Reionization set from .true. to .false. after CAMB_GetResults is called

Use of Cobaya. camb, CLASS, cosmomc, compilers, etc.
Post Reply
Anais Rassat
Posts: 24
Joined: October 20 2004
Affiliation: LASTRO (EPFL)
Contact:

CAMB: Reionization set from .true. to .false. after CAMB_Get

Post by Anais Rassat » October 23 2008

Hi,

I'm having a small problem setting the Reionization in CAMB.
In my params.ini file I set:

---
#Reionization (assumed sharp), ignored unless reionization = T
reionization = T
---

Then I've modified the Inidriver.F90 to output the value of P%Reionization such that (and these are the only modifications to inidriver.F90 that I've done):

---
write(*,*) 'Reionization', P%Reionization
write(*,*) 'Starting CAMB_GetResults'
call CAMB_GetResults(P)
write(*,*) 'Finished CAMB_GetResults'
write(*,*) 'Reionization', P%Reionization
---

When I run ./inidriver params.ini the screen output is:

---
Reionization T
Starting CAMB_GetResults
Finished CAMB_GetResults
Reionization F
---

Is this normal? Does this mean I need to reset P%Reionization = .true. everytime I run CAMB_GetResults?

Cheers, Anais.

Anais Rassat
Posts: 24
Joined: October 20 2004
Affiliation: LASTRO (EPFL)
Contact:

CAMB: Reionization set from .true. to .false. after CAMB_Get

Post by Anais Rassat » October 23 2008

In the above post I am using CAMB version from February 2008

Post Reply