USER_PROBLEM: Normalize the magnetic flux by the X area.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2024-08-16 16:37:11 -03:00
parent ab481646b6
commit b24ab38389

View File

@ -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)