From b24ab3838938c350d264f9f61b65b25998403169 Mon Sep 17 00:00:00 2001 From: Grzegorz Kowal Date: Fri, 16 Aug 2024 16:37:11 -0300 Subject: [PATCH] USER_PROBLEM: Normalize the magnetic flux by the X area. Signed-off-by: Grzegorz Kowal --- sources/user_problem.F90 | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/sources/user_problem.F90 b/sources/user_problem.F90 index b19c6c5..c37c089 100644 --- a/sources/user_problem.F90 +++ b/sources/user_problem.F90 @@ -618,7 +618,7 @@ module user_problem #else /* NDIMS == 3 */ use coordinates, only : periodic, ymin, ymax #endif /* NDIMS == 3 */ - use coordinates, only : ax, adx, ady, adz + use coordinates, only : ax, adx, ady, adz, xarea #if NDIMS == 3 use equations , only : ivx, ivy, ivz, ibx, iby, ibz, ibp #else /* NDIMS == 3 */ @@ -915,6 +915,7 @@ module user_problem call reduce_sum(rterms(:)) #endif /* MPI */ + rterms(1) = rterms(1) / xarea rterms(6) = eta * rterms(6) rterms(7) = eta * rterms(7)