PPF tweak for better consistency in some extreme models

cmbant · Jun 14, 2024 · 9a4c01a · 9a4c01a
1 parent 3ad6173
commit 9a4c01a
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/fortran/DarkEnergyPPF.f90 b/fortran/DarkEnergyPPF.f90
@@ -129,7 +129,9 @@ subroutine TDarkEnergyPPF_PerturbedStressEnergy(this, dgrhoe, dgqe, &
     S_Gamma = grhov_t * (1 + w) * (vT + sigma) * k / adotoa / 2._dl / k2
     ckH = this%c_Gamma_ppf * k / adotoa
 
-    if (ckH * ckH .gt. 1e3) then !3.d1) then ! ckH^2 > 30 ?????????
+    if (ckH * ckH > 1000) then
+        ! Was ckH^2 > 30 originally, but this is better behaved (closer to fluid)
+        ! for some extreme models (thanks Yanhui Yang, Simeon Bird 2024)
         Gamma = 0
         Gammadot = 0.d0
     else