diff --git a/src/MNH/adv_boundaries.f90 b/src/MNH/adv_boundaries.f90
index c5e4df7525046cbb35ad06b6d235a428ef7f30ec..43cbe4735f02658bafefd0366b1a7d71140df2b7 100644
--- a/src/MNH/adv_boundaries.f90
+++ b/src/MNH/adv_boundaries.f90
@@ -24,7 +24,6 @@ END INTERFACE
 !
 CHARACTER(LEN=4), DIMENSION(2), INTENT(IN)      :: HLBCX,HLBCY   ! X and Y-direc. LBC type
 REAL, DIMENSION(:,:,:),   INTENT(INOUT)         :: PFIELD
-!$acc declare present(PFIELD)
 !
 END SUBROUTINE ADV_BOUNDARIES_DEVICE1
 !
@@ -33,7 +32,6 @@ END SUBROUTINE ADV_BOUNDARIES_DEVICE1
 CHARACTER(LEN=4), DIMENSION(2), INTENT(IN)      :: HLBCX,HLBCY   ! X and Y-direc. LBC type
 REAL, DIMENSION(:,:,:),   INTENT(INOUT)         :: PFIELD
 REAL, DIMENSION(:,:,:),   INTENT(IN)            :: PFIELDI
-!$acc declare present(PFIELD,PFIELDI)
 !
 END SUBROUTINE ADV_BOUNDARIES_DEVICE2
 !
@@ -43,7 +41,6 @@ CHARACTER(LEN=4), DIMENSION(2), INTENT(IN)      :: HLBCX,HLBCY   ! X and Y-direc
 REAL, DIMENSION(:,:,:),   INTENT(INOUT)         :: PFIELD
 REAL, DIMENSION(:,:,:),   INTENT(IN)            :: PFIELDI
 CHARACTER(LEN=1),         INTENT(IN)            :: HFIELD  ! Field type
-!$acc declare present(PFIELD,PFIELDI)
 !
 END SUBROUTINE ADV_BOUNDARIES_DEVICE3
 !
diff --git a/src/MNH/advec_4th_order_aux.f90 b/src/MNH/advec_4th_order_aux.f90
index 4313d08eb63bdf9ea54f6b1000cda61964ac8122..c06782fd5f3fe1e999d4ff6d298dd6d393b1bcdb 100644
--- a/src/MNH/advec_4th_order_aux.f90
+++ b/src/MNH/advec_4th_order_aux.f90
@@ -18,9 +18,7 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCX ! X direction LBC type
 CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCY ! Y direction LBC type
 !
 REAL, DIMENSION(:,:,:), INTENT(OUT)   :: PMEANX, PMEANY ! fluxes
-!$acc declare present(PMEANX,PMEANY)
 REAL, DIMENSION(:,:,:), INTENT(IN)    :: PFIELDT  ! variable at t
-!$acc declare present(PFIELDT)
 INTEGER,                INTENT(IN)    :: KGRID    ! C grid localisation
 !
 TYPE(HALO2_ll), OPTIONAL, POINTER :: TPHALO2      ! halo2 for the field at t
diff --git a/src/MNH/advec_ppm_algo.f90 b/src/MNH/advec_ppm_algo.f90
index 6fc8e9945a301411b945cf56c5289aaed8e9915a..196b546cfcf7a5b75bc80dbfdb90cb46cec546d0 100644
--- a/src/MNH/advec_ppm_algo.f90
+++ b/src/MNH/advec_ppm_algo.f90
@@ -21,18 +21,14 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCY  ! Y direction LBC type
 CHARACTER (LEN=6),               INTENT(IN) :: HMET_ADV_SCHEME
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PFIELDT      ! variable at t
-!$acc declare present(PFIELDT)
 INTEGER,                INTENT(IN)  :: KGRID        ! C grid localisation
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PCRU, PCRV, PCRW ! Courant numbers 
-!$acc declare present(PCRU,PCRV,PCRW)
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRHODJ  ! density
-!$acc declare present(PRHODJ)
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRHOX1, PRHOX2
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRHOY1, PRHOY2
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRHOZ1, PRHOZ2
-!$acc declare present(PRHOX1,PRHOX2,PRHOY1,PRHOY2,PRHOZ1,PRHOZ2)
 REAL,                   INTENT(IN)  :: PTSTEP ! Time step model  
 REAL,                   INTENT(IN)  :: PTSTEP_PPM ! Time Step PPM
 TYPE (DATE_TIME),       INTENT(IN)  :: TPDTCUR ! current date and time
diff --git a/src/MNH/advec_weno_k_1_aux.f90 b/src/MNH/advec_weno_k_1_aux.f90
index 0af47fef5511a5272f47878153a3c32d8aa50043..b12e73b4535742c694e646d68f72fd4a57a6fcb1 100644
--- a/src/MNH/advec_weno_k_1_aux.f90
+++ b/src/MNH/advec_weno_k_1_aux.f90
@@ -1,6 +1,6 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !     ##############################
       MODULE MODI_ADVEC_WENO_K_1_AUX
@@ -20,7 +20,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on U grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRUCT ! contrav. comp. on MASS GRID
 !
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR ! output src term
-!$acc declare present(PSRC,PRUCT,PR)
 END SUBROUTINE UP_UX_DEVICE
 !
 FUNCTION UP_MX(PSRC, PRUCT) RESULT(PR)
@@ -35,7 +34,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on U grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRUCT ! contrav. comp. on MASS GRID
 !
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR ! output src term
-!$acc declare present(PSRC,PRUCT,PR)
 END SUBROUTINE UP_MX_DEVICE
 !
 FUNCTION UP_VY(PSRC, PRVCT) RESULT(PR)
@@ -50,7 +48,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on U grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRVCT ! contrav. comp. on MASS GRID
 !
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR ! output src term
-!$acc declare present(PSRC,PRVCT,PR)
 END SUBROUTINE UP_VY_DEVICE
 !
 FUNCTION UP_MY(PSRC, PRVCT) RESULT(PR)
@@ -65,7 +62,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on U grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRVCT ! contrav. comp. on MASS GRID
 !
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR ! output src term
-!$acc declare present(PSRC,PRVCT,PR)
 END SUBROUTINE UP_MY_DEVICE
 !
 FUNCTION UP_WZ(PSRC, PRWCT) RESULT(PR)
@@ -80,7 +76,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on U grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRWCT ! contrav. comp. on MASS GRID
 !
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR ! output src term
-!$acc declare present(PSRC,PRWCT,PR)
 END SUBROUTINE UP_WZ_DEVICE
 !
 FUNCTION UP_MZ(PSRC, PRWCT) RESULT(PR)
@@ -95,7 +90,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on U grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRWCT ! contrav. comp. on MASS GRID
 !
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR ! output src term
-!$acc declare present(PSRC,PRWCT,PR)
 END SUBROUTINE UP_MZ_DEVICE
 !
 END INTERFACE
diff --git a/src/MNH/advec_weno_k_2_aux.f90 b/src/MNH/advec_weno_k_2_aux.f90
index d6dfbe12a92fc51c419e0d6c3152124ea15c4536..7e355f164a71bbd99e51f283a4b0cc6d011aeb99 100644
--- a/src/MNH/advec_weno_k_2_aux.f90
+++ b/src/MNH/advec_weno_k_2_aux.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
 !MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
@@ -26,27 +26,21 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCX  ! X direction LBC type
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on U grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRUCT ! contrav. comp. on MASS GRID
-!$acc declare present(PSRC,PRUCT)
 !
 ! output source term
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR
-!$acc declare present(PR)
 #ifndef _OPENACC
 TYPE(HALO2_ll), OPTIONAL, POINTER :: TPHALO2      ! halo2 for the field at t
 #else
 REAL, DIMENSION(:,:), INTENT(IN) :: TWEST,TEAST
-!$acc declare copyin(TWEST,TEAST)
 !
 ! Work arrays
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZFPOS1, ZFPOS2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZFNEG1, ZFNEG2
-!$acc declare present(ZFPOS1,ZFPOS2,ZFNEG1,ZFNEG2)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZBPOS1, ZBPOS2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZBNEG1, ZBNEG2
-!$acc declare present(ZBPOS1,ZBPOS2,ZBNEG1,ZBNEG2)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZOMP1, ZOMP2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZOMN1, ZOMN2
-!$acc declare present(ZOMP1,ZOMP2,ZOMN1,ZOMN2)
 #endif
 !
 END SUBROUTINE ADVEC_WENO_K_2_UX
@@ -70,27 +64,21 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCX  ! X direction LBC type
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on U grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRUCT ! contrav. comp. on MASS GRID
-!$acc declare present(PSRC,PRUCT)
 !
 ! output source term
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR
-!$acc declare present(PR)
 #ifndef _OPENACC
 TYPE(HALO2_ll), OPTIONAL, POINTER :: TPHALO2      ! halo2 for the field at t
 #else
 REAL, DIMENSION(:,:), INTENT(IN) :: TWEST,TEAST
-!$acc declare copyin(TWEST,TEAST)
 !
 ! Work arrays
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZFPOS1, ZFPOS2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZFNEG1, ZFNEG2
-!$acc declare present(ZFPOS1,ZFPOS2,ZFNEG1,ZFNEG2)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZBPOS1, ZBPOS2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZBNEG1, ZBNEG2
-!$acc declare present(ZBPOS1,ZBPOS2,ZBNEG1,ZBNEG2)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZOMP1, ZOMP2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZOMN1, ZOMN2
-!$acc declare present(ZOMP1,ZOMP2,ZOMN1,ZOMN2)
 #endif
 !
 END SUBROUTINE ADVEC_WENO_K_2_MX
@@ -114,28 +102,22 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCY  ! Y direction LBC type
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on U grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRVCT ! contrav. comp. on MASS GRID
-!$acc declare present(PSRC,PRVCT)
 !
 !
 ! output source term
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR
-!$acc declare present(PR)
 #ifndef _OPENACC
 TYPE(HALO2_ll), OPTIONAL, POINTER :: TPHALO2      ! halo2 for the field at t
 #else
 REAL, DIMENSION(:,:), INTENT(IN) :: TNORTH,TSOUTH
-!$acc declare copyin(TNORTH,TSOUTH)
 !
 ! Work arrays
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZFPOS1, ZFPOS2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZFNEG1, ZFNEG2
-!$acc declare present(ZFPOS1,ZFPOS2,ZFNEG1,ZFNEG2)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZBPOS1, ZBPOS2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZBNEG1, ZBNEG2
-!$acc declare present(ZBPOS1,ZBPOS2,ZBNEG1,ZBNEG2)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZOMP1, ZOMP2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZOMN1, ZOMN2
-!$acc declare present(ZOMP1,ZOMP2,ZOMN1,ZOMN2)
 #endif
 !
 END SUBROUTINE ADVEC_WENO_K_2_VY
