gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

INTEGRALS: Use primitive variables to calculate reconnection rate terms.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2020-12-25 11:47:02 -03:00
parent 62d9a6f97b
commit a30bfa0dd8
1 changed files with 12 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

24
sources/integrals.F90
View File

@ -684,9 +684,9 @@ module integrals
! the integral of |Bx| ! the integral of |Bx|
! !
#if NDIMS == 3 #if NDIMS == 3
inarr(11) = inarr(11) + sum(abs(pdata%u(ibx,nb:ne,nb:ne,nb:ne))) * dvol inarr(11) = inarr(11) + sum(abs(pdata%q(ibx,nb:ne,nb:ne,nb:ne))) * dvol
#else /* NDIMS == 3 */ #else /* NDIMS == 3 */
inarr(11) = inarr(11) + sum(abs(pdata%u(ibx,nb:ne,nb:ne, : ))) * dvol inarr(11) = inarr(11) + sum(abs(pdata%q(ibx,nb:ne,nb:ne, : ))) * dvol
#endif /* NDIMS == 3 */ #endif /* NDIMS == 3 */
! magnetic energy ! magnetic energy
@ -701,10 +701,10 @@ module integrals
! !
#if NDIMS == 3 #if NDIMS == 3
inarr(12) = inarr(12) & inarr(12) = inarr(12) &
+ 0.5d+00 * sum(abs(pdata%u(ibx,nb:ne,nb,nb:ne))) * dxz + 0.5d+00 * sum(abs(pdata%q(ibx,nb:ne,nb,nb:ne))) * dxz
#else /* NDIMS == 3 */ #else /* NDIMS == 3 */
inarr(12) = inarr(12) & inarr(12) = inarr(12) &
+ 0.5d+00 * sum(abs(pdata%u(ibx,nb:ne,nb, : ))) * dxz + 0.5d+00 * sum(abs(pdata%q(ibx,nb:ne,nb, : ))) * dxz
#endif /* NDIMS == 3 */ #endif /* NDIMS == 3 */
! advection of Bx along Y ! advection of Bx along Y
@ -765,12 +765,12 @@ module integrals
! !
#if NDIMS == 3 #if NDIMS == 3
inarr(22) = inarr(22) & inarr(22) = inarr(22) &
+ 2.5d-01 * sum((pdata%u(ibx,nb:ne,nb,nb:ne)**2 & + 2.5d-01 * sum((pdata%q(ibx,nb:ne,nb,nb:ne)**2 &
+ pdata%u(ibz,nb:ne,nb,nb:ne)**2)) * dxz + pdata%q(ibz,nb:ne,nb,nb:ne)**2)) * dxz
#else /* NDIMS == 3 */ #else /* NDIMS == 3 */
inarr(22) = inarr(22) & inarr(22) = inarr(22) &
+ 2.5d-01 * sum((pdata%u(ibx,nb:ne,nb, : )**2 & + 2.5d-01 * sum((pdata%q(ibx,nb:ne,nb, : )**2 &
+ pdata%u(ibz,nb:ne,nb, : )**2)) * dxz + pdata%q(ibz,nb:ne,nb, : )**2)) * dxz
#endif /* NDIMS == 3 */ #endif /* NDIMS == 3 */
! advection of magnetic field along Y ! advection of magnetic field along Y
@ -861,12 +861,12 @@ module integrals
! !
#if NDIMS == 3 #if NDIMS == 3
inarr(22) = inarr(22) & inarr(22) = inarr(22) &
+ 2.5d-01 * sum((pdata%u(ibx,nb:ne,ne,nb:ne)**2 & + 2.5d-01 * sum((pdata%q(ibx,nb:ne,ne,nb:ne)**2 &
+ pdata%u(ibz,nb:ne,ne,nb:ne)**2)) * dxz + pdata%q(ibz,nb:ne,ne,nb:ne)**2)) * dxz
#else /* NDIMS == 3 */ #else /* NDIMS == 3 */
inarr(22) = inarr(22) & inarr(22) = inarr(22) &
+ 2.5d-01 * sum((pdata%u(ibx,nb:ne,ne, : )**2 & + 2.5d-01 * sum((pdata%q(ibx,nb:ne,ne, : )**2 &
+ pdata%u(ibz,nb:ne,ne, : )**2)) * dxz + pdata%q(ibz,nb:ne,ne, : )**2)) * dxz
#endif /* NDIMS == 3 */ #endif /* NDIMS == 3 */
! advection of magnetic field along Y ! advection of magnetic field along Y