oops

PTsolvers · Sep 20, 2024 · 3f53a50 · 3f53a50
1 parent 45681d6
commit 3f53a50
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/docs/src/man/subduction2D/subduction2D.md b/docs/src/man/subduction2D/subduction2D.md
@@ -78,7 +78,7 @@ particle_args       = (pT, pPhases)
 Now we assign the material phases from the arrays we computed with help of [GeophysicalModelGenerator.jl](https://github.com/JuliaGeodynamics/GeophysicalModelGenerator.jl)
 ```julia
 phases_device    = PTArray(backend)(phases_GMG)
-phase_ratios     = phase_ratios = PhaseRatios(backend, length(rheology), ni);
+phase_ratios     = PhaseRatios(backend, length(rheology), ni);
 init_phases!(pPhases, phases_device, particles, xvi)
 phase_ratios_center!(phase_ratios, particles, xci, pPhases)
 ```

diff --git a/miniapps/benchmarks/stokes3D/shear_heating/Shearheating3D.jl b/miniapps/benchmarks/stokes3D/shear_heating/Shearheating3D.jl
@@ -73,8 +73,8 @@ function main3D(igg; ar=8, ny=16, nx=ny*8, nz=ny*8, figdir="figs3D", do_vtk =fal
     yc_anomaly       = ly/2   # origin of thermal anomaly
     zc_anomaly       = 40e3  # origin of thermal anomaly
     r_anomaly        = 3e3    # radius of perturbation
-    init_phases!(backend_JR, pPhases, particles, xc_anomaly, yc_anomaly, zc_anomaly, r_anomaly)
-    phase_ratios     = PhaseRatio(ni, length(rheology))
+    init_phases!(pPhases, particles, xc_anomaly, yc_anomaly, zc_anomaly, r_anomaly)
+    phase_ratios     = PhaseRatio(backend, length(rheology), ni)
     phase_ratios_center!(phase_ratios, particles, xci, pPhases)
     # ----------------------------------------------------