diff --git a/src/inital.F90 b/src/inital.F90
index f6a7f2f8d06f15206a84f6cafb1121b37c90142b..19690e50fc97b59f71ef7240276605492bbae737 100644
--- a/src/inital.F90
+++ b/src/inital.F90
@@ -211,7 +211,7 @@ SUBROUTINE init_phi
         ENDDO
       ENDDO
 
-    ELSEIF(INIT_ZF_PHI .GT. 0)
+    ELSEIF(INIT_ZF_PHI .GT. 0) THEN
 
       !**** Zonal Flow initialization *******************************************
       ! Every mode is zero
@@ -222,7 +222,7 @@ SUBROUTINE init_phi
       END DO
       ! Except at ikr = mode number + 1, symmetry is already included since kr>=0
       phi(INIT_ZF_PHI+1,ikz_0) = -init_background * imagu/2._dp
-      
+
     ELSE ! we compute phi from noisy moments and poisson
 
       CALL poisson