In [1]:
%pylab inline
import numpy as np
import matplotlib.pyplot as plt
import analytic_p1d_PD2013 as aP1D
In [2]:
# s/km units
k=np.logspace(-4,-1,1000)
pk_z2=aP1D.P1D_z_kms_PD2013(2.0,k)
pk_z3=aP1D.P1D_z_kms_PD2013(3.0,k)
pk_z4=aP1D.P1D_z_kms_PD2013(4.0,k)
In [3]:
# something related to fonts in the plot
plt.rc('text', usetex=True)
plt.rcParams['text.latex.preamble'] = [
r'\usepackage{siunitx}', # i need upright \micro symbols, but you need
r'\sisetup{detect-all}', # this to force siunitx to actually use your fonts
r'\usepackage{helvet}', # set the normal font here
r'\usepackage{sansmath}', # load up the sansmath so that math -> helvet
r'\sansmath' # <- tricky! -- gotta actually tell tex to use!
]
plt.rc('xtick', labelsize=15)
plt.rc('ytick', labelsize=15)
In [4]:
plt.xlim(1e-4,0.1)
plt.xlabel(r'k [s $\rm{km}^{-1}$]',fontsize=15)
plt.ylabel(r'k P(k) / $\pi$',fontsize=15)
plt.title(r'1D flux power',fontsize=15)
plt.loglog(k,k*pk_z2/np.pi,label=r'z=2.0')
plt.loglog(k,k*pk_z3/np.pi,label=r'z=3.0')
plt.loglog(k,k*pk_z4/np.pi,label=r'z=4.0')
plt.legend(fontsize=15)
plt.show()
In [5]:
plt.xlim(1e-4,0.1)
plt.xlabel(r'k [s $\rm{km}^{-1}$]',fontsize=15)
plt.ylabel(r'k P(k) / $\pi$',fontsize=15)
plt.title(r'1D flux power',fontsize=15)
plt.loglog(k,pk_z2,label=r'z=2.0')
plt.loglog(k,pk_z3,label=r'z=3.0')
plt.loglog(k,pk_z4,label=r'z=4.0')
plt.legend(fontsize=15)
plt.show()