Merge pull request #69 from JuliaGeodynamics/albert-de-montserrat-pat…

…ch-4 Cap charecteristic dt for subgrid diffusion
JuliaGeodynamics · Mar 10, 2024 · 0faa439 · 0faa439
2 parents 2ebc10c + 27952a6
commit 0faa439
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/Physics/subgrid_diffusion.jl b/src/Physics/subgrid_diffusion.jl
@@ -65,7 +65,7 @@ end
         pTᵢ = @cell pT[ip, I...]
 
         # subgrid diffusion of the i-th particle
-        pΔTᵢ = (pTᵢ - pT0ᵢ) * (1 - exp(-d * dt / @cell(dt₀[ip, I...])))
+        pΔTᵢ = (pTᵢ - pT0ᵢ) * (1 - exp(-d * dt / max(@cell(dt₀[ip, I...]), 1e-9)))
         @cell pT0[ip, I...] = pT0ᵢ + pΔTᵢ
         @cell pΔT[ip, I...] = pΔTᵢ
     end