def thetamc(**cosmo):
omegab = cosmo["omega_b_0"]
omegac = cosmo["omega_M_0"] - cosmo["omega_b_0"] - cosmo["omega_n_0"]
omegav = cosmo["omega_lambda_0"]
omegan = cosmo["omega_n_0"]
H0 = cosmo["h"] * 100
Num_Nu_massless = cosmo["Num_Nu_massless"]
Num_Nu_massive = cosmo["Num_Nu_massive"]
omegak = 1.0 - omegav - omegac - omegab - omegan
w_lam = cosmo["w"]
thetamc = pycamb.cosmomctheta(
omegab=omegab,
omegac=omegac,
omegav=omegav,
omegan=omegan,
H0=H0,
Num_Nu_massless=Num_Nu_massless,
Num_Nu_massive=Num_Nu_massive,
omegak=omegak,
w_lam=w_lam,
)
return thetamc
