CAMB: Reionization set from .true. to .false. after CAMB_Get
Posted: 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.
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.