diff --git a/src/mgcamb.f90 b/src/mgcamb.f90 index b42fd5e..f84dad7 100644 --- a/src/mgcamb.f90 +++ b/src/mgcamb.f90 @@ -799,9 +799,7 @@ function MGCAMB_Gammadot( a, mg_par_cache, mg_cache ) if ( muSigma_par == 1 ) then ! DES parametrization ! changing the following - !omegaDE_t = mg_cache%grhov_t / a**2 / 3._dl / mg_par_cache%h0_Mpc**2 - !omegaDEdot = - 3._dl * mg_cache%adotoa * (mg_cache%grhov_t + mg_cache%gpresv_t) & - ! & / a**2 / 3._dl / mg_par_cache%h0_Mpc**2 + omegaDE_t = mg_cache%grhov_t / 3._dl / mg_cache%adotoa**2 omegaDEdot =-(mg_cache%grhov_t+3._dl*mg_cache%gpresv_t)/3._dl/mg_cache%adotoa & & - 2._dl*mg_cache%Hdot/3._dl/mg_cache%adotoa**3*mg_cache%grhov_t sigma_t = 1._dl + sigma0 * omegaDE_t / mg_par_cache%omegav