@@ -159,27 +141,21 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCY  ! Y direction LBC type
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on U grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRVCT ! contrav. comp. on MASS GRID
-!$acc declare present(PSRC,PRVCT)
 !
 ! output source term
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR
-!$acc declare present(PR)
 #ifndef _OPENACC
 TYPE(HALO2_ll), OPTIONAL, POINTER :: TPHALO2      ! halo2 for the field at t
 #else
 REAL, DIMENSION(:,:), INTENT(IN) :: TNORTH,TSOUTH
-!$acc declare copyin(TNORTH,TSOUTH)
 !
 ! Work arrays
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZFPOS1, ZFPOS2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZFNEG1, ZFNEG2
-!$acc declare present(ZFPOS1,ZFPOS2,ZFNEG1,ZFNEG2)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZBPOS1, ZBPOS2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZBNEG1, ZBNEG2
-!$acc declare present(ZBPOS1,ZBPOS2,ZBNEG1,ZBNEG2)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZOMP1, ZOMP2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZOMN1, ZOMN2
-!$acc declare present(ZOMP1,ZOMP2,ZOMN1,ZOMN2)
 #endif
 !
 END SUBROUTINE ADVEC_WENO_K_2_MY
@@ -191,7 +167,6 @@ FUNCTION WENO_K_2_WZ(PSRC, PRWCT) RESULT(PR)
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on W grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRWCT ! contrav. comp. on MASS GRID
-!$acc declare present(PSRC,PRWCT)
 !
 ! output source term
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: PR
@@ -204,22 +179,17 @@ SUBROUTINE WENO_K_2_WZ(PSRC, PRWCT, PR,                                &
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on W grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRWCT ! contrav. comp. on MASS GRID
-!$acc declare present(PSRC,PRWCT)
 !
 ! output source term
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR
-!$acc declare present(PR)
 !
 ! Work arrays
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZFPOS1, ZFPOS2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZFNEG1, ZFNEG2
-!$acc declare present(ZFPOS1,ZFPOS2,ZFNEG1,ZFNEG2)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZBPOS1, ZBPOS2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZBNEG1, ZBNEG2
-!$acc declare present(ZBPOS1,ZBPOS2,ZBNEG1,ZBNEG2)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZOMP1, ZOMP2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZOMN1, ZOMN2
-!$acc declare present(ZOMP1,ZOMP2,ZOMN1,ZOMN2)
 !
 END SUBROUTINE WENO_K_2_WZ
 #endif
@@ -231,7 +201,6 @@ FUNCTION WENO_K_2_MZ(PSRC, PRWCT) RESULT(PR)
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on MASS grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRWCT ! contrav. comp. on W grid
-!$acc declare present(PSRC,PRWCT)
 !
 ! output source term
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: PR
@@ -244,22 +213,17 @@ SUBROUTINE WENO_K_2_MZ(PSRC, PRWCT, PR,                                &
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on MASS grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRWCT ! contrav. comp. on W grid
-!$acc declare present(PSRC,PRWCT)
 !
 ! output source term
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR
-!$acc declare present(PR)
 !
 ! Work arrays
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZFPOS1, ZFPOS2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZFNEG1, ZFNEG2
-!$acc declare present(ZFPOS1,ZFPOS2,ZFNEG1,ZFNEG2)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZBPOS1, ZBPOS2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZBNEG1, ZBNEG2
-!$acc declare present(ZBPOS1,ZBPOS2,ZBNEG1,ZBNEG2)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZOMP1, ZOMP2
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: ZOMN1, ZOMN2
-!$acc declare present(ZOMP1,ZOMP2,ZOMN1,ZOMN2)
 !
 END SUBROUTINE WENO_K_2_MZ
 #endif
diff --git a/src/MNH/advec_weno_k_3_aux.f90 b/src/MNH/advec_weno_k_3_aux.f90
index d775355120e9495e8465c7479c74cb7989108f71..440a55e4323b8e26c4ce7f2c38f83a5da1cf6901 100644
--- a/src/MNH/advec_weno_k_3_aux.f90
+++ b/src/MNH/advec_weno_k_3_aux.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
 !MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
@@ -29,23 +29,17 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCX  ! X direction LBC type
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on U grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRUCT ! contrav. comp. on MASS GRID
-!$acc declare present(PSRC,PRUCT)
 !
 ! output source term
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR
-!$acc declare present(PR)
 #ifdef _OPENACC
 ! Work arrays
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZFPOS1, ZFPOS2, ZFPOS3
-!$acc declare present(ZFPOS1, ZFPOS2, ZFPOS3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZFNEG1, ZFNEG2, ZFNEG3
-!$acc declare present(ZFNEG1, ZFNEG2, ZFNEG3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZBPOS1, ZBPOS2, ZBPOS3
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZBNEG1, ZBNEG2, ZBNEG3
-!$acc declare present(ZBPOS1, ZBPOS2, ZBPOS3, ZBNEG1, ZBNEG2, ZBNEG3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZOMP1, ZOMP2, ZOMP3
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZOMN1, ZOMN2, ZOMN3
-!$acc declare present(ZOMP1, ZOMP2, ZOMP3, ZOMN1, ZOMN2, ZOMN3)
 #endif
 !
 END SUBROUTINE ADVEC_WENO_K_3_UX
@@ -69,23 +63,17 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCX  ! X direction LBC type
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on U grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRUCT ! contrav. comp. on MASS GRID
-!$acc declare present(PSRC,PRUCT)
 !
 ! output source term
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR
-!$acc declare present(PR)
 #ifdef _OPENACC
 ! Work arrays
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZFPOS1, ZFPOS2, ZFPOS3
-!$acc declare present(ZFPOS1, ZFPOS2, ZFPOS3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZFNEG1, ZFNEG2, ZFNEG3
-!$acc declare present(ZFNEG1, ZFNEG2, ZFNEG3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZBPOS1, ZBPOS2, ZBPOS3
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZBNEG1, ZBNEG2, ZBNEG3
-!$acc declare present(ZBPOS1, ZBPOS2, ZBPOS3, ZBNEG1, ZBNEG2, ZBNEG3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZOMP1, ZOMP2, ZOMP3
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZOMN1, ZOMN2, ZOMN3
-!$acc declare present(ZOMP1, ZOMP2, ZOMP3, ZOMN1, ZOMN2, ZOMN3)
 #endif
 !
 END SUBROUTINE ADVEC_WENO_K_3_MX
@@ -109,24 +97,18 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCY  ! Y direction LBC type
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on U grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRVCT ! contrav. comp. on MASS GRID
-!$acc declare present(PSRC,PRVCT)
 !
 !
 ! output source term
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR
-!$acc declare present(PR)
 #ifdef _OPENACC
 ! Work arrays
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZFPOS1, ZFPOS2, ZFPOS3
-!$acc declare present(ZFPOS1, ZFPOS2, ZFPOS3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZFNEG1, ZFNEG2, ZFNEG3
-!$acc declare present(ZFNEG1, ZFNEG2, ZFNEG3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZBPOS1, ZBPOS2, ZBPOS3
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZBNEG1, ZBNEG2, ZBNEG3
-!$acc declare present(ZBPOS1, ZBPOS2, ZBPOS3, ZBNEG1, ZBNEG2, ZBNEG3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZOMP1, ZOMP2, ZOMP3
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZOMN1, ZOMN2, ZOMN3
-!$acc declare present(ZOMP1, ZOMP2, ZOMP3, ZOMN1, ZOMN2, ZOMN3)
 #endif
 !
 END SUBROUTINE ADVEC_WENO_K_3_VY
@@ -150,23 +132,17 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCY  ! Y direction LBC type
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on U grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRVCT ! contrav. comp. on MASS GRID
-!$acc declare present(PSRC,PRVCT)
 !
 ! output source term
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR
-!$acc declare present(PR)
 #ifdef _OPENACC
 ! Work arrays
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZFPOS1, ZFPOS2, ZFPOS3
-!$acc declare present(ZFPOS1, ZFPOS2, ZFPOS3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZFNEG1, ZFNEG2, ZFNEG3
-!$acc declare present(ZFNEG1, ZFNEG2, ZFNEG3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZBPOS1, ZBPOS2, ZBPOS3
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZBNEG1, ZBNEG2, ZBNEG3
-!$acc declare present(ZBPOS1, ZBPOS2, ZBPOS3, ZBNEG1, ZBNEG2, ZBNEG3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZOMP1, ZOMP2, ZOMP3
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZOMN1, ZOMN2, ZOMN3
-!$acc declare present(ZOMP1, ZOMP2, ZOMP3, ZOMN1, ZOMN2, ZOMN3)
 #endif
 !
 END SUBROUTINE ADVEC_WENO_K_3_MY
@@ -178,7 +154,6 @@ FUNCTION WENO_K_3_WZ(PSRC, PRWCT) RESULT(PR)
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on W grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRWCT ! contrav. comp. on MASS GRID
-!$acc declare present(PSRC,PRVCT)
 !
 ! output source term
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: PR
@@ -192,23 +167,17 @@ SUBROUTINE WENO_K_3_WZ(PSRC, PRWCT, PR, &
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on W grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRWCT ! contrav. comp. on MASS GRID
-!$acc declare present(PSRC,PRVCT)
 !
 ! output source term
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR
-!$acc declare present(PR)
 !
 ! Work arrays
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZFPOS1, ZFPOS2, ZFPOS3
-!$acc declare present(ZFPOS1, ZFPOS2, ZFPOS3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZFNEG1, ZFNEG2, ZFNEG3
-!$acc declare present(ZFNEG1, ZFNEG2, ZFNEG3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZBPOS1, ZBPOS2, ZBPOS3
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZBNEG1, ZBNEG2, ZBNEG3
-!$acc declare present(ZBPOS1, ZBPOS2, ZBPOS3, ZBNEG1, ZBNEG2, ZBNEG3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZOMP1, ZOMP2, ZOMP3
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZOMN1, ZOMN2, ZOMN3
-!$acc declare present(ZOMP1, ZOMP2, ZOMP3, ZOMN1, ZOMN2, ZOMN3)
 !
 END SUBROUTINE WENO_K_3_WZ
 #endif
@@ -220,7 +189,6 @@ FUNCTION WENO_K_3_MZ(PSRC, PRWCT) RESULT(PR)
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on MASS grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRWCT ! contrav. comp. on W grid
-!$acc declare present(PSRC,PRWCT)
 !
 ! output source term
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)) :: PR
@@ -234,23 +202,17 @@ SUBROUTINE WENO_K_3_MZ(PSRC, PRWCT, PR, &
 !
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PSRC  ! variable on MASS grid at t
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PRWCT ! contrav. comp. on W grid
-!$acc declare present(PSRC,PRWCT)
 !
 ! output source term
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR
-!$acc declare present(PR)
 !
 ! Work arrays
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZFPOS1, ZFPOS2, ZFPOS3
-!$acc declare present(ZFPOS1, ZFPOS2, ZFPOS3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZFNEG1, ZFNEG2, ZFNEG3
-!$acc declare present(ZFNEG1, ZFNEG2, ZFNEG3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZBPOS1, ZBPOS2, ZBPOS3
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZBNEG1, ZBNEG2, ZBNEG3
-!$acc declare present(ZBPOS1, ZBPOS2, ZBPOS3, ZBNEG1, ZBNEG2, ZBNEG3)
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZOMP1, ZOMP2, ZOMP3
 REAL, DIMENSION(SIZE(PSRC,1),SIZE(PSRC,2),SIZE(PSRC,3)):: ZOMN1, ZOMN2, ZOMN3
-!$acc declare present(ZOMP1, ZOMP2, ZOMP3, ZOMN1, ZOMN2, ZOMN3)
 !
 END SUBROUTINE WENO_K_3_MZ
 #endif
diff --git a/src/MNH/advecuvw_4th.f90 b/src/MNH/advecuvw_4th.f90
index 5e14c78447af68936df2403c7340cf010e90fb08..64f4059b737b510f5e6bbc1bf53ea2a6bf524f44 100644
--- a/src/MNH/advecuvw_4th.f90
+++ b/src/MNH/advecuvw_4th.f90
@@ -20,13 +20,10 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCY ! Y direction LBC type
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PRUCT ! contravariant
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PRVCT !  components
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PRWCT ! of momentum
-!$acc declare present(PRUCT,PRVCT,PRWCT)
 !
 REAL, DIMENSION(:,:,:),   INTENT(IN) :: PUT, PVT, PWT        ! U,V,W at t
-!$acc declare present(PUT,PVT,PWT)
 !
 REAL, DIMENSION(:,:,:),   INTENT(INOUT) :: PRUS, PRVS, PRWS     ! Source terms
-!$acc declare present(PRUS, PRVS, PRWS)
 !
 TYPE(HALO2LIST_ll), POINTER :: TPHALO2LIST ! list for diffusion
 !
diff --git a/src/MNH/advecuvw_rk.f90 b/src/MNH/advecuvw_rk.f90
index 9563e0dacfb01126980d2e5f2c358e2a8b72bf37..83c97ae5482793e359a28284374b035c961e3911 100644
--- a/src/MNH/advecuvw_rk.f90
+++ b/src/MNH/advecuvw_rk.f90
@@ -34,31 +34,23 @@ CHARACTER(LEN=4),DIMENSION(2),INTENT(IN):: HLBCX, HLBCY  ! X- and Y-direc LBC
 REAL,                     INTENT(IN)    :: PTSTEP
 !
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PU , PV  , PW
-!$acc declare present(PU,PV,PW)
                                                   ! Variables to advect
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PUT, PVT , PWT
-!$acc declare present(PUT,PVT,PWT)
                                                   ! Variables at t
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PMXM_RHODJ
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PMYM_RHODJ
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PMZM_RHODJ
-!$acc declare present(PMXM_RHODJ,PMYM_RHODJ,PMZM_RHODJ)
                                                   !  metric coefficients
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PRUCT , PRVCT, PRWCT
-!$acc declare present(PRUCT,PRVCT,PRWCT)
                                                   ! Contravariant wind components
 REAL, DIMENSION(:,:,:),   INTENT(OUT)   :: PRUS_ADV , PRVS_ADV, PRWS_ADV
-!$acc declare present(PRUS_ADV,PRVS_ADV,PRWS_ADV)
                                                   ! Tendency due to advection
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PRUS_OTHER , PRVS_OTHER, PRWS_OTHER
-!$acc declare present(PRUS_OTHER,PRVS_OTHER,PRWS_OTHER)
 !                                                 ! tendencies from other processes
 #ifdef _OPENACC
 ! Work arrays
 REAL, DIMENSION(:,:,:) :: ZUT, ZVT, ZWT
-!$acc declare present(ZUT,ZVT,ZWT)
 REAL, DIMENSION(:,:,:,:) :: ZRUS,ZRVS,ZRWS
-!$acc declare present(ZRUS,ZRVS,ZRWS)
 #endif
 !
 !
diff --git a/src/MNH/advecuvw_weno_k.f90 b/src/MNH/advecuvw_weno_k.f90
index 2984522835ea0efa0b71bb47e8792197697a954e..01b66a841435d9c7f3db27fbfd986862a5cb629b 100644
--- a/src/MNH/advecuvw_weno_k.f90
+++ b/src/MNH/advecuvw_weno_k.f90
@@ -1,4 +1,4 @@
-!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 2013-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
 !MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
@@ -27,20 +27,16 @@ INTEGER,                  INTENT(IN)    :: KWENO_ORDER   ! Order of the WENO
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PRUCT ! contravariant
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PRVCT !  components
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PRWCT ! of momentum
-!$acc declare present(PRUCT,PRVCT,PRWCT)
 !
 REAL, DIMENSION(:,:,:),   INTENT(IN) :: PUT, PVT, PWT        ! U,V,W at t
-!$acc declare present(PUT,PVT,PWT)
 !
 REAL, DIMENSION(:,:,:),   INTENT(INOUT) :: PRUS, PRVS, PRWS     ! Source terms
-!$acc declare present(PRUS,PRVS,PRWS)
 !
 TYPE(HALO2LIST_ll), POINTER :: TPHALO2LIST ! list for diffusion
 !
 #ifdef _OPENACC
 ! Work arrays
 REAL, DIMENSION(:,:,:) :: ZMEAN, ZWORK
-!$acc declare present(ZMEAN,ZWORK)
 #endif
 !
 END SUBROUTINE ADVECUVW_WENO_K
diff --git a/src/MNH/contrav.f90 b/src/MNH/contrav.f90
index e4556009d3fc48d1f55fda0c24d28a5b4319fa67..c5d4157f199c1ff8d4b6c1397efbc5642a9a2faf 100644
--- a/src/MNH/contrav.f90
+++ b/src/MNH/contrav.f90
@@ -40,20 +40,16 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCY  ! Y direction LBC type
 REAL, DIMENSION(:,:,:),  INTENT(IN)    ::  PRUT       ! Cartesian comp along x
 REAL, DIMENSION(:,:,:),  INTENT(IN)    ::  PRVT       ! Cartesian comp along y
 REAL, DIMENSION(:,:,:),  INTENT(IN)    ::  PRWT       ! Cartesian comp along z
-!$acc declare present(PRUT,PRVT,PRWT)
 REAL, DIMENSION(:,:,:),  INTENT(IN)    ::  PDXX       ! Metric coefficients
 REAL, DIMENSION(:,:,:),  INTENT(IN)    ::  PDYY       ! Metric coefficients
 REAL, DIMENSION(:,:,:),  INTENT(IN)    ::  PDZZ       ! Metric coefficients
 REAL, DIMENSION(:,:,:),  INTENT(IN)    ::  PDZX       ! Metric coefficients
 REAL, DIMENSION(:,:,:),  INTENT(IN)    ::  PDZY       ! Metric coefficients
-!$acc declare present(PDXX,PDYY,PDZZ,PDZX,PDZY)
 REAL, DIMENSION(:,:,:),  INTENT(OUT)   ::  PRUCT      ! Contrav comp along x-bar
 REAL, DIMENSION(:,:,:),  INTENT(OUT)   ::  PRVCT      ! Contrav comp along y-bar
 REAL, DIMENSION(:,:,:),  INTENT(OUT)   ::  PRWCT      ! Contrav comp along z-bar
-!$acc declare present(PRUCT,PRVCT,PRWCT)
 INTEGER,                 INTENT(IN)    ::  KADV_ORDER ! Order of the advection scheme
 REAL, DIMENSION(:,:,:),  INTENT(OUT)   :: Z1,Z2       ! Work arrays
-!$acc declare present(Z1,Z2)
 LOGICAL, OPTIONAL,       INTENT(IN)    :: ODATA_ON_DEVICE ! Is some of the data on the accelerator device
 !
 !
diff --git a/src/MNH/emoist.f90 b/src/MNH/emoist.f90
index becf7b6d96dc43b344dbd6774be90adcbb7e6d74..baebb8a4f611679d02104bff7fadabc32fb01a4d 100644
--- a/src/MNH/emoist.f90
+++ b/src/MNH/emoist.f90
@@ -1,13 +1,8 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1995-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 turb 2006/05/18 13:07:25
-!-----------------------------------------------------------------
 !#################
 MODULE MODI_EMOIST
 !#################
@@ -35,7 +30,6 @@ REAL, DIMENSION(:,:,:),  INTENT(IN)  ::   PSRCM     ! Normalized 2dn_order
 REAL,DIMENSION(SIZE(PTHLM,1),SIZE(PTHLM,2),SIZE(PTHLM,3)):: PEMOIST ! result
 #else
 REAL,DIMENSION(SIZE(PTHLM,1),SIZE(PTHLM,2),SIZE(PTHLM,3)),  INTENT(OUT):: PEMOIST ! result
-!$acc declare present(PTHLM,PRM,PLOCPEXNM,PAMOIST,PSRCM,PEMOIST)
 #endif
 !
 #ifndef _OPENACC
diff --git a/src/MNH/etheta.f90 b/src/MNH/etheta.f90
index 59456b2adf77a197a62e8fd2364d061b85749e82..ac103d6680e24e57c5c906c88470ca0368631c51 100644
--- a/src/MNH/etheta.f90
+++ b/src/MNH/etheta.f90
@@ -1,13 +1,8 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1995-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 turb 2006/05/18 13:07:25
-!-----------------------------------------------------------------
 !#################
 MODULE MODI_ETHETA
 !#################
@@ -36,7 +31,6 @@ REAL, DIMENSION(:,:,:),  INTENT(IN)  ::   PSRCM     ! Normalized 2dn_order
 REAL,DIMENSION(SIZE(PTHLM,1),SIZE(PTHLM,2),SIZE(PTHLM,3)):: PETHETA ! result
 #else
 REAL,DIMENSION(SIZE(PTHLM,1),SIZE(PTHLM,2),SIZE(PTHLM,3)),  INTENT(OUT):: PETHETA ! result
-!$acc declare present(PTHLM,PRM,PLOCPEXNM,PATHETA,PSRCM,PETHETA)
 #endif
 !
 !
diff --git a/src/MNH/get_halo.f90 b/src/MNH/get_halo.f90
index 4d478ad730445f6398ed0bfcea08accae58e1f0b..b41bc5c05eeef16596af8f9ff5f621801840062a 100644
--- a/src/MNH/get_halo.f90
+++ b/src/MNH/get_halo.f90
@@ -38,7 +38,6 @@ INTERFACE
    SUBROUTINE GET_HALO_D(PSRC, HDIR, HNAME)
      !
      REAL, DIMENSION(:,:,:), INTENT(INOUT)   :: PSRC    ! variable at t
-     !$acc declare present (PSRC)
      CHARACTER(len=4), OPTIONAL :: HDIR ! to send only halo on X or Y direction
      character(len=*), optional, intent(in) :: HNAME ! Name of the field to be added
      !
diff --git a/src/MNH/gradient_m.f90 b/src/MNH/gradient_m.f90
index 7431b5121e0251de93177e13e8f64e2e1d125ea7..79590507a04ad3826629ec2c2daff00ad0ff2ffa 100644
--- a/src/MNH/gradient_m.f90
+++ b/src/MNH/gradient_m.f90
@@ -33,7 +33,6 @@ REAL, DIMENSION(:,:,:),  INTENT(IN)  :: PDZZ    ! metric coefficient dzz
 REAL, DIMENSION(:,:,:),  INTENT(IN)  :: PDZX    ! metric coefficient dzx
 !
 REAL, DIMENSION(SIZE(PA,1),SIZE(PA,2),SIZE(PA,3)),  INTENT(OUT) :: PGX_M_M ! result mass point
-!$acc declare present(PA,PDXX,PDZZ,PDZX,PGX_M_M)
 !
 END SUBROUTINE GX_M_M_DEVICE
 #endif
@@ -62,7 +61,6 @@ REAL, DIMENSION(:,:,:),  INTENT(IN)  :: PDZZ    ! metric coefficient dzz
 REAL, DIMENSION(:,:,:),  INTENT(IN)  :: PDZY    ! metric coefficient dzy
 !
 REAL, DIMENSION(SIZE(PA,1),SIZE(PA,2),SIZE(PA,3)),  INTENT(OUT) :: PGY_M_M ! result mass point
-!$acc declare present(PA,PDYY,PDZZ,PDZY,PGY_M_M)
 !
 END SUBROUTINE GY_M_M_DEVICE
 #endif
@@ -88,7 +86,6 @@ REAL, DIMENSION(:,:,:),  INTENT(IN)  :: PA      ! variable at the mass point
 REAL, DIMENSION(:,:,:),  INTENT(IN)  :: PDZZ    ! metric coefficient dzz
 !
 REAL, DIMENSION(SIZE(PA,1),SIZE(PA,2),SIZE(PA,3)),  INTENT(OUT) :: PGZ_M_M ! result mass point
-!$acc declare present(PA,PDZZ,PGY_M_M)
 !
 END SUBROUTINE GZ_M_M_DEVICE
 #endif
@@ -119,7 +116,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN)  :: PDZZ     ! d*zz
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PY       ! variable at mass localization
 !
 REAL, DIMENSION(SIZE(PY,1),SIZE(PY,2),SIZE(PY,3)),INTENT(OUT) :: PGX_M_U ! result at flux side
-!$acc declare present(PY,PDXX,PDZZ,PDZX,PGX_M_U)
 !
 END SUBROUTINE GX_M_U_DEVICE
 #endif
@@ -149,7 +145,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN)  :: PDZZ     !d*zz
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PY       ! variable at mass localization
 !
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PGY_M_V ! result at flux side
-!$acc declare present(PY,PDYY,PDZZ,PDZY,PGY_M_V)
 !
 END SUBROUTINE GY_M_V_DEVICE
 #endif
@@ -177,7 +172,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN)  :: PDZZ     !d*zz
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PY       ! variable at mass localization
 !
 REAL, DIMENSION(SIZE(PY,1),SIZE(PY,2),SIZE(PY,3)),INTENT(OUT) :: PGZ_M_W ! result at flux side
-!$acc declare present(PY,PDZZ,PGZ_M_W)
 !
 END SUBROUTINE GZ_M_W_DEVICE
 #endif
diff --git a/src/MNH/gradient_u.f90 b/src/MNH/gradient_u.f90
index 9cd43ea9c10def9c3319449a72b3b8de99ba6450..4878e8137682ed5af0b80fd7852b6cd760844998 100644
--- a/src/MNH/gradient_u.f90
+++ b/src/MNH/gradient_u.f90
@@ -1,13 +1,8 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 operators 2006/05/18 13:07:25
-!-----------------------------------------------------------------
 !     ######################
       MODULE MODI_GRADIENT_U
 !     ######################
@@ -38,7 +33,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN) :: PDZZ     ! metric coefficient dzz
 REAL, DIMENSION(:,:,:), INTENT(IN) :: PDZX     ! metric coefficient dzx
 !
 REAL, DIMENSION(SIZE(PA,1),SIZE(PA,2),SIZE(PA,3)), INTENT(OUT) :: PGX_U_M_DEVICE ! result mass point
-!$acc declare present(PA,PDXX,PDZZ,PDZX,PGX_U_M_DEVICE)
 !
 END SUBROUTINE GX_U_M_DEVICE
 #endif
@@ -68,7 +62,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN) :: PDZZ     ! metric coefficient dzz
 REAL, DIMENSION(:,:,:), INTENT(IN) :: PDZY     ! metric coefficient dzy
 !
 REAL, DIMENSION(SIZE(PA,1),SIZE(PA,2),SIZE(PA,3)), INTENT(OUT) :: PGY_U_UV_DEVICE ! result UV point
-!$acc declare present(PA,PDYY,PDZZ,PDZY,PGY_U_UV_DEVICE)
 !
 END SUBROUTINE GY_U_UV_DEVICE
 #endif
@@ -94,7 +87,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN) :: PA       ! variable at the U point
 REAL, DIMENSION(:,:,:), INTENT(IN) :: PDZZ     ! metric coefficient dzz
 !
 REAL, DIMENSION(SIZE(PA,1),SIZE(PA,2),SIZE(PA,3)), INTENT(OUT) :: PGZ_U_UW_DEVICE ! result UW point
-!$acc declare present(PA,PDZZ,PGZ_U_UW_DEVICE)
 !
 END SUBROUTINE GZ_U_UW_DEVICE
 #endif
diff --git a/src/MNH/gradient_v.f90 b/src/MNH/gradient_v.f90
index be2265b84e0383c46596b6b333d5d5c008d2278a..19bcf1a3716612d21f2cb3e86caa91dc887f5586 100644
--- a/src/MNH/gradient_v.f90
+++ b/src/MNH/gradient_v.f90
@@ -1,13 +1,8 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 operators 2006/05/18 13:07:25
-!-----------------------------------------------------------------
 !     ######################
       MODULE MODI_GRADIENT_V
 !     ######################
@@ -40,7 +35,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN) :: PDZZ     ! metric coefficient dzz
 REAL, DIMENSION(:,:,:), INTENT(IN) :: PDZY     ! metric coefficient dzy
 !
 REAL, DIMENSION(SIZE(PA,1),SIZE(PA,2),SIZE(PA,3)), INTENT(OUT) :: PGY_V_M_DEVICE ! result mass point
-!$acc declare present(PA,PDYY,PDZZ,PDZY,PGY_V_M_DEVICE)
 !
 END SUBROUTINE GY_V_M_DEVICE
 #endif
@@ -71,7 +65,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN) :: PDZZ     ! metric coefficient dzz
 REAL, DIMENSION(:,:,:), INTENT(IN) :: PDZX     ! metric coefficient dzx
 !
 REAL, DIMENSION(SIZE(PA,1),SIZE(PA,2),SIZE(PA,3)), INTENT(OUT) :: PGX_V_UV_DEVICE ! result UV point
-!$acc declare present(PA,PDXX,PDZZ,PDZX,PGX_V_UV_DEVICE)
 !
 END SUBROUTINE GX_V_UV_DEVICE
 #endif
@@ -98,7 +91,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN) :: PA       ! variable at the V point
 REAL, DIMENSION(:,:,:), INTENT(IN) :: PDZZ     ! metric coefficient dzz
 !
 REAL, DIMENSION(SIZE(PA,1),SIZE(PA,2),SIZE(PA,3)), INTENT(OUT) :: PGZ_V_VW_DEVICE ! result VW point
-!$acc declare present(PA,PDZZ,PGZ_V_VW_DEVICE)
 !
 END SUBROUTINE GZ_V_VW_DEVICE
 #endif
diff --git a/src/MNH/gradient_w.f90 b/src/MNH/gradient_w.f90
index 7d8b2dc6ea8b91aa15a1ec69c15515fe9b2055ba..d3daba6c4f3864e4a754c526971ee65a4b213aa9 100644
--- a/src/MNH/gradient_w.f90
+++ b/src/MNH/gradient_w.f90
@@ -1,13 +1,8 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 operators 2006/05/18 13:07:25
-!-----------------------------------------------------------------
 !     ######################
       MODULE MODI_GRADIENT_W
 !     ######################
@@ -36,7 +31,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN) :: PA       ! variable at the W point
 REAL, DIMENSION(:,:,:), INTENT(IN) :: PDZZ     ! metric coefficient dzz
 !
 REAL, DIMENSION(SIZE(PA,1),SIZE(PA,2),SIZE(PA,3)),  INTENT(OUT) :: PGZ_W_M_DEVICE ! result mass point
-!$acc declare present(PA,PDZZ,PGZ_W_M_DEVICE)
 !
 END SUBROUTINE GZ_W_M_DEVICE
 #endif
@@ -67,7 +61,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN) :: PDZZ     ! metric coefficient dzz
 REAL, DIMENSION(:,:,:), INTENT(IN) :: PDZX     ! metric coefficient dzx
 !
 REAL, DIMENSION(SIZE(PA,1),SIZE(PA,2),SIZE(PA,3)),  INTENT(OUT) :: PGX_W_UW_DEVICE ! result UW point
-!$acc declare present(PA,PDXX,PDZZ,PDZX,PGX_W_UW_DEVICE)
 !
 END SUBROUTINE GX_W_UW_DEVICE
 #endif
@@ -98,7 +91,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN) :: PDZZ     ! metric coefficient dzz
 REAL, DIMENSION(:,:,:), INTENT(IN) :: PDZY     ! metric coefficient dzy
 !
 REAL, DIMENSION(SIZE(PA,1),SIZE(PA,2),SIZE(PA,3)),  INTENT(OUT) :: PGY_W_VW_DEVICE ! result VW point
-!$acc declare present(PA,PDYY,PDZZ,PDZY,PGY_W_VW_DEVICE)
 !
 END SUBROUTINE GY_W_VW_DEVICE
 #endif
diff --git a/src/MNH/les_mean_subgrid.f90 b/src/MNH/les_mean_subgrid.f90
index 2c07a5c1e2c02edd3630470fb1665288ee138ab5..d54880e5b4beec2e7c834089778044953b427a51 100644
--- a/src/MNH/les_mean_subgrid.f90
+++ b/src/MNH/les_mean_subgrid.f90
@@ -14,7 +14,6 @@ INTERFACE LES_MEAN_SUBGRID
 REAL,    DIMENSION(:,:,:), INTENT(IN)    :: PA
 !
 REAL,    DIMENSION(:,:,:), INTENT(INOUT) :: PA_MEAN
-!$acc declare present(PA,PA_MEAN)
 !
 LOGICAL, OPTIONAL,         INTENT(IN)    :: OSUM
 !
@@ -25,7 +24,6 @@ END SUBROUTINE LES_MEAN_SUBGRID_3D
 REAL,    DIMENSION(:,:), INTENT(IN)    :: PA
 !
 REAL,    DIMENSION(:),   INTENT(INOUT) :: PA_MEAN
-!$acc declare present(PA,PA_MEAN)
 !
 LOGICAL, OPTIONAL,       INTENT(IN)    :: OSUM
 !
diff --git a/src/MNH/les_ver_int.f90 b/src/MNH/les_ver_int.f90
index f87e87d153e3caa7b904778ca4ea975eb8192cde..935bd87116b85a3d605f98582abcf674a1729e3c 100644
--- a/src/MNH/les_ver_int.f90
+++ b/src/MNH/les_ver_int.f90
@@ -27,7 +27,6 @@ INTERFACE LES_VER_INT_DEVICE
 REAL,    DIMENSION(:,:,:), INTENT(IN)  :: PA_MNH
 !
 REAL,    DIMENSION(:,:,:), INTENT(OUT) :: PA_LES
-!$acc declare present(PA_MNH,PA_LES)
 !
 END SUBROUTINE LES_VER_INT_DEVICE
 !
diff --git a/src/MNH/ppm.f90 b/src/MNH/ppm.f90
index 2085af4a1c141d2c7be514272a260e4646a2024b..5a9a3659393fbc42d932d69a0932536112ac814e 100644
--- a/src/MNH/ppm.f90
+++ b/src/MNH/ppm.f90
@@ -62,7 +62,6 @@ REAL,                   INTENT(IN)  :: PTSTEP  ! Time step
 REAL, DIMENSION(SIZE(PCR,1),SIZE(PCR,2),SIZE(PCR,3)) :: PR
 #else
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PR
-!$acc declare present(PR)
 #endif
 !
 #ifndef _OPENACC
diff --git a/src/MNH/ppm_met.f90 b/src/MNH/ppm_met.f90
index 54ae39457e1c39fd8d7b1a454a8635dc6031bbe8..c14c37871dd2519b081d0c84b1542f5224b89cbf 100644
--- a/src/MNH/ppm_met.f90
+++ b/src/MNH/ppm_met.f90
@@ -28,36 +28,22 @@ TYPE (DATE_TIME),         INTENT(IN)    :: TPDTCUR ! current date and time
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PCRU  ! Courant
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PCRV  ! numbers
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PCRW  ! 
-!$acc declare present(PCRU,PCRV,PCRW)
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PRHODJ ! density
-!$acc declare present(PRHODJ)
 ! Temporary advected rhodj
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PRHOX1,PRHOX2
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PRHOY1,PRHOY2
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PRHOZ1,PRHOZ2
-!$acc declare present(PRHOX1,PRHOX2,PRHOY1,PRHOY2,PRHOZ1,PRHOZ2)
 !
 REAL,                     INTENT(IN)    :: PTSTEP ! Time step model  
 REAL,                     INTENT(IN)    :: PTSTEP_PPM ! Time Step PPM
 !
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PTHT, PTKET        ! Vars at t
 REAL, DIMENSION(:,:,:,:), INTENT(IN)    :: PRT 
-!PW: bug workaround with PGI (tested up to 16.10)
-!PTKET is present but if zero-size => not detected as present
-!!$acc declare present(PTHT,PTKET,PRT)
-!$acc declare present(PTHT,PRT)
-!$acc declare pcopyin(PTKET)
 !
 REAL, DIMENSION(:,:,:),   INTENT(OUT) :: PRTHS, PRTKES! Source terms
-!PW: bug workaround with PGI (tested up to 16.10)
-!PRTKES is present but if zero-size => not detected as present
-!!$acc declare present(PRTHS,PRTKES)
-!$acc declare present(PRTHS)
-!$acc declare pcopyout(PRTKES)
 REAL, DIMENSION(:,:,:,:), INTENT(OUT) :: PRRS 
-!$acc declare present(PRRS)
 !
-END SUBROUTINE PPM_MET   
+END SUBROUTINE PPM_MET
 !
 END INTERFACE
 !
diff --git a/src/MNH/ppm_rhodj.f90 b/src/MNH/ppm_rhodj.f90
index d99ef1db98e41fa1298d46ef78f03a86040fe18d..47e269c1f2206ec90325346ec472fe97276f4bc3 100644
--- a/src/MNH/ppm_rhodj.f90
+++ b/src/MNH/ppm_rhodj.f90
@@ -19,18 +19,15 @@ CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCX ! X direction LBC type
 CHARACTER (LEN=4), DIMENSION(2), INTENT(IN) :: HLBCY ! Y direction LBC type
 !
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PCRU  ! Contravariants compon.
-REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PCRV  ! 
-REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PCRW  ! 
-!$acc declare present(PCRU,PCRV,PCRW)
+REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PCRV  !
+REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PCRW  !
 REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PRHODJ ! density
-!$acc declare present(PRHODJ)
 !
-REAL,                     INTENT(IN)    :: PTSTEP ! Single Time step 
+REAL,                     INTENT(IN)    :: PTSTEP ! Single Time step
 ! Temporary advected rhodj
 REAL, DIMENSION(:,:,:),   INTENT(OUT)   :: PRHOX1,PRHOX2
 REAL, DIMENSION(:,:,:),   INTENT(OUT)   :: PRHOY1,PRHOY2
 REAL, DIMENSION(:,:,:),   INTENT(OUT)   :: PRHOZ1,PRHOZ2
-!$acc declare present(PRHOX1,PRHOX2, PRHOY1,PRHOY2, PRHOZ1,PRHOZ2)
 !
 END SUBROUTINE PPM_RHODJ   
 !
diff --git a/src/MNH/prandtl.f90 b/src/MNH/prandtl.f90
index c1565caacb0d87765249ddd00c53509834cc0416..ba6a205024c8798e994f114db69cd5f52a2cebfc 100644
--- a/src/MNH/prandtl.f90
+++ b/src/MNH/prandtl.f90
@@ -68,11 +68,6 @@ REAL, DIMENSION(:,:,:,:), INTENT(OUT)::  PRED2RS3 ! Redelsperger number R*2_qsv
 REAL, DIMENSION(:,:,:), INTENT(OUT)  ::  PBLL_O_E! beta*Lk*Leps/tke
 REAL, DIMENSION(:,:,:), INTENT(OUT)  ::  PETHETA ! coefficient E_theta
 REAL, DIMENSION(:,:,:), INTENT(OUT)  ::  PEMOIST ! coefficient E_moist
-!$acc declare present(PDXX,PDYY,PDZZ,PDZX,PDZY,                                            &
-!$acc &               PTHVREF,PLOCPEXNM,PATHETA,PAMOIST,PLM,PLEPS,PTHLM,PTKEM,             &
-!$acc &               PRM,PSVM,PSRCM,                                                      &
-!$acc &               PREDTH1,PREDR1,PRED2TH3,PRED2R3,PRED2THR3,PREDS1,PRED2THS3,PRED2RS3, &
-!$acc &               PBLL_O_E,PETHETA,PEMOIST)
 !
 END SUBROUTINE PRANDTL
 !
diff --git a/src/MNH/rotate_wind.f90 b/src/MNH/rotate_wind.f90
index c71acc211e84a661521ed964dae2542a5786c462..0297b6cfd4f9283a971fd42fab48cd83edea7d5b 100644
--- a/src/MNH/rotate_wind.f90
+++ b/src/MNH/rotate_wind.f90
@@ -1,13 +1,8 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1995-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 turb 2006/05/18 13:07:25
-!-----------------------------------------------------------------
 !    #######################  
      MODULE MODI_ROTATE_WIND
 !    #######################
@@ -36,8 +31,6 @@ REAL, DIMENSION(:,:),   INTENT(OUT)  ::  PUSLOPE         ! wind component along
                                  ! the maximum slope direction
 REAL, DIMENSION(:,:),   INTENT(OUT)  ::  PVSLOPE         ! wind component along
                                  !  the direction normal to the maximum slope one
-!$acc declare present(PU,PV,PW,PDIRCOSXW, PDIRCOSYW, PDIRCOSZW,  &
-!$acc &               PCOSSLOPE,PSINSLOPE,PDXX,PDYY,PDZZ,PUSLOPE,PVSLOPE)
 !
 !-------------------------------------------------------------------------------
 !
diff --git a/src/MNH/shuman_device.f90 b/src/MNH/shuman_device.f90
index 46a2227806186f5e512a0745d5ecd08ee7392614..1d40497fe925afd05162e48c22304c1d1326b63a 100644
--- a/src/MNH/shuman_device.f90
+++ b/src/MNH/shuman_device.f90
@@ -1,13 +1,8 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 operators 2006/05/18 13:07:25
-!-----------------------------------------------------------------
 #ifdef _OPENACC
 !     #########################
       MODULE MODI_SHUMAN_DEVICE
@@ -18,25 +13,21 @@ INTERFACE
 SUBROUTINE DXF_DEVICE(PA,PDXF)
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PA     ! variable at flux side
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PDXF   ! result at mass localization
-!$acc declare present(PA,PDXM)
 END SUBROUTINE DXF_DEVICE
 !
 SUBROUTINE DXM_DEVICE(PA,PDXM)
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PA     ! variable at mass localization
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PDXM   ! result at flux side
-!$acc declare present(PA,PDXM)
 END SUBROUTINE DXM_DEVICE
 !
 SUBROUTINE DYF_DEVICE(PA,PDYF)
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PA     ! variable at flux side
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PDYF   ! result at mass localization
-!$acc declare present(PA,PDYF)
 END SUBROUTINE DYF_DEVICE
 !
 SUBROUTINE DYM_DEVICE(PA,PDYM)
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PA     ! variable at mass localization
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PDYM   ! result at flux side
-!$acc declare present(PA,PDYM)
 END SUBROUTINE DYM_DEVICE
 !
 SUBROUTINE DZF_DEVICE(KKA,KKU,KL,PA,PDZF)
@@ -44,7 +35,6 @@ INTEGER,                INTENT(IN)  :: KKA, KKU ! near ground and uppest atmosph
 INTEGER,                INTENT(IN)  :: KL       ! +1 if grid goes from ground to atmosphere top, -1 otherwise
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PA       ! variable at flux side
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PDZF     ! result at mass localization
-!$acc declare present(PA,PDZF)
 END SUBROUTINE DZF_DEVICE
 !
 SUBROUTINE DZM_DEVICE(KKA,KKU,KL,PA,PDZM)
@@ -52,31 +42,26 @@ INTEGER,                INTENT(IN)  :: KKA, KKU ! near ground and uppest atmosph
 INTEGER,                INTENT(IN)  :: KL       ! +1 if grid goes from ground to atmosphere top, -1 otherwise
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PA       ! variable at mass localization
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PDZM     ! result at flux side
-!$acc declare present(PA,PDZM)
 END SUBROUTINE DZM_DEVICE
 !
 SUBROUTINE MXF_DEVICE(PA,PMXF)
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PA     ! variable at flux side
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PMXF   ! result at mass localization
-!$acc declare present(PA,PMXF)
 END SUBROUTINE MXF_DEVICE
 !
 SUBROUTINE MXM_DEVICE(PA,PMXM)
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PA     ! variable at mass localization
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PMXM   ! result at flux localization
-!$acc declare present(PA,PMXM)
 END SUBROUTINE MXM_DEVICE
 !
 SUBROUTINE MYF_DEVICE(PA,PMYF)
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PA     ! variable at flux side
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PMYF   ! result at mass localization
-!$acc declare present(PA,PMYF)
 END SUBROUTINE MYF_DEVICE
 !
 SUBROUTINE MYM_DEVICE(PA,PMYM)
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PA     ! variable at mass localization
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PMYM   ! result at flux localization
-!$acc declare present(PA,PMYM)
 END SUBROUTINE MYM_DEVICE
 !
 SUBROUTINE MZF_DEVICE(KKA,KKU,KL,PA,PMZF)
@@ -84,13 +69,11 @@ INTEGER,                INTENT(IN)  :: KKA, KKU ! near ground and uppest atmosph
 INTEGER,                INTENT(IN)  :: KL       ! +1 if grid goes from ground to atmosphere top, -1 otherwise
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PA       ! variable at flux side
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PMZF     ! result at mass localization
-!$acc declare present(PA,PMZF)
 END SUBROUTINE MZF_DEVICE
 !
 SUBROUTINE MZM_DEVICE(PA,PMZM)
 REAL, DIMENSION(:,:,:), INTENT(IN)  :: PA     ! variable at mass localization
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PMZM   ! result at flux localization
-!$acc declare present(PA,PMZM)
 END SUBROUTINE MZM_DEVICE
 !
 END INTERFACE
diff --git a/src/MNH/tke_eps_sources.f90 b/src/MNH/tke_eps_sources.f90
index 8153adf2d1ccb5be911b25abea7eab02b15f52b9..2960671bd4d62378ddd459fa96dedae112589ccf 100644
--- a/src/MNH/tke_eps_sources.f90
+++ b/src/MNH/tke_eps_sources.f90
@@ -48,10 +48,6 @@ REAL, DIMENSION(:,:,:),  INTENT(INOUT)::  PRTHLS       ! Source of Theta_l
 REAL, DIMENSION(:,:,:),  INTENT(IN)   ::  PCOEF_DISS   ! 1/(Cph*Exner)
 REAL, DIMENSION(:,:,:),  INTENT(OUT)  ::  PTR          ! Transport prod. of TKE
 REAL, DIMENSION(:,:,:),  INTENT(OUT)  ::  PDISS        ! Dissipati prod. of TKE
-!$acc declare present(PTKEM,PLM,PLEPS,PDP,PTRH,              &
-!$acc &               PRHODJ,PDZZ,PDXX,PDYY,PDZX,PDZY,PZZ,   &
-!$acc &               PTP,PRTKES,PRTKESM, PRTHLS,PCOEF_DISS, &
-!$acc &               PTR,PDISS)
 !
 !
 !
diff --git a/src/MNH/tm06_h.f90 b/src/MNH/tm06_h.f90
index 85a368f444f4410d98ff86a96cced296f0ab8278..5fb42f61763dc675b4321d0e2da423e4835193b0 100644
--- a/src/MNH/tm06_h.f90
+++ b/src/MNH/tm06_h.f90
@@ -20,8 +20,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN)    :: PZZ       ! altitude of flux levels
 REAL, DIMENSION(:,:,:), INTENT(IN)    :: PFLXZ     ! heat flux
 REAL, DIMENSION(:,:),   INTENT(INOUT) :: PBL_DEPTH ! boundary layer height
 !
-!$acc declare present(PZZ,PFLXZ,PBL_DEPTH)
-!
 !-------------------------------------------------------------------------------
 !
 END SUBROUTINE TM06_H
diff --git a/src/MNH/tridiag_thermo.f90 b/src/MNH/tridiag_thermo.f90
index 11888cfd758ed969208d6337e7e03cc834ea2d59..34618dfd2f37a08d37fffa33123111772168c53b 100644
--- a/src/MNH/tridiag_thermo.f90
+++ b/src/MNH/tridiag_thermo.f90
@@ -1,13 +1,8 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 newsrc 2006/06/02 17:32:30
-!-----------------------------------------------------------------
 !     ###################
       MODULE MODI_TRIDIAG_THERMO
 !     ###################
@@ -29,8 +24,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN) :: PRHODJ  ! (dry rho)*J          at mass poi
 !
 REAL, DIMENSION(:,:,:), INTENT(OUT):: PVARP   ! variable at t+1      at mass point
 !
-!$acc declare present(PVARM,PF,PDFDDTDZ,PDZZ,PRHODJ,PVARP)
-!
 END SUBROUTINE TRIDIAG_THERMO
 !
 END INTERFACE
diff --git a/src/MNH/tridiag_tke.f90 b/src/MNH/tridiag_tke.f90
index a91fb9d1552b6b6bedf711e31f1b1d3c57c53a1e..7e9e12808567204a5582edcf8c1bf8f95d01ee01 100644
--- a/src/MNH/tridiag_tke.f90
+++ b/src/MNH/tridiag_tke.f90
@@ -1,13 +1,8 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 turb 2006/06/06 10:00:49
-!-----------------------------------------------------------------
 !     ##########################
       MODULE MODI_TRIDIAG_TKE
 !     ##########################
@@ -30,8 +25,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN)  :: PDIAG       ! diagonal term linked to
 !
 REAL, DIMENSION(:,:,:), INTENT(OUT) :: PVARP       ! variable at t+1       
 !
-!$acc declare present(PVARM,PA,PRHODJ,PSOURCE,PDIAG,PVARP)
-!
 END SUBROUTINE TRIDIAG_TKE
 !
 END INTERFACE
diff --git a/src/MNH/tridiag_w.f90 b/src/MNH/tridiag_w.f90
index 57eaef202f43ce5637ed8781f3874c082a04e30f..b3ac682f9f9fbd54a27ecd198828bf14610216a6 100644
--- a/src/MNH/tridiag_w.f90
+++ b/src/MNH/tridiag_w.f90
@@ -1,7 +1,8 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
+!-----------------------------------------------------------------
 !     ###################
       MODULE MODI_TRIDIAG_W
 !     ###################
@@ -19,8 +20,6 @@ REAL, DIMENSION(:,:,:), INTENT(IN) :: PRHODJ  ! (dry rho)*J          at mass poi
 !
 REAL, DIMENSION(:,:,:), INTENT(OUT):: PVARP   ! variable at t+1      at flux point
 !
-!$acc declare present(PVARM,PF,PDFDDWDZ,PMZF_DZZ,PRHODJ,PVARP)
-!
 END SUBROUTINE TRIDIAG_W
 !
 END INTERFACE
diff --git a/src/MNH/tridiag_wind.f90 b/src/MNH/tridiag_wind.f90
index 9f769b152e3547eaa2b6118a72a9d84b09a59b35..bda39c9331b18267f684bbd28ac2a9e09acaf577 100644
--- a/src/MNH/tridiag_wind.f90
+++ b/src/MNH/tridiag_wind.f90
@@ -1,13 +1,8 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 turb 2006/06/06 09:36:38
-!-----------------------------------------------------------------
 !     ########################
       MODULE MODI_TRIDIAG_WIND
 !     ########################
@@ -30,8 +25,6 @@ REAL, DIMENSION(:,:,:),    INTENT(IN)  :: PSOURCE     ! source term of PVAR
 !
 REAL, DIMENSION(:,:,:),    INTENT(OUT) :: PVARP       ! variable at t+1        
 !
-!$acc declare present(PVARM,PA,PCOEFS,PRHODJA,PSOURCE,PVARP )
-!
 END SUBROUTINE TRIDIAG_WIND
 !
 END INTERFACE
diff --git a/src/MNH/turb_hor.f90 b/src/MNH/turb_hor.f90
index a68194934cb88840011944590c6be2162e1d84c0..121b3ec0e4e16b41504e7db18e3caf03f9e93918 100644
--- a/src/MNH/turb_hor.f90
+++ b/src/MNH/turb_hor.f90
@@ -102,18 +102,6 @@ REAL, DIMENSION(:,:,:),   INTENT(INOUT) ::  PSIGS
                                   ! IN: Vertical part of Sigma_s at t
                                   ! OUT: Total Sigma_s at t
 !
-!$acc declare present(PK,PINV_PDXX,PINV_PDYY,PINV_PDZZ,PMZM_PRHODJ,  &
-!$acc &               PDXX,PDYY,PDZZ,PDZX,PDZY,PZZ,                  &
-!$acc &               PDIRCOSXW,PDIRCOSYW,PDIRCOSZW,                 &
-!$acc &               PCOSSLOPE,PSINSLOPE,                           &
-!$acc &               PRHODJ,PTHVREF,PSFTHM,PSFRM,                   &
-!$acc &               PCDUEFF,PTAU11M,PTAU12M,PTAU22M,PTAU33M,       &
-!$acc &               PUM,PVM,PWM,PTHLM,PRM,PSVM,PTKEM,PLM,PLEPS,    &
-!$acc &               PLOCPEXNM,PATHETA,PAMOIST,PSRCM,PFRAC_ICE,     &
-!$acc &               PUSLOPEM,PVSLOPEM,PSFSVM,                      &
-!$acc &               PDP,PTP,PRTHLS,PRSVS,PRRS,                     &
-!$acc &               PRUS,PRVS,PRWS,PSIGS)
-!
 END SUBROUTINE TURB_HOR
 !
 END INTERFACE
diff --git a/src/MNH/turb_hor_dyn_corr.f90 b/src/MNH/turb_hor_dyn_corr.f90
index cab69003bf88f99c51bb2ba0059481145cf0c8c0..5b01a5d144f3cae34664cd900627c768a114e893 100644
--- a/src/MNH/turb_hor_dyn_corr.f90
+++ b/src/MNH/turb_hor_dyn_corr.f90
@@ -70,13 +70,6 @@ REAL, DIMENSION(:,:,:),   INTENT(IN)    ::  PLM          ! Turb. mixing length
 REAL, DIMENSION(:,:,:),   INTENT(INOUT) ::  PRUS, PRVS, PRWS
 REAL, DIMENSION(:,:,:),   INTENT(INOUT) ::  PDP,PTP      ! TKE production terms
 !
-!$acc declare present(PINV_PDZZ,PDXX,PDYY,PDZZ,PDZX,PDZY,PZZ,PDIRCOSZW,       &
-!$acc &               PCOSSLOPE,PSINSLOPE,PRHODJ,PCDUEFF,                     &
-!$acc &               PTAU11M,PTAU12M,PTAU22M,PTAU33M,                        &
-!$acc &               PUM,PVM,PWM,PTHLM,PRM,PSVM,PUSLOPEM,PVSLOPEM,PTKEM,PLM, &
-!$acc &               PRUS,PRVS,PRWS,PDP,PTP)
-!
-!
 END SUBROUTINE TURB_HOR_DYN_CORR
 !
 END INTERFACE
diff --git a/src/MNH/turb_hor_splt.f90 b/src/MNH/turb_hor_splt.f90
index e74b0c3f87f0f0f93b230143d5f1c4394efc04c7..0b8d32d44f9653989f281a70e708bc7a804e7dd1 100644
--- a/src/MNH/turb_hor_splt.f90
+++ b/src/MNH/turb_hor_splt.f90
@@ -97,19 +97,6 @@ REAL, DIMENSION(:,:,:),   INTENT(INOUT) ::  PSIGS
                                   ! IN: Vertical part of Sigma_s at t
                                   ! OUT: Total Sigma_s at t
 !
-!$acc declare present(PDXX,PDYY,PDZZ,PDZX,PDZY,PZZ,                  &
-!$acc &               PDIRCOSXW,PDIRCOSYW,PDIRCOSZW,                 &
-!$acc &               PCOSSLOPE,PSINSLOPE,                           &
-!$acc &               PRHODJ,PTHVREF,                                & 
-!$acc &               PSFTHM,PSFRM,PSFSVM,                           &
-!$acc &               PCDUEFF,PTAU11M,PTAU12M,PTAU22M,PTAU33M,       &
-!$acc &               PUM,PVM,PWM,PUSLOPEM,PVSLOPEM,PTHLM,PRM,PSVM,  &
-!$acc &               PTKEM,PLM,PLEPS,                               &
-!$acc &               PLOCPEXNM,PATHETA,PAMOIST,PSRCM,PFRAC_ICE,     &
-!$acc &               PDP,PTP,PSIGS,                                 &
-!$acc &               PTRH,                                          &
-!$acc &               PRUS,PRVS,PRWS,PRTHLS,PRRS,PRSVS               )
-!
 END SUBROUTINE TURB_HOR_SPLT
 !
 END INTERFACE
diff --git a/src/MNH/turb_hor_sv_flux.f90 b/src/MNH/turb_hor_sv_flux.f90
index 4148d7f780dbf8a2966401ca7d3624122371427f..b892fd00b1a7b891ecc342274e7ccfec59fce4be 100644
--- a/src/MNH/turb_hor_sv_flux.f90
+++ b/src/MNH/turb_hor_sv_flux.f90
@@ -50,14 +50,6 @@ REAL, DIMENSION(:,:,:,:), INTENT(IN)    ::  PSVM         ! scalar var. at t-1
 !
 REAL, DIMENSION(:,:,:,:), INTENT(INOUT) ::  PRSVS        ! var. at t+1 -split-
 !
-!$acc declare present(PK,PINV_PDXX,PINV_PDYY,PINV_PDZZ,PMZM_PRHODJ,  &
-!$acc &               PDXX,PDYY,PDZZ,PDZX,PDZY,                      &
-!$acc &               PDIRCOSXW,PDIRCOSYW,                           &
-!$acc &               PRHODJ,PWM,                                    &
-!$acc &               PSFSVM,                                        &
-!$acc &               PSVM,                                          &
-!$acc &               PRSVS)
-!
 END SUBROUTINE TURB_HOR_SV_FLUX
 !
 END INTERFACE
diff --git a/src/MNH/turb_hor_thermo_corr.f90 b/src/MNH/turb_hor_thermo_corr.f90
index 551945145af507066cb79fa9e97080a5973ee8c8..19006775795f4bbe294be170e6c5b6e252c569e7 100644
--- a/src/MNH/turb_hor_thermo_corr.f90
+++ b/src/MNH/turb_hor_thermo_corr.f90
@@ -60,15 +60,6 @@ REAL, DIMENSION(:,:,:),   INTENT(INOUT) ::  PSIGS
                                   ! IN: Vertical part of Sigma_s at t
                                   ! OUT: Total Sigma_s at t
 !
-!$acc declare present(PINV_PDXX,PINV_PDYY,                           &
-!$acc &               PDXX,PDYY,PDZZ,PDZX,PDZY,                      &
-!$acc &               PTHVREF,                                       &
-!$acc &               PWM,PTHLM,PRM,                                 &
-!$acc &               PTKEM,PLM,PLEPS,                               &
-!$acc &               PLOCPEXNM,PATHETA,PAMOIST,PSRCM,               &
-!$acc &               PSIGS)
-!
-!
 END SUBROUTINE TURB_HOR_THERMO_CORR
 !
 END INTERFACE
diff --git a/src/MNH/turb_hor_thermo_flux.f90 b/src/MNH/turb_hor_thermo_flux.f90
index 896d992b5fecbe1caba87d4d7726874b6655779a..7af4d30f5056f4802b2d8999673e9e9274598b1e 100644
--- a/src/MNH/turb_hor_thermo_flux.f90
+++ b/src/MNH/turb_hor_thermo_flux.f90
@@ -67,15 +67,6 @@ REAL, DIMENSION(:,:,:),   INTENT(IN)    ::  PFRAC_ICE    ! ri fraction of rc+ri
 REAL, DIMENSION(:,:,:),   INTENT(INOUT) ::  PRTHLS
 REAL, DIMENSION(:,:,:,:), INTENT(INOUT) ::  PRRS         ! var. at t+1 -split-
 !
-!$acc declare present(PK,PINV_PDXX,PINV_PDYY,PINV_PDZZ,PMZM_PRHODJ,  &
-!$acc &               PDXX,PDYY,PDZZ,PDZX,PDZY,                      &
-!$acc &               PDIRCOSXW,PDIRCOSYW,                           &
-!$acc &               PRHODJ,                                        &
-!$acc &               PSFTHM,PSFRM,                                  &
-!$acc &               PWM,PTHLM,PRM,                                 &
-!$acc &               PATHETA,PAMOIST,PSRCM,PFRAC_ICE,               &
-!$acc &               PRTHLS,PRRS                                    )
-!
 END SUBROUTINE TURB_HOR_THERMO_FLUX
 !
 END INTERFACE
diff --git a/src/MNH/turb_hor_tke.f90 b/src/MNH/turb_hor_tke.f90
index 736392149ac18e1e534d0d5ee41a40241f60ed99..355023335742b3112a539487fab3320ee40f955c 100644
--- a/src/MNH/turb_hor_tke.f90
+++ b/src/MNH/turb_hor_tke.f90
@@ -1,13 +1,8 @@
-!MNH_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
+!MNH_LIC Copyright 1994-2019 CNRS, Meteo-France and Universite Paul Sabatier
 !MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
+!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
 !MNH_LIC for details. version 1.
 !-----------------------------------------------------------------
-!--------------- special set of characters for RCS information
-!-----------------------------------------------------------------
-! $Source$ $Revision$
-! MASDEV4_7 turb 2006/05/18 13:07:25
-!-----------------------------------------------------------------
 !    ####################  
      MODULE MODI_TURB_HOR_TKE
 !    ####################  
@@ -37,11 +32,6 @@ REAL, DIMENSION(:,:,:),   INTENT(IN) :: PRHODJ       ! density * grid volume
 REAL, DIMENSION(:,:,:),   INTENT(IN)    ::  PTKEM    ! TKE at time t- dt
 REAL, DIMENSION(:,:,:),   INTENT(OUT)   ::  PTRH     ! horizontal transport of Tke 
 !
-!$acc declare present(PDXX,PDYY,PDZZ,PDZX,PDZY,                      &
-!$acc &               PINV_PDXX, PINV_PDYY, PINV_PDZZ, PMZM_PRHODJ,  &
-!$acc &               PK, PRHODJ, PTKEM,                             &
-!$acc &               PTRH                                           )
-!
 END SUBROUTINE TURB_HOR_TKE
 !
 END INTERFACE
diff --git a/src/MNH/turb_hor_uv.f90 b/src/MNH/turb_hor_uv.f90
index fb1481548a886ffbe49f40bb461da34d966e26c5..abd31ff25be4f827d0e1aaebd3b432e5791ff5b0 100644
--- a/src/MNH/turb_hor_uv.f90
+++ b/src/MNH/turb_hor_uv.f90
@@ -66,11 +66,6 @@ REAL, DIMENSION(:,:),      INTENT(IN)   ::  PVSLOPEM     ! wind component along
 REAL, DIMENSION(:,:,:),   INTENT(INOUT) ::  PRUS, PRVS   ! var. at t+1 -split-
 REAL, DIMENSION(:,:,:),   INTENT(INOUT) ::  PDP          ! TKE production terms
 !
-!$acc declare present(PK,PINV_PDXX,PINV_PDYY,PINV_PDZZ,PMZM_PRHODJ,PDXX,PDYY,PDZZ,PDZX,PDZY,        &
-!$acc &               PDIRCOSZW,PCOSSLOPE,PSINSLOPE,PRHODJ,PCDUEFF,PTAU11M,PTAU12M,PTAU22M,PTAU33M, &
-!$acc &               PUM,PVM,PUSLOPEM,PVSLOPEM,                                                    &
-!$acc &               PRUS,PRVS,PDP)
-!
 END SUBROUTINE TURB_HOR_UV
 !
 END INTERFACE
diff --git a/src/MNH/turb_hor_uw.f90 b/src/MNH/turb_hor_uw.f90
index 0b43466fa3961d3ca3ee732008735c533d0765fa..bfab5451144341db9d124c1c9f710b92fb7c4129 100644
--- a/src/MNH/turb_hor_uw.f90
+++ b/src/MNH/turb_hor_uw.f90
@@ -52,11 +52,6 @@ REAL, DIMENSION(:,:,:),   INTENT(IN)    ::  PLM          ! Turb. mixing length
 REAL, DIMENSION(:,:,:),   INTENT(INOUT) ::  PRUS, PRWS
 REAL, DIMENSION(:,:,:),   INTENT(INOUT) ::  PDP          ! TKE production terms
 !
-!$acc declare present(PK,PINV_PDXX,PINV_PDZZ,PMZM_PRHODJ,PDXX,PDZZ,PDZX, &
-!$acc &               PRHODJ,PTHVREF,                                    &
-!$acc &               PUM,PWM,PTHLM,PRM,PSVM,PTKEM,PLM,                  &
-!$acc &               PRUS,PRWS,PDP)
-!
 END SUBROUTINE TURB_HOR_UW
 !
 END INTERFACE
diff --git a/src/MNH/turb_hor_vw.f90 b/src/MNH/turb_hor_vw.f90
index 8e8bb89539fbfd31adcb3adc129146a39f156b97..9e1b3be389e5bcf33526a0d75430521b8c234b14 100644
--- a/src/MNH/turb_hor_vw.f90
+++ b/src/MNH/turb_hor_vw.f90
@@ -37,17 +37,13 @@ REAL, DIMENSION(:,:,:),   INTENT(IN)    ::  PINV_PDZZ   ! 1./PDZZ
 REAL, DIMENSION(:,:,:),   INTENT(IN)    ::  PMZM_PRHODJ ! MZM(PRHODJ)
 REAL, DIMENSION(:,:,:),   INTENT(IN)    ::  PDYY, PDZZ, PDZY 
                                                          ! Metric coefficients
-!$acc declare present(PDYY,PDZZ,PDZY)
 REAL, DIMENSION(:,:,:),   INTENT(IN)    ::  PRHODJ       ! density * grid volume
 REAL, DIMENSION(:,:,:),   INTENT(IN)    ::  PTHVREF      ! ref. state VPT       
 !
 ! Variables at t-1
 REAL, DIMENSION(:,:,:),   INTENT(IN)    ::  PVM,PWM,PTHLM
-!$acc declare copyin(PVM,PTHLM)
 REAL, DIMENSION(:,:,:,:), INTENT(IN)    ::  PRM
-!$acc declare copyin(PRM)
 REAL, DIMENSION(:,:,:,:), INTENT(IN)    ::  PSVM
-!$acc declare copyin(PSVM)
 !
 REAL, DIMENSION(:,:,:),   INTENT(IN)    ::  PTKEM        ! TKE at time t- dt
 REAL, DIMENSION(:,:,:),   INTENT(IN)    ::  PLM          ! Turb. mixing length
@@ -55,11 +51,6 @@ REAL, DIMENSION(:,:,:),   INTENT(IN)    ::  PLM          ! Turb. mixing length
 REAL, DIMENSION(:,:,:),   INTENT(INOUT) ::  PRVS, PRWS   ! var. at t+1 -split-
 REAL, DIMENSION(:,:,:),   INTENT(INOUT) ::  PDP          ! TKE production terms
 !
-!$acc declare present(PK,PINV_PDYY,PINV_PDZZ,PMZM_PRHODJ,PDYY,PDZZ,PDZY, &
-!$acc &               PRHODJ,PTHVREF,                                    &
-!$acc &               PVM,PWM,PTHLM,PRM,PSVM,PTKEM,PLM,                  &
-!$acc &               PRVS,PRWS,PDP)
-!
 END SUBROUTINE TURB_HOR_VW
 !
 END INTERFACE
diff --git a/src/MNH/turb_ver.f90 b/src/MNH/turb_ver.f90
index 8c3926d9f06cc8434c00fef956c5535678d0a854..29bcf15f8cc82cef3c5c330e0b04e579ce68192a 100644
--- a/src/MNH/turb_ver.f90
+++ b/src/MNH/turb_ver.f90
@@ -109,19 +109,6 @@ REAL, DIMENSION(:,:,:), INTENT(OUT)  :: PWTH       ! heat flux
 REAL, DIMENSION(:,:,:), INTENT(OUT)  :: PWRC       ! cloud water flux
 REAL, DIMENSION(:,:,:,:),INTENT(OUT) :: PWSV       ! scalar flux
 !
-!$acc declare present(PDXX,PDYY,PDZZ,PDZX,PDZY,PDIRCOSZW,PZZ,       &
-!$acc &               PCOSSLOPE,PSINSLOPE,                          &
-!$acc &               PRHODJ,PTHVREF,                               &
-!$acc &               PSFTHM,PSFRM,PSFSVM,PSFTHP,PSFRP,PSFSVP,      &
-!$acc &               PCDUEFF,PTAU11M,PTAU12M,PTAU33M,              &
-!$acc &               PUM,PVM,PWM,PUSLOPEM,PVSLOPEM,PTHLM,PRM,PSVM, &
-!$acc &               PTKEM,PLM,PLEPS,                              &
-!$acc &               PLOCPEXNM,PATHETA,PAMOIST,PSRCM,PFRAC_ICE,    &
-!$acc &               PFWTH,PFWR,PFTH2,PFR2,PFTHR,PBL_DEPTH,        &
-!$acc &               PSBL_DEPTH,PLMO,                              &
-!$acc &               PRUS,PRVS,PRWS,PRTHLS,PRRS,PRSVS,             &
-!$acc &               PDP,PTP,PSIGS,PWTH,PWRC,PWSV                  )
-!
 END SUBROUTINE TURB_VER
 !
 END INTERFACE
diff --git a/src/MNH/turb_ver_dyn_flux.f90 b/src/MNH/turb_ver_dyn_flux.f90
index abbc04e77a6f17561965a2375e9426ab0131beb9..531fe8781d37896d02c0325eacccd9d89a37beed 100644
--- a/src/MNH/turb_ver_dyn_flux.f90
+++ b/src/MNH/turb_ver_dyn_flux.f90
@@ -77,15 +77,6 @@ REAL, DIMENSION(:,:,:), INTENT(INOUT)   ::  PRUS, PRVS, PRWS
 !
 REAL, DIMENSION(:,:,:), INTENT(OUT)  ::  PDP          ! Dynamic TKE production term
 !
-!$acc declare present(PDXX,PDYY,PDZZ,PDZX,PDZY,PDIRCOSZW,PZZ,       &
-!$acc &               PCOSSLOPE,PSINSLOPE,                          &
-!$acc &               PRHODJ,                                       &
-!$acc &               PCDUEFF,PTAU11M,PTAU12M,PTAU33M,              &
-!$acc &               PTHLM,PRM,PSVM,PUM,PVM,PWM,PUSLOPEM,PVSLOPEM, &
-!$acc &               PTKEM,PLM,PWU,PWV,                            &
-!$acc &               PRUS,PRVS,PRWS,                               &
-!$acc &               PDP                                           )
-!
 END SUBROUTINE TURB_VER_DYN_FLUX
 !
 END INTERFACE
diff --git a/src/MNH/turb_ver_thermo_corr.f90 b/src/MNH/turb_ver_thermo_corr.f90
index 372cc9f32d829f372b80f28c3c4abac7bd9426ce..1d532b081ba1afda21dce4ba5aed54817238a4c9 100644
--- a/src/MNH/turb_ver_thermo_corr.f90
+++ b/src/MNH/turb_ver_thermo_corr.f90
@@ -100,18 +100,6 @@ REAL, DIMENSION(:,:,:),   INTENT(IN)    :: PRP        ! guess of r at t+ deltat
 !
 REAL, DIMENSION(:,:,:),   INTENT(OUT)  ::  PSIGS     ! Vert. part of Sigma_s at t
 !
-!$acc declare present(PDXX,PDYY,PDZZ,PDZX,PDZY,PDIRCOSZW,           &
-!$acc &               PRHODJ,PTHVREF,                               &
-!$acc &               PSFTHM,PSFRM,PSFTHP,PSFRP,                    &
-!$acc &               PWM,PTHLM,PRM,PSVM,                           &
-!$acc &               PTKEM,PLM,PLEPS,                              &
-!$acc &               PLOCPEXNM,PATHETA,PAMOIST,PSRCM,              &
-!$acc &               PBETA, PSQRT_TKE, PDTH_DZ, PDR_DZ, PRED2TH3,  &
-!$acc &               PRED2R3, PRED2THR3, PBLL_O_E, PETHETA,        &
-!$acc &               PEMOIST, PREDTH1, PREDR1, PPHI3, PPSI3, PD,   &
-!$acc &               PFWTH,PFWR,PFTH2,PFR2,PFTHR,                  &
-!$acc &               PTHLP,PRP,PSIGS                          )
-!
 END SUBROUTINE TURB_VER_THERMO_CORR
 !
 END INTERFACE
diff --git a/src/MNH/turb_ver_thermo_flux.f90 b/src/MNH/turb_ver_thermo_flux.f90
index 78e1f44fb06cc27a79bcf13c988b159b6e444638..1828d660553886a3408840a446bc20fb6d01f95f 100644
--- a/src/MNH/turb_ver_thermo_flux.f90
+++ b/src/MNH/turb_ver_thermo_flux.f90
@@ -112,18 +112,6 @@ REAL, DIMENSION(:,:,:),   INTENT(OUT)   :: PTP       ! Dynamic and thermal
 REAL, DIMENSION(:,:,:),   INTENT(OUT)   :: PWTH       ! heat flux
 REAL, DIMENSION(:,:,:),   INTENT(OUT)   :: PWRC       ! cloud water flux
 !
-!$acc declare present(PDXX,PDYY,PDZZ,PDZX,PDZY,PDIRCOSZW,PZZ,       &
-!$acc &               PRHODJ,PTHVREF,                               &
-!$acc &               PSFTHM,PSFRM,PSFTHP,PSFRP,                    &
-!$acc &               PWM,PTHLM,PRM,PSVM,                           &
-!$acc &               PTKEM,PLM,PLEPS,                              &
-!$acc &               PLOCPEXNM,PATHETA,PAMOIST,PSRCM,PFRAC_ICE,    &
-!$acc &               PBETA, PSQRT_TKE, PDTH_DZ, PDR_DZ, PRED2TH3,  &
-!$acc &               PRED2R3, PRED2THR3, PBLL_O_E, PETHETA,        &
-!$acc &               PEMOIST, PREDTH1, PREDR1, PPHI3, PPSI3, PD,   &
-!$acc &               PFWTH,PFWR,PFTH2,PFR2,PFTHR,PBL_DEPTH,        &
-!$acc &               PWTHV,PRTHLS,PRRS,PTHLP,PRP,PTP,PWTH,PWRC     )
-!
 END SUBROUTINE TURB_VER_THERMO_FLUX
 !
 END INTERFACE
diff --git a/src/MNH/update_lm.f90 b/src/MNH/update_lm.f90
index 9edc121e329a7069f678ef9392b5968026953967..562afddd5e9ac7632179bc38b1600a94a15a4a89 100644
--- a/src/MNH/update_lm.f90
+++ b/src/MNH/update_lm.f90
@@ -15,7 +15,6 @@ CHARACTER(LEN=4), DIMENSION(2), INTENT(IN) :: HLBCY ! Y boundary type
 !
 REAL, DIMENSION(:,:,:), INTENT(INOUT) :: PLM   ! mixing length
 REAL, DIMENSION(:,:,:), INTENT(INOUT) :: PLEPS ! dissipative length
-!$acc declare present(PLM,PLEPS)
 !
 END SUBROUTINE UPDATE_LM
 !