From db8bd3359ad52d5ad64c837537fde70e22cf82f8 Mon Sep 17 00:00:00 2001
From: Pascal Aellig <pascal.aellig@hotmail.com>
Date: Mon, 9 Sep 2024 15:31:07 +0200
Subject: [PATCH] =?UTF-8?q?adds=20pressure=20function=20w/=20melt=20fracti?=
 =?UTF-8?q?on=20&=20w/out=20=CE=94Tc?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

---
 src/stokes/PressureKernels.jl | 19 +++++++++++++++++++
 1 file changed, 19 insertions(+)

diff --git a/src/stokes/PressureKernels.jl b/src/stokes/PressureKernels.jl
index b42f7f901..d423c2b89 100644
--- a/src/stokes/PressureKernels.jl
+++ b/src/stokes/PressureKernels.jl
@@ -88,6 +88,25 @@ end
     return nothing
 end
 
+@parallel_indices (I...) function compute_P_kernel!(
+    P,
+    P0,
+    RP,
+    ∇V,
+    η,
+    rheology::NTuple{N,MaterialParams},
+    phase_ratio::C,
+    dt,
+    r,
+    θ_dτ,
+    ::Nothing,
+    ϕ,
+) where {N,C<:JustRelax.CellArray}
+    K = fn_ratio(get_bulk_modulus, rheology, phase_ratio[I...])
+    RP[I...], P[I...] = _compute_P!(P[I...], P0[I...], ∇V[I...], η[I...], K, dt, r, θ_dτ)
+    return nothing
+end
+
 @parallel_indices (I...) function compute_P_kernel!(
     P,
     P0,