From 40507846c90d53465d2008e8bc857ebae3705ea2 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Wed, 30 Oct 2019 14:47:42 +0100 Subject: [PATCH] cell displacement are between nodes which start at (0,0,0) --- src/grid/spectral_utilities.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/grid/spectral_utilities.f90 b/src/grid/spectral_utilities.f90 index 1a3aefa35..b36820ec8 100644 --- a/src/grid/spectral_utilities.f90 +++ b/src/grid/spectral_utilities.f90 @@ -1102,7 +1102,7 @@ subroutine utilities_updateCoords(F) ! calculate cell center displacements do k = 1,grid3; do j = 1,grid(2); do i = 1,grid(1) IPcoords(1:3,i,j,k) = vectorField_real(1:3,i,j,k) & - + matmul(Favg,step*real([i,j,k+grid3Offset]-1,pReal)) + + matmul(Favg,step*real([i,j,k+grid3Offset]-0.5_pReal,pReal)) enddo; enddo; enddo call discretization_setNodeCoords(reshape(NodeCoords,[3,(grid(1)+1)*(grid(2)+1)*(grid3+1)]))