diff --git a/src/SURFEX/lib_mpp.F90 b/src/SURFEX/lib_mpp.F90
index 254409eac7ceb6532457ea2c8801ef5d9269add1..004bdad0946c0548789f6bc5a0564e471ba4af39 100644
--- a/src/SURFEX/lib_mpp.F90
+++ b/src/SURFEX/lib_mpp.F90
@@ -1,6 +1,6 @@
-!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!SFX_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
-!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !SFX_LIC for details. version 1.
 !GLT_LIC The GELATO model is a seaice model used in stand-alone or embedded mode. 
 !GLT_LIC  It has been developed by Meteo-France. The holder of GELATO is Meteo-France.
@@ -36,6 +36,9 @@
 !GLT_LIC  Attempts to use this software in a way not complying with CeCILL-C license
 !GLT_LIC  may lead to prosecution. 
 !GLT_LIC 
+!!    MODIFICATIONS
+!!    --------------
+!!       P. Wautelet 22/01/2019  use standard FLUSH statement instead of non standard intrinsics
 MODULE lib_mpp
    !!======================================================================
    !!                       ***  MODULE  lib_mpp  ***
@@ -3084,10 +3087,10 @@ CONTAINS
          IF( PRESENT(cd9 ) )   WRITE(numout,*) cd9
          IF( PRESENT(cd10) )   WRITE(numout,*) cd10
       ENDIF
-                               CALL FLUSH(numout    )
-      IF( numstp     /= -1 )   CALL FLUSH(numstp    )
-      IF( numsol     /= -1 )   CALL FLUSH(numsol    )
-      IF( numevo_ice /= -1 )   CALL FLUSH(numevo_ice)
+                               FLUSH(unit=numout    )
+      IF( numstp     /= -1 )   FLUSH(unit=numstp    )
+      IF( numsol     /= -1 )   FLUSH(unit=numsol    )
+      IF( numevo_ice /= -1 )   FLUSH(unit=numevo_ice)
       !
       IF( cd1 == 'STOP' ) THEN
          IF(lwp) WRITE(numout,*)  'huge E-R-R-O-R : immediate stop'
@@ -3123,7 +3126,7 @@ CONTAINS
          IF( PRESENT(cd9 ) ) WRITE(numout,*) cd9
          IF( PRESENT(cd10) ) WRITE(numout,*) cd10
       ENDIF
-      CALL FLUSH(numout)
+      FLUSH(unit=numout)
       !
    END SUBROUTINE ctl_warn
 
diff --git a/src/SURFEX/mode_pos_surf.F90 b/src/SURFEX/mode_pos_surf.F90
index 4f199ec59472960cd6d348961edf95f294f9ea20..478ee2034e3b350f6b3e2d92e8770f506c3fba0b 100644
--- a/src/SURFEX/mode_pos_surf.F90
+++ b/src/SURFEX/mode_pos_surf.F90
@@ -1,6 +1,6 @@
-!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!SFX_LIC Copyright 2001-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
-!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !SFX_LIC for details. version 1.
 !!    ###############
       MODULE MODE_POS_SURF
@@ -52,6 +52,7 @@ END INTERFACE
 !!    MODIFICATIONS
 !!    --------------
 !!       I. Mallet  15/10/01     adaptation to MesoNH (F90 norm)
+!!       P. Wautelet 22/01/2019  use standard FLUSH statement instead of non standard intrinsics
 !------------------------------------------------------------------------------
 !
 IMPLICIT NONE
@@ -101,7 +102,7 @@ DO JFILE=1,2
         IF (PRESENT(KLUOUT)) THEN
           WRITE(KLUOUT,FMT=*) 'MODE_POS_SURF : error reading from unit ',&
                 KULNAM,' file ',HDNAML,' line ',YLINE
-          CALL FLUSH(KLUOUT)
+          FLUSH(unit=KLUOUT)
         ENDIF        
         CALL ABOR1_SFX('MODE_POS_SURF: read error in namelist file') 
       ELSE
diff --git a/src/SURFEX/modi_gltools_chkout.F90 b/src/SURFEX/modi_gltools_chkout.F90
index 12f3d1ee0a2c3f274f8f6a7fda6e8e3b126f43d1..66f5a44ffc9bf93e98d1d3600c681b3763c505f9 100644
--- a/src/SURFEX/modi_gltools_chkout.F90
+++ b/src/SURFEX/modi_gltools_chkout.F90
@@ -1,6 +1,6 @@
-!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!SFX_LIC Copyright 2003-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
-!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !SFX_LIC for details. version 1.
 !GLT_LIC The GELATO model is a seaice model used in stand-alone or embedded mode. 
 !GLT_LIC  It has been developed by Meteo-France. The holder of GELATO is Meteo-France.
@@ -55,6 +55,7 @@
 ! Modified: 2009/08 (D. Salas y Melia) Adapted to Gelato's new interface
 ! Modified: 2012/11 (D. Salas y Melia) parallelism & super-type
 ! Modified: 2014/01 (D. Salas y Melia) use standard wrivai for writing
+! Modified: 2019/01 (P. Wautelet) use standard FLUSH statement instead of non standard intrinsics
 !
 ! ------------------- BEGIN MODULE modi_gltools_chkout --------------------
 !
@@ -349,7 +350,7 @@ SUBROUTINE gltools_chkout( kdate,tpglt )
       '                 END OF gelato TIME STEP Nr =',tpglt%ind%cur
     WRITE(noutlu,*)  &
       '                ************************************'
-    CALL FLUSH(noutlu)
+    FLUSH(unit=noutlu)
   ENDIF
 !
 !