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

INTEGRALS: Add forcing statistics.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2020-04-20 13:14:35 -03:00
parent b76eec7022
commit a50a4232ed
1 changed files with 10 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
sources/integrals.F90
View File

@ -179,9 +179,10 @@ module integrals
! write the integral file header
!
write(funit,"('#',a8,10(1x,a18))") 'step', 'time', 'dt' &
write(funit,"('#',a8,13(1x,a18))") 'step', 'time', 'dt' &
, 'mass', 'momx', 'momy', 'momz' &
, 'ener', 'ekin', 'emag', 'eint'
, 'ener', 'ekin', 'emag', 'eint' &
, 'einj', 'erat', 'arms'
write(funit,"('#')")
! depending on the append parameter, choose the right file initialization for
@ -321,6 +322,7 @@ module integrals
use equations , only : ien, imx, imy, imz
use equations , only : magnetized, gamma, csnd
use evolution , only : step, time, dtn
use forcing , only : einj, rinj, arms
use mpitools , only : master
#ifdef MPI
use mpitools , only : reduce_sum_real_array
@ -458,6 +460,11 @@ module integrals
#endif /* NDIMS == 3 */
end if
! sum up the injected energy and injection rate
!
inarr( 9) = einj
inarr(10) = rinj
! get average, minimum and maximum values of density
!
#if NDIMS == 3
@ -578,7 +585,7 @@ module integrals
! write down the integrals and statistics to appropriate files
!
if (master) then
write(funit,"(i9,10(1x,1es18.8e3))") step, time, dtn, inarr(1:8)
write(funit,"(i9,13(1x,1es18.8e3))") step, time, dtn, inarr(1:10), arms
write(sunit,"(i9,23(1x,1es18.8e3))") step, time &
, avarr(1), mnarr(1), mxarr(1) &
, avarr(2), mnarr(2), mxarr(2) &