cap charecteristic dt for subgrid diffusion

JuliaGeodynamics · Mar 10, 2024 · aa1a486 · aa1a486
1 parent 2ebc10c
commit aa1a486
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