fix the num_perturb_eqs in DarkEnergy

sfu-cosmo · Aug 15, 2024 · cd0c51c · cd0c51c
1 parent 22f0444
commit cd0c51c
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 6 deletions.
diff --git a/fortran/DarkEnergyFluid.f90 b/fortran/DarkEnergyFluid.f90
@@ -91,10 +91,10 @@ subroutine TDarkEnergyFluid_Init(this, State)
             this%num_perturb_equations = 2
         end if
 	else
-		if (MGDE_const) then
-			this%num_perturb_equations = 0
-		else
+		if (MGDE_const .and. MGDE_pert) then
 			this%num_perturb_equations = 2
+		else
+			this%num_perturb_equations = 0
 		end if
 	end if
     end subroutine TDarkEnergyFluid_Init

diff --git a/fortran/DarkEnergyPPF.f90 b/fortran/DarkEnergyPPF.f90
@@ -65,10 +65,10 @@ subroutine TDarkEnergyPPF_Init(this, State)
             this%num_perturb_equations = 1
         end if
 	else
-		if (MGDE_const) then
-			this%num_perturb_equations = 0
-		else
+		if (MGDE_const .and. MGDE_pert) then
 			this%num_perturb_equations = 1
+		else
+			this%num_perturb_equations = 0
 		end if
 	end if
     if (this%cs2_lam /= 1._dl) &