update miniapps

PTsolvers · Jan 7, 2025 · 4ec4b2a · 4ec4b2a
1 parent 5905ef5
commit 4ec4b2a
Show file tree

Hide file tree

Showing 4 changed files with 9 additions and 7 deletions.
diff --git a/miniapps/benchmarks/stokes2D/shear_band/ShearBand2D.jl b/miniapps/benchmarks/stokes2D/shear_band/ShearBand2D.jl
@@ -74,24 +74,22 @@ function main(igg; nx=64, ny=64, figdir="model_figs")
             Phase             = 1,
             Density           = ConstantDensity(; ρ = 0.0),
             Gravity           = ConstantGravity(; g = 0.0),
-            # CompositeRheology = CompositeRheology((visc, el_bg, )),
                         CompositeRheology = CompositeRheology((visc, el_bg, pl)),
             Elasticity        = el_bg,
 
         ),
         # High density phase
         SetMaterialParams(;
-            Phase             = 1,
+            Phase             = 2,
             Density           = ConstantDensity(; ρ = 0.0),
             Gravity           = ConstantGravity(; g = 0.0),
-            # CompositeRheology = CompositeRheology((visc, el_inc, )),
-                        CompositeRheology = CompositeRheology((visc, el_inc, pl)),
+            CompositeRheology = CompositeRheology((visc, el_inc, pl)),
             Elasticity        = el_inc,
         ),
     )
 
     # perturbation array for the cohesion
-    perturbation_C = @rand(ni...)
+    perturbation_C = @zeros(ni...)
 
     # Initialize phase ratios -------------------------------
     radius       = 0.1

diff --git a/miniapps/benchmarks/stokes2D/shear_band/ShearBand2D_MPI.jl b/miniapps/benchmarks/stokes2D/shear_band/ShearBand2D_MPI.jl
@@ -68,6 +68,7 @@ function main(igg; nx=64, ny=64, figdir="model_figs")
         η_vp = η_reg,
         Ψ    = 0,
     )
+
     rheology = (
         # Low density phase
         SetMaterialParams(;
@@ -80,6 +81,7 @@ function main(igg; nx=64, ny=64, figdir="model_figs")
         ),
         # High density phase
         SetMaterialParams(;
+            Phase             = 2,
             Density           = ConstantDensity(; ρ = 0.0),
             Gravity           = ConstantGravity(; g = 0.0),
             CompositeRheology = CompositeRheology((visc, el_inc, pl)),

diff --git a/miniapps/benchmarks/stokes2D/shear_band/ShearBand2D_displacement.jl b/miniapps/benchmarks/stokes2D/shear_band/ShearBand2D_displacement.jl
@@ -81,7 +81,7 @@ function main(igg; nx=64, ny=64, figdir="model_figs")
         ),
         # High density phase
         SetMaterialParams(;
-            Phase             = 1,
+            Phase             = 2,
             Density           = ConstantDensity(; ρ = 0.0),
             Gravity           = ConstantGravity(; g = 0.0),
             CompositeRheology = CompositeRheology((visc, el_inc, pl)),
@@ -90,7 +90,7 @@ function main(igg; nx=64, ny=64, figdir="model_figs")
     )
 
     # perturbation array for the cohesion
-    perturbation_C = @rand(ni...)
+    perturbation_C = @zeros(ni...)
 
     # Initialize phase ratios -------------------------------
     radius       = 0.1

diff --git a/miniapps/benchmarks/stokes2D/shear_band/ShearBand2D_softening.jl b/miniapps/benchmarks/stokes2D/shear_band/ShearBand2D_softening.jl
@@ -71,6 +71,7 @@ function main(igg; nx=64, ny=64, figdir="model_figs")
         Ψ    = 0,
         softening_C = soft_C
     )
+
     rheology = (
         # Low density phase
         SetMaterialParams(;
@@ -83,6 +84,7 @@ function main(igg; nx=64, ny=64, figdir="model_figs")
         ),
         # High density phase
         SetMaterialParams(;
+            Phase             = 2,
             Density           = ConstantDensity(; ρ = 0.0),
             Gravity           = ConstantGravity(; g = 0.0),
             CompositeRheology = CompositeRheology((visc, el_inc, pl)),