Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/com_io_dyn_mod.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/com_io_dyn_mod.F90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/com_io_dyn_mod.F90	(revision 5113)
@@ -4,5 +4,5 @@
 module com_io_dyn_mod
 
-  implicit none 
+  IMPLICIT NONE
 
 ! Names of various files for outputs (in the dynamics)
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/conf_dat_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/conf_dat_m.F90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/conf_dat_m.F90	(revision 5113)
@@ -24,5 +24,5 @@
 !  * interbar is TRUE for barycentric interpolation.
 !-------------------------------------------------------------------------------
-  USE assert_eq_m, ONLY: assert_eq
+  USE lmdz_assert_eq, ONLY: assert_eq
   IMPLICIT NONE
 !-------------------------------------------------------------------------------
@@ -128,5 +128,5 @@
 !  * interbar is TRUE for barycentric interpolation.
 !-------------------------------------------------------------------------------
-  USE assert_eq_m, ONLY: assert_eq
+  USE lmdz_assert_eq, ONLY: assert_eq
   IMPLICIT NONE
 !-------------------------------------------------------------------------------
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert.F90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert.F90	(revision 5113)
@@ -5,5 +5,5 @@
   use ioipsl, only: getin
   use new_unit_m, only: new_unit
-  use assert_m, only: assert
+  use lmdz_assert, only: assert
   USE comvert_mod, ONLY: ap, bp, aps, bps, nivsigs, nivsig, dpres, presnivs, &
                          pseudoalt, pa, preff, scaleheight, presinter
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert_noterre.f90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert_noterre.f90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert_noterre.f90	(revision 5113)
@@ -295,5 +295,5 @@
 
 
-  implicit none
+  IMPLICIT NONE
   real :: x1, x2, sig,pa,preff, newsig, F
   integer :: j
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/fxhyp_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/fxhyp_m.F90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/fxhyp_m.F90	(revision 5113)
@@ -18,5 +18,5 @@
     ! 1., taux=0., clon=0.) est à - 180 degrés.
 
-    use arth_m, only: arth
+    use lmdz_arth, only: arth
     use invert_zoom_x_m, only: invert_zoom_x, nmax
     use nrtype, only: pi, pi_d, twopi, twopi_d, k8
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gr_int_dyn.f90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gr_int_dyn.f90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gr_int_dyn.f90	(revision 5113)
@@ -3,5 +3,5 @@
 
 SUBROUTINE gr_int_dyn(champin,champdyn,iim,jp1)
-  implicit none
+  IMPLICIT NONE
   !=======================================================================
   !   passage d'un champ interpole a un champ sur grille scalaire
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inigrads.f90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inigrads.f90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inigrads.f90	(revision 5113)
@@ -5,5 +5,5 @@
         , dt, file, titlel)
 
-  implicit none
+  IMPLICIT NONE
 
   integer :: if, im, jm, lm, i, j, l
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/initdynav.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/initdynav.F90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/initdynav.F90	(revision 5113)
@@ -11,5 +11,5 @@
   USE temps_mod, ONLY: itau_dyn
   
-  implicit none
+  IMPLICIT NONE
 
 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/initfluxsto.f90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/initfluxsto.f90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/initfluxsto.f90	(revision 5113)
@@ -11,5 +11,5 @@
   USE temps_mod, ONLY: annee_ref, day_ref, itau_dyn
 
-  implicit none
+  IMPLICIT NONE
 
   !
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inithist.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inithist.F90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inithist.F90	(revision 5113)
@@ -11,5 +11,5 @@
   USE temps_mod, ONLY: itau_dyn
 
-  implicit none
+  IMPLICIT NONE
 
   !
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inter_barxy_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inter_barxy_m.F90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inter_barxy_m.F90	(revision 5113)
@@ -6,7 +6,7 @@
   ! Authors: Robert SADOURNY, Phu LE VAN, Lionel GUEZ
 
-  implicit none
-
-  private
+  IMPLICIT NONE
+
+  PRIVATE
   public inter_barxy
 
