def test_plot_integrate_ion_recomb_collapse(): """Plot results of integrate_ion_recomb_collapse (no quantitative tests). """ cosmo = cparam.WMAP5_mean(flat=True) dz = 0.5 z = numpy.arange(25., 8. - 1.5 * dz, -1. * dz) T_min = 1e4 #K c_ion = numpy.array([[500.], [40.], [12.]]) # calculate ionized fraction from collapse fraction x_fcol = cr.ionization_from_collapse(z, c_ion, T_min, **cosmo) x_rec = numpy.empty(x_fcol.shape) w_rec = numpy.empty(x_fcol.shape) for i in range(x_fcol.shape[0]): # calculate ionized fraction, including recombinations x_rec[i], w_rec[i], t = cr.integrate_ion_recomb_collapse(z, c_ion[i,0], temp_min = T_min, **cosmo) #linestyle = ['-', ':', '--'] color = ['r', 'g', 'b'] pylab.figure() pylab.subplot(2,1,1) for i in range(len(color)): pylab.plot(z, x_fcol[i], ls='--', color=color[i]) pylab.plot(z, x_rec[i], ls='-', color=color[i]) pylab.plot(z, w_rec[i], ls=':', color=color[i]) pylab.axhline(y=0.75) pylab.yscale('log') pylab.xlim(8,25) pylab.ylim(1e-4, 1) pylab.subplot(2,1,2) for i in range(len(color)): pylab.plot(z, x_fcol[i], ls='--', color=color[i]) pylab.plot(z, x_rec[i], ls='-', color=color[i]) pylab.plot(z, w_rec[i], ls=':', color=color[i]) pylab.axhline(y=0.75) pylab.xlim(8,25) pylab.ylim(0, 1)
def test_plot_FZH(): """Plot figure 3 from FZH (2004ApJ...613....1F) (no quantitative tests). """ cosmo = {} cosmo['omega_M_0'] = 0.3 cosmo['omega_k_0'] = 0.0 cosmo['omega_lambda_0'] = 0.7 cosmo['h'] = 0.7 cosmo['n'] = 1.0 cosmo['sigma_8'] = 0.9 cosmo['omega_b_0'] = 0.0462 # not specified in paper cosmo['omega_n_0'] = 0.0 # not specified in paper cosmo['N_nu'] = 0 # not specified in paper cosmo['Y_He'] = 0.24 # not specified in paper dz = 0.1 z = numpy.arange(25., 8. - 1.5 * dz, -1. * dz) T_min = 1e4 #K c_ion = numpy.array([[500.], [40.], [12.]]) # calculate ionized fraction from collapse fraction x_fcol = cr.ionization_from_collapse(z, c_ion, T_min, **cosmo) #linestyle = ['-', ':', '--'] color = ['r', 'g', 'b'] pylab.figure() pylab.subplot(2,1,1) for i in range(len(color)): pylab.plot(z, x_fcol[i], ls='--', color=color[i]) pylab.axhline(y=0.75) pylab.yscale('log') pylab.xlim(8,25) pylab.ylim(1e-4, 1) pylab.title("Compare to figure 3 from FZH (2004ApJ...613....1F)") pylab.subplot(2,1,2) for i in range(len(color)): pylab.plot(z, x_fcol[i], ls='--', color=color[i]) pylab.axhline(y=0.75) pylab.xlim(8,25) pylab.ylim(0, 1)
def test_plot_FZH(): """Plot figure 3 from FZH (2004ApJ...613....1F) (no quantitative tests). """ cosmo = {} cosmo['omega_M_0'] = 0.3 cosmo['omega_k_0'] = 0.0 cosmo['omega_lambda_0'] = 0.7 cosmo['h'] = 0.7 cosmo['n'] = 1.0 cosmo['sigma_8'] = 0.9 cosmo['omega_b_0'] = 0.0462 # not specified in paper cosmo['omega_n_0'] = 0.0 # not specified in paper cosmo['N_nu'] = 0 # not specified in paper cosmo['Y_He'] = 0.24 # not specified in paper cosmo['baryonic_effects'] = False dz = 0.1 z = numpy.arange(25., 8. - 1.5 * dz, -1. * dz) T_min = 1e4 #K c_ion = numpy.array([[500.], [40.], [12.]]) # calculate ionized fraction from collapse fraction x_fcol = cr.ionization_from_collapse(z, c_ion, T_min, **cosmo) #linestyle = ['-', ':', '--'] color = ['r', 'g', 'b'] pylab.figure() pylab.subplot(2,1,1) for i in range(len(color)): pylab.plot(z, x_fcol[i], ls='--', color=color[i]) pylab.axhline(y=0.75) pylab.yscale('log') pylab.xlim(8,25) pylab.ylim(1e-4, 1) pylab.title("Compare to figure 3 from FZH (2004ApJ...613....1F)") pylab.subplot(2,1,2) for i in range(len(color)): pylab.plot(z, x_fcol[i], ls='--', color=color[i]) pylab.axhline(y=0.75) pylab.xlim(8,25) pylab.ylim(0, 1)
def test_plot_FZH(): """Plot figure 3 from FZH (2004ApJ...613....1F) (no quantitative tests). """ cosmo = {} cosmo["omega_M_0"] = 0.3 cosmo["omega_k_0"] = 0.0 cosmo["omega_lambda_0"] = 0.7 cosmo["h"] = 0.7 cosmo["n"] = 1.0 cosmo["sigma_8"] = 0.9 cosmo["omega_b_0"] = 0.0462 # not specified in paper cosmo["omega_n_0"] = 0.0 # not specified in paper cosmo["N_nu"] = 0 # not specified in paper cosmo["Y_He"] = 0.24 # not specified in paper cosmo["baryonic_effects"] = False dz = 0.1 z = numpy.arange(25.0, 8.0 - 1.5 * dz, -1.0 * dz) T_min = 1e4 # K c_ion = numpy.array([[500.0], [40.0], [12.0]]) # calculate ionized fraction from collapse fraction x_fcol = cr.ionization_from_collapse(z, c_ion, T_min, **cosmo) # linestyle = ['-', ':', '--'] color = ["r", "g", "b"] pylab.figure() pylab.subplot(2, 1, 1) for i in range(len(color)): pylab.plot(z, x_fcol[i], ls="--", color=color[i]) pylab.axhline(y=0.75) pylab.yscale("log") pylab.xlim(8, 25) pylab.ylim(1e-4, 1) pylab.title("Compare to figure 3 from FZH (2004ApJ...613....1F)") pylab.subplot(2, 1, 2) for i in range(len(color)): pylab.plot(z, x_fcol[i], ls="--", color=color[i]) pylab.axhline(y=0.75) pylab.xlim(8, 25) pylab.ylim(0, 1)