From bd1b32d049428403e7f05d65c7d46ee7842224f8 Mon Sep 17 00:00:00 2001 From: Geoffroy Lesur Date: Tue, 5 Nov 2024 23:02:27 +0100 Subject: [PATCH] fix field computation for non-axisymmetric vector potential in spherical geometry that led to non-zero divergence --- src/dataBlock/dataBlockHost.cpp | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/dataBlock/dataBlockHost.cpp b/src/dataBlock/dataBlockHost.cpp index b5dff915..cfc83a3a 100644 --- a/src/dataBlock/dataBlockHost.cpp +++ b/src/dataBlock/dataBlockHost.cpp @@ -223,8 +223,8 @@ void DataBlockHost::MakeVsFromAmag(IdefixHostArray4D &Ain) { + 1/(x1m(i)*(cos(x2m(j)) - cos(x2m(j+1)))) * (sin(x2m(j+1))*Ain(KDIR,k,j+1,i) - sin(x2m(j))*Ain(KDIR,k,j,i) ) , - - 1/(x1m(i)*sin(x2(j))*dx3(k)) * (Ain(JDIR,k+1,j,i) - - Ain(JDIR,k,j,i) ) ); + - dx2(j)/(x1m(i)*((cos(x2m(j))- cos(x2m(j+1))))*dx3(k)) + * (Ain(JDIR,k+1,j,i) - Ain(JDIR,k,j,i) ) ); real Ax2m = fabs(sin(x2m(j))); // Regularisation along the axis