From 1192570d08be1fce479f28bc175e1e2145314e02 Mon Sep 17 00:00:00 2001 From: cflag Date: Mon, 9 Dec 2024 10:52:17 +0100 Subject: [PATCH] Fix force estimation --- src/forces.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/forces.f90 b/src/forces.f90 index 7ce681a2..f6656a9d 100644 --- a/src/forces.f90 +++ b/src/forces.f90 @@ -195,8 +195,8 @@ subroutine forces_unst(rhs_dux, rhs_duy, px, py, dt) do k = 1, xsize(3) do j = 1, xsize(2) do i = 1, xsize(1) - dux(i,j,k) = rhs_dux(i,j,k) - px(i,j,k) / dt - duy(i,j,k) = rhs_duy(i,j,k) - py(i,j,k) / dt + dux(i,j,k) = dt*rhs_dux(i,j,k) - px(i,j,k) + duy(i,j,k) = dt*rhs_duy(i,j,k) - py(i,j,k) end do end do end do