From e5be3acc0e29b202602f33968a4d6289fa2747f3 Mon Sep 17 00:00:00 2001
From: Antoine Cyril David Hoffmann <ahoffman@spcpc606.epfl.ch>
Date: Tue, 13 Apr 2021 15:33:05 +0200
Subject: [PATCH] comments and correction of Sapj computed before dnjs
---
src/inital.F90 | 20 +++++++++++---------
1 file changed, 11 insertions(+), 9 deletions(-)
diff --git a/src/inital.F90 b/src/inital.F90
index 7aa10a34..a42b90c6 100644
--- a/src/inital.F90
+++ b/src/inital.F90
@@ -25,21 +25,20 @@ SUBROUTINE inital
!!!!!! Set the moments arrays Nepj, Nipj and phi!!!!!!
IF ( RESTART ) THEN
IF (my_id .EQ. 0) WRITE(*,*) 'Load moments'
- CALL load_moments
+ CALL load_moments ! get N_0
IF (my_id .EQ. 0) WRITE(*,*) 'Init phi with Poisson'
- CALL poisson
- CALL MPI_BARRIER(MPI_COMM_WORLD,ierr);
+ CALL poisson ! compute phi_0=phi(N_0)
ELSE
IF (INIT_NOISY_PHI) THEN
IF (my_id .EQ. 0) WRITE(*,*) 'Init noisy phi'
- CALL init_phi
+ CALL init_phi ! init noisy phi_0, N_0 = 0
ELSE
IF (my_id .EQ. 0) WRITE(*,*) 'Init noisy moments and ghosts'
- CALL init_moments
+ CALL init_moments ! init noisy N_0
IF (my_id .EQ. 0) WRITE(*,*) 'Init phi with Poisson'
- CALL poisson
+ CALL poisson ! get phi_0 = phi(N_0)
ENDIF
ENDIF
@@ -52,15 +51,18 @@ SUBROUTINE inital
!!!!!! Set Sepj, Sipj and dnjs coeff table !!!!!!
IF ( NON_LIN ) THEN;
- IF (my_id .EQ. 0) WRITE(*,*) 'Init Sapj'
- CALL compute_Sapj
- ! WRITE(*,*) 'Building Dnjs table'
+ WRITE(*,*) 'Building Dnjs table'
CALL build_dnjs_table
+
+ IF (my_id .EQ. 0) WRITE(*,*) 'Init Sapj'
+ CALL compute_Sapj ! compute S_0 = S(phi_0,N_0)
ENDIF
!!!!!! Load the COSOlver collision operator coefficients !!!!!!
IF (ABS(CO) .GT. 1) THEN
CALL load_COSOlver_mat
+ ! Compute collision
+ CALL compute_TColl ! compute C_0 = C(N_0)
ENDIF
END SUBROUTINE inital
--
GitLab