@@ -15,6 +15,6 @@
   SUBROUTINE inter_barxy(dlonid, dlatid, champ, rlonimod, rlatimod, champint)
 
-    use assert_eq_m, only: assert_eq
-    use assert_m, only: assert
+    use lmdz_assert_eq, only: assert_eq
+    use lmdz_assert, only: assert
 
     include "dimensions.h"
@@ -117,5 +117,5 @@
     !      ( Les abscisses sont exprimees en degres)
 
-    use assert_eq_m, only: assert_eq
+    use lmdz_assert_eq, only: assert_eq
 
     IMPLICIT NONE
@@ -299,5 +299,5 @@
     ! L'indice 1 correspond à l'interface maille 1 -- maille 2.
 
-    use assert_m, only: assert
+    use lmdz_assert, only: assert
 
     IMPLICIT NONE
@@ -373,5 +373,5 @@
     ! Finally, the procedure adds 90° as the last value of the array.
 
-    use assert_eq_m, only: assert_eq
+    use lmdz_assert_eq, only: assert_eq
     use comconst_mod, only: pi
 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/interpost.f90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/interpost.f90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/interpost.f90	(revision 5113)
@@ -4,5 +4,5 @@
   SUBROUTINE interpost(q,qppm)
 
-   implicit none
+   IMPLICIT NONE
 
 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/interpre.f90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/interpre.f90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/interpre.f90	(revision 5113)
@@ -9,5 +9,5 @@
   USE comvert_mod, ONLY: ap, bp
 
-   implicit none
+   IMPLICIT NONE
 
   include "dimensions.h"
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/invert_zoom_x_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/invert_zoom_x_m.F90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/invert_zoom_x_m.F90	(revision 5113)
@@ -1,5 +1,5 @@
 module invert_zoom_x_m
 
-  implicit none
+  IMPLICIT NONE
 
   INTEGER, PARAMETER:: nmax = 30000
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/iso_verif_dyn.f90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/iso_verif_dyn.f90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/iso_verif_dyn.f90	(revision 5113)
@@ -1,12 +1,12 @@
   function iso_verif_noNaN_nostop(x,err_msg)
-    implicit none
-    ! ! si x est NaN, on affiche message
-    ! ! d'erreur et return 1 si erreur
+    IMPLICIT NONE
+    ! si x est NaN, on affiche message
+    ! d'erreur et return 1 si erreur
 
-    ! ! input:
+    ! input:
     real :: x
     character(len=*) :: err_msg ! message d''erreur à afficher
 
-    ! ! output
+    ! output
     real :: borne
     parameter (borne=1e19)
@@ -27,14 +27,14 @@
   function iso_verif_egalite_nostop &
           (a,b,err_msg)
-    implicit none
-    ! ! compare a et b. Si pas egal, on affiche message
-    ! ! d'erreur et stoppe
-    ! ! pour egalite, on verifie erreur absolue et arreur relative
+    IMPLICIT NONE
+    ! compare a et b. Si pas egal, on affiche message
+    ! d'erreur et stoppe
+    ! pour egalite, on verifie erreur absolue et arreur relative
 
-    ! ! input:
+    ! input:
     real :: a, b
     character(len=*) :: err_msg ! message d''erreur à afficher
 
-    ! ! locals
+    ! locals
     real :: errmax ! erreur maximale en absolu.
     real :: errmaxrel ! erreur maximale en relatif autorisée
@@ -42,5 +42,5 @@
     parameter (errmaxrel=1e-3)
 
-    ! ! output
+    ! output
     integer :: iso_verif_egalite_nostop
 
@@ -65,12 +65,12 @@
           (x,iso,q,err_msg)
     USE infotrac, ONLY: isoName, getKey
-    implicit none
+    IMPLICIT NONE
 
-    ! ! input:
+    ! input:
     real :: x,q
     integer :: iso ! 2=HDO, 1=O18
     character(len=*) :: err_msg ! message d''erreur à afficher
 
-    ! ! locals
+    ! locals
     real :: qmin,deltaD
     real :: deltaDmax,deltaDmin,tnat
@@ -78,10 +78,10 @@
     parameter (deltaDmax=200.0,deltaDmin=-999.9)
 
