fixed bug in mgcamb.f90 for mu-sigma

sfu-cosmo · Jan 15, 2020 · a156a07 · a156a07
1 parent 2397037
commit a156a07
Showing 1 changed file with 1 addition and 3 deletions.
diff --git 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