Model universes

import numpy as np
import simple_python_package_template

Einsten-de Sitter universe

param_EdS = simple_python_package_template.param()
param_EdS.cosmo.Om  = 1.
param_EdS.cosmo.Or  = 0
param_EdS.cosmo.Ode = 0
param_EdS.cosmo.Ok  = 0
param_EdS.cosmo.h   = 0.68
param_EdS.code.Nz   = 100

print('Cosmological parameters')
print(param_EdS.cosmo.__dict__)
print('Code parameters')
print(param_EdS.code.__dict__)

Cosmological parameters
{'Om': 1.0, 'Or': 0, 'Ok': 0, 'Ode': 0, 'h': 0.68}
Code parameters
{'zmin': 0.01, 'zmax': 9.0, 'Nz': 100, 'verbose': True}

D_EdS = simple_python_package_template.CosmoDistances(param_EdS)

de Sitter universe

param_dS = simple_python_package_template.param()
param_dS.cosmo.Om  = 0.
param_dS.cosmo.Or  = 0
param_dS.cosmo.Ode = 1.
param_dS.cosmo.Ok  = 0
param_dS.cosmo.h   = 0.68
param_dS.code.Nz   = 100

print('Cosmological parameters')
print(param_dS.cosmo.__dict__)
print('Code parameters')
print(param_dS.code.__dict__)

Cosmological parameters
{'Om': 0.0, 'Or': 0, 'Ok': 0, 'Ode': 1.0, 'h': 0.68}
Code parameters
{'zmin': 0.01, 'zmax': 9.0, 'Nz': 100, 'verbose': True}

D_dS = simple_python_package_template.CosmoDistances(param_dS)

Planck universe

param = simple_python_package_template.param()
param.cosmo.Om  = 0.31
param.cosmo.Or  = 0.0
param.cosmo.Ode = 0.69
param.cosmo.Ok  = 0.0
param.cosmo.h   = 0.68
param.code.Nz   = 100

print('Cosmological parameters')
print(param.cosmo.__dict__)
print('Code parameters')
print(param.code.__dict__)

Cosmological parameters
{'Om': 0.31, 'Or': 0.0, 'Ok': 0.0, 'Ode': 0.69, 'h': 0.68}
Code parameters
{'zmin': 0.01, 'zmax': 9.0, 'Nz': 100, 'verbose': True}

D = simple_python_package_template.CosmoDistances(param)

Plot models

import matplotlib.pyplot as plt

zs = np.linspace(param.code.zmin,param.code.zmax,param.code.Nz)
fig, ax = plt.subplots(1,1,figsize=(5,4))
ax.loglog(zs, D_EdS.comoving_dist(zs), ls='--', label='EdS')
ax.loglog(zs, D_dS.comoving_dist(zs), ls='-.', label='dS')
ax.loglog(zs, D.comoving_dist(zs), ls='-', label='Planck')
ax.set_xlabel('Redshift')
ax.set_ylabel('Comoving distance')
ax.axis([0.01,10,0])
plt.show()