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

STATISTICS: Calculate current density only if magnetic field present.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2022-12-17 15:09:57 -03:00
parent 4f8c9219b2
commit d9a7cc1b9e
1 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
sources/statistics.F90
View File

@ -772,10 +772,6 @@ module statistics
dh(2) = ady(pmeta%level) dh(2) = ady(pmeta%level)
dh(3) = adz(pmeta%level) dh(3) = adz(pmeta%level)
! calculate current density (J = ∇xB)
!
call curl(dh(:), pdata%q(ibx:ibz,:,:,:), jc(1:3,:,:,:))
! total mass ! total mass
! !
#if NDIMS == 3 #if NDIMS == 3
@ -820,6 +816,11 @@ module statistics
+ pdata%q(iby,nb:ne,nb:ne, : )**2 & + pdata%q(iby,nb:ne,nb:ne, : )**2 &
+ pdata%q(ibz,nb:ne,nb:ne, : )**2) * dvolh + pdata%q(ibz,nb:ne,nb:ne, : )**2) * dvolh
#endif /* NDIMS == 3 */ #endif /* NDIMS == 3 */
! calculate current density (J = ∇xB)
!
call curl(dh(:), pdata%q(ibx:ibz,:,:,:), jc(1:3,:,:,:))
end if end if
if (.not. periodic(1)) then if (.not. periodic(1)) then