-    ! ! output
+    ! output
     integer :: iso_verif_aberrant_nostop
 
     iso_verif_aberrant_nostop=0
 
-    ! ! verifier que HDO est raisonable
+    ! verifier que HDO est raisonable
      if (q>qmin) then
          IF(getKey('tnat', tnat, isoName(iso))) THEN
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/ppm3d.f90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/ppm3d.f90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/ppm3d.f90	(revision 5113)
@@ -65,5 +65,5 @@
         JNP,j1,NLAY,AP,BP,PT,AE,fill,dum,Umax)
 
-  implicit none
+  IMPLICIT NONE
 
   ! rajout de déclarations
@@ -791,5 +791,5 @@
 SUBROUTINE FZPPM(IMR,JNP,NLAY,j1,DQ,WZ,P,DC,DQDT,AR,AL,A6, &
         flux,wk1,wk2,wz2,delp,KORD)
-  implicit none
+  IMPLICIT NONE
   integer,parameter :: kmax = 150
   real,parameter :: R23 = 2./3., R3 = 1./3.
@@ -966,5 +966,5 @@
 SUBROUTINE xtp(IMR,JNP,IML,j1,j2,JN,JS,PU,DQ,Q,UC, &
         fx1,xmass,IORD)
-  implicit none
+  IMPLICIT NONE
   integer :: IMR,JNP,IML,j1,j2,JN,JS,IORD
   real :: PU,DQ,Q,UC,fx1,xmass
@@ -1088,5 +1088,5 @@
 !
 SUBROUTINE fxppm(IMR,IML,UT,P,DC,flux,IORD)
-  implicit none
+  IMPLICIT NONE
   integer :: IMR,IML,IORD
   real :: UT,P,DC,flux
@@ -1151,5 +1151,5 @@
 !
 SUBROUTINE xmist(IMR,IML,P,DC)
-  implicit none
+  IMPLICIT NONE
   integer :: IMR,IML
   real,parameter :: R24 = 1./24.
@@ -1169,5 +1169,5 @@
 SUBROUTINE ytp(IMR,JNP,j1,j2,acosp,RCAP,DQ,P,VC,DC2 &
         ,ymass,fx,A6,AR,AL,JORD)
-  implicit none
+  IMPLICIT NONE
   integer :: IMR,JNP,j1,j2,JORD
   real :: acosp,RCAP,DQ,P,VC,DC2,ymass,fx,A6,AR,AL
@@ -1239,5 +1239,5 @@
 !
 subroutine  ymist(IMR,JNP,j1,P,DC,ID)
-  implicit none
+  IMPLICIT NONE
   integer :: IMR,JNP,j1,ID
   real,parameter :: R24 = 1./24.
@@ -1321,5 +1321,5 @@
 !
 SUBROUTINE fyppm(VC,P,DC,flux,IMR,JNP,j1,j2,A6,AR,AL,JORD)
-  implicit none
+  IMPLICIT NONE
   integer :: IMR,JNP,j1,j2,JORD
   real,parameter :: R3 = 1./3., R23 = 2./3.
@@ -1401,5 +1401,5 @@
 !
   SUBROUTINE yadv(IMR,JNP,j1,j2,p,VA,ady,wk,IAD)
-    implicit none
+    IMPLICIT NONE
     integer :: IMR,JNP,j1,j2,IAD
     REAL :: p(IMR,JNP),ady(IMR,JNP),VA(IMR,JNP)
@@ -1491,5 +1491,5 @@
 !
   SUBROUTINE xadv(IMR,JNP,j1,j2,p,UA,JS,JN,IML,adx,IAD)
-    implicit none
+    IMPLICIT NONE
     INTEGER :: IMR,JNP,j1,j2,JS,JN,IML,IAD
     REAL :: p(IMR,JNP),adx(IMR,JNP),qtmp(-IMR:IMR+IMR),UA(IMR,JNP)
@@ -1584,5 +1584,5 @@
 !
 SUBROUTINE lmtppm(DC,A6,AR,AL,P,IM,LMT)
