From 7ddc4acc256cd1f9aaafba5dc521f4527b569e0d Mon Sep 17 00:00:00 2001 From: albert-de-montserrat Date: Tue, 23 Apr 2024 13:47:39 +0200 Subject: [PATCH] format --- src/IO/VTK.jl | 2 +- src/stokes/VelocityKernels.jl | 4 ++-- 2 files changed, 3 insertions(+), 3 deletions(-) diff --git a/src/IO/VTK.jl b/src/IO/VTK.jl index 8aa11358..33d61017 100644 --- a/src/IO/VTK.jl +++ b/src/IO/VTK.jl @@ -89,4 +89,4 @@ function save_vtk(fname::String, xi, data::NamedTuple) end return nothing -end \ No newline at end of file +end diff --git a/src/stokes/VelocityKernels.jl b/src/stokes/VelocityKernels.jl index 47bc5a89..349ceace 100644 --- a/src/stokes/VelocityKernels.jl +++ b/src/stokes/VelocityKernels.jl @@ -206,12 +206,12 @@ end Vx[i + 1, j + 1, k + 1] += Rx_ijk * ηdτ / av_x(ητ) end if all((i, j, k) .< size(Vy) .- 1) - Ry_ijk = + Ry_ijk = Ry[i, j, k] = _dx * (τxy[i + 1, j + 1, k] - τxy[i, j + 1, k]) + _dy * (τyy[i, j + 1, k] - τyy[i, j, k]) + _dz * (τyz[i, j + 1, k + 1] - τyz[i, j + 1, k]) - d_ya(P) - av_y(fy) - Vy[i + 1, j + 1, k + 1] += Ry_ijk * ηdτ / av_y(ητ) + Vy[i + 1, j + 1, k + 1] += Ry_ijk * ηdτ / av_y(ητ) end if all((i, j, k) .< size(Vz) .- 1) Rz_ijk =