From 055fe34b283a29e44580ef39cf3aec1f8c6752b9 Mon Sep 17 00:00:00 2001 From: Pascal Aellig Date: Wed, 25 Sep 2024 12:04:27 +0200 Subject: [PATCH] udpate JP --- Sill_2D_OM.jl | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/Sill_2D_OM.jl b/Sill_2D_OM.jl index 63daab77..46ed2585 100644 --- a/Sill_2D_OM.jl +++ b/Sill_2D_OM.jl @@ -153,7 +153,7 @@ function main2D(igg; ar=8, ny=16, nx=ny*8, figdir="figs2D", do_vtk =false) phases_dev = PTArray(backend_JR)(phases_GMG) phase_ratios = PhaseRatios(backend, length(rheology), ni); init_phases!(pPhases, phases_dev, particles, xvi) - phase_ratios_center!(phase_ratios, particles, xci, pPhases) + update_phase_ratios!(phase_ratios, particles, xci, xvi, pPhases) thermal = ThermalArrays(backend_JR, ni) @views thermal.T[2:end-1, :] .= PTArray(backend_JR)(nondimensionalize(T_GMG.*K, CharDim)) @@ -299,7 +299,7 @@ function main2D(igg; ar=8, ny=16, nx=ny*8, figdir="figs2D", do_vtk =false) # check if we need to inject particles inject_particles_phase!(particles, pPhases, (pT, ), (T_WENO,), xvi) # update phase ratios - phase_ratios_center!(phase_ratios, particles, xci, pPhases) + update_phase_ratios!(phase_ratios, particles, xci, xvi, pPhases) @show it += 1 t += dt @show extrema(thermal.T)