From 9c057a15efac6030d19ef7bc70f0c44f43146271 Mon Sep 17 00:00:00 2001 From: Antoine Hoffmann <antoine.hoffmann@epfl.ch> Date: Wed, 5 Jul 2023 09:21:08 +0200 Subject: [PATCH] compute transport at init avoid incoherent value in terminal --- src/inital.F90 | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) diff --git a/src/inital.F90 b/src/inital.F90 index edae542b..cfc12ecb 100644 --- a/src/inital.F90 +++ b/src/inital.F90 @@ -10,7 +10,7 @@ SUBROUTINE inital USE closure, ONLY: apply_closure_model USE ghosts, ONLY: update_ghosts_moments, update_ghosts_EM USE restarts, ONLY: load_moments, job2load - USE processing, ONLY: compute_fluid_moments + USE processing, ONLY: compute_fluid_moments, compute_radial_heatflux, compute_radial_transport USE model, ONLY: LINEARITY USE nonlinear, ONLY: compute_Sapj, nonlinear_init IMPLICIT NONE @@ -88,14 +88,17 @@ SUBROUTINE inital !! Preparing auxiliary arrays at initial state ! particle density, fluid velocity and temperature (used in diagnose) - CALL speak('Computing fluid moments') + CALL speak('Computing fluid moments and transport') CALL compute_fluid_moments + CALL compute_radial_transport + CALL compute_radial_heatflux ! init auxval for nonlinear CALL nonlinear_init ! compute nonlinear for t=0 diagnostic CALL compute_Sapj ! compute S_0 = S(phi_0,N_0) + END SUBROUTINE inital !******************************************************************************! -- GitLab