-  implicit none
+  IMPLICIT NONE
   !
   ! A6 =  CURVATURE OF THE TEST PARABOLA
@@ -1665,5 +1665,5 @@
 !
 SUBROUTINE A2C(U,V,IMR,JMR,j1,j2,CRX,CRY,dtdx5,DTDY5)
-  implicit none
+  IMPLICIT NONE
   integer :: IMR,JMR,j1,j2
   real :: U(IMR,*),V(IMR,*),CRX(IMR,*),CRY(IMR,*),DTDX5(*),DTDY5
@@ -1687,5 +1687,5 @@
 !
 SUBROUTINE cosa(cosp,cose,JNP,PI,DP)
-  implicit none
+  IMPLICIT NONE
   integer :: JNP
   real :: cosp(*),cose(*),PI,DP
@@ -1720,5 +1720,5 @@
 !
 SUBROUTINE cosc(cosp,cose,JNP,PI,DP)
-  implicit none
+  IMPLICIT NONE
   integer :: JNP
   real :: cosp(*),cose(*),PI,DP
@@ -1853,5 +1853,5 @@
 !
 SUBROUTINE filcr(q,IMR,JNP,j1,j2,cosp,acosp,icr,tiny)
-  implicit none
+  IMPLICIT NONE
   integer :: IMR,JNP,j1,j2,icr
   real :: q(IMR,*),cosp(*),acosp(*),tiny
@@ -1956,5 +1956,5 @@
 !
 SUBROUTINE filns(q,IMR,JNP,j1,j2,cosp,acosp,ipy,tiny)
-  implicit none
+  IMPLICIT NONE
   integer :: IMR,JNP,j1,j2,ipy
   real :: q(IMR,*),cosp(*),acosp(*),tiny
@@ -2043,5 +2043,5 @@
 !
 SUBROUTINE filew(q,qtmp,IMR,JNP,j1,j2,ipx,tiny)
-  implicit none
+  IMPLICIT NONE
   integer :: IMR,JNP,j1,j2,ipx
   real :: q(IMR,*),qtmp(JNP,IMR),tiny
@@ -2125,5 +2125,5 @@
 !
 SUBROUTINE zflip(q,im,km,nc)
-  implicit none
+  IMPLICIT NONE
   ! This routine flip the array q (in the vertical).
   integer :: im,km,nc
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/principal_cshift_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/principal_cshift_m.F90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/principal_cshift_m.F90	(revision 5113)
@@ -1,5 +1,5 @@
 module principal_cshift_m
 
-  implicit none
+  IMPLICIT NONE
 
 contains
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/sortvarc.f90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/sortvarc.f90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/sortvarc.f90	(revision 5113)
@@ -160,6 +160,6 @@
   END IF
 
-  ! ! compute relative changes in etot,... (except if 'reference' values
-  ! ! are zero, which can happen when using iniacademic)
+  ! compute relative changes in etot,... (except if 'reference' values
+  ! are zero, which can happen when using iniacademic)
   if (etot0/=0) then
     etot= etot/etot0
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/traceurpole.f90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/traceurpole.f90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/traceurpole.f90	(revision 5113)
@@ -4,5 +4,5 @@
     SUBROUTINE traceurpole(q,masse)
 
-      implicit none
+      IMPLICIT NONE
 
   include "dimensions.h"
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/ugeostr.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/ugeostr.F90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/ugeostr.F90	(revision 5113)
@@ -13,5 +13,5 @@
   use comconst_mod, only: omeg, rad
   
-  implicit none
+  IMPLICIT NONE
 
   include "dimensions.h"
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/writedynav.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/writedynav.F90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/writedynav.F90	(revision 5113)
@@ -9,5 +9,5 @@
   USE temps_mod, ONLY: itau_dyn
 
-  implicit none
+  IMPLICIT NONE
 
   !   Ecriture du fichier histoire au format IOIPSL
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/writehist.f90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/writehist.f90	(revision 5108)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/writehist.f90	(revision 5113)
@@ -9,5 +9,5 @@
   USE temps_mod, ONLY: itau_dyn
 
-  implicit none
+  IMPLICIT NONE
 
   !