Merge pull request #233 from PTsolvers/paellig-patch-2

Add missing Pressure function

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,