diff --git a/src/inital.F90 b/src/inital.F90
index 7aa10a3497a839b44929a4eb2c66aeccdbe4b0ca..a42b90c67175f7b6d9b63aae3eef9a9c7ba0dd41 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