CosmoCoffee Forum Index CosmoCoffee

 
 FAQFAQ   SearchSearch  MemberlistSmartFeed   MemberlistMemberlist    RegisterRegister 
   ProfileProfile   Log inLog in 
Arxiv New Filter | Bookmarks & clubs | Arxiv ref/author:

clxg for lcdm
 
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software
View previous topic :: View next topic  
Author Message
Omar Ruiz Macias



Joined: 15 Feb 2014
Posts: 10
Affiliation: Cinvestav

PostPosted: May 30 2016  Reply with quote

Hi everyone,

I have a question about clxg (photons) perturbations. I was comparing the CAMB outcome using w=−1 (labeled 'clxg_lcdm ') and w=−0.99999999999 (labeled 'clxg_lcdm_false'). Of course, I was expecting the spectra to be indistinguishable and they are. However, checking the clxg as function of the scale factor (a), I have found the difference shown in the figure.

clxg_a_lcdm.png

What may be the reason for this difference? And, Why it takes place at this particular value of a?
Back to top
View user's profile  
Antony Lewis



Joined: 23 Sep 2004
Posts: 1203
Affiliation: University of Sussex

PostPosted: May 31 2016  Reply with quote

Is this using equations.f90 or equations_ppf.f90?
Back to top
View user's profile [ Hidden ] Visit poster's website
Omar Ruiz Macias



Joined: 15 Feb 2014
Posts: 10
Affiliation: Cinvestav

PostPosted: June 13 2016  Reply with quote

It's using 'equations.f90' and 'fixq = 1._dl'.
Back to top
View user's profile  
Antony Lewis



Joined: 23 Sep 2004
Posts: 1203
Affiliation: University of Sussex

PostPosted: June 13 2016  Reply with quote

I think you are probably plotting different things here. With do_late_rad_truncation = T, the photon perturbations are not calculated at late time, so the array element you are outputting will be something else (dark energy perturbation or something else, depending on whether w=−1 and hence whether there are any dark energy perturbations).

The best way to plot parameter evolution is to use the python wrapper
Code:

from matplotlib import pyplot as plt
import camb
import numpy as np

pars = camb.set_params(H0=67.5, ombh2=0.022, omch2=0.122, As=2e-9, ns=0.95)
data= camb.get_background(pars)

a=10**(np.linspace(-8, 0, 500))
z = 1/a -1
eta = data.conformal_time(z)
ev = data.get_time_evolution(1., eta, ['delta_photon'])
plt.loglog(a, np.abs(ev))
Back to top
View user's profile [ Hidden ] Visit poster's website
Omar Ruiz Macias



Joined: 15 Feb 2014
Posts: 10
Affiliation: Cinvestav

PostPosted: June 15 2016  Reply with quote

thank you, Antony. Now it calculates the photon perturbations at late time.
Back to top
View user's profile  
Display posts from previous:   
Post new topic   Reply to topic    CosmoCoffee Forum Index -> Computers and software All times are GMT + 5 Hours
Page 1 of 1

 
Jump to:  
You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot vote in polls in this forum


Powered by phpBB © 2001, 2005 phpBB Group. Sponsored by WordWeb online dictionary and dictionary software.