Merge pull request #103 from scog1234/timestep_bugfix

Avoid use of uninitialized variables when setting the timestep in dense gas
grackle-project · Apr 25, 2022 · 121e980 · 121e980
2 parents abb942c + ac0d095
commit 121e980
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/clib/solve_rate_cool_g.F b/src/clib/solve_rate_cool_g.F
@@ -591,7 +591,8 @@ subroutine solve_rate_cool_g(icool, d, e, u, v, w, de,
      &              dt-ttot(i), 0.5_DKIND*dt)
 
                if (d(i,j,k)*dom .gt. 1e8_DKIND .and. 
-     &              edot(i) .gt. 0._DKIND) then
+     &              edot(i) .gt. 0._DKIND      .and.
+     &             ispecies .gt. 1) then
 !              Equilibrium value for H is:
 !              H = (-1._DKIND / (4*k22)) * (k13 - sqrt(8 k13 k22 rho + k13^2))
 !              We now want this to change by 10% or less, but we're only