Index: /LMDZ6/trunk/libf/dyn3d/addfi.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/addfi.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/addfi.f90	(revision 5281)
@@ -11,4 +11,5 @@
   USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
+  USE comgeom_mod_h
   IMPLICIT NONE
   !
@@ -45,9 +46,4 @@
   !
   !-----------------------------------------------------------------------
-  !
-  !    0.  Declarations :
-  !    ------------------
-  !
-  include "comgeom.h"
   !
   !    Arguments :
Index: /LMDZ6/trunk/libf/dyn3d/advect.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/advect.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/advect.f90	(revision 5281)
@@ -3,4 +3,5 @@
 !
 SUBROUTINE advect(ucov,vcov,teta,w,massebx,masseby,du,dv,dteta)
+  USE comgeom_mod_h
   USE comconst_mod, ONLY: daysec
   USE logic_mod, ONLY: conser
@@ -28,5 +29,4 @@
   !   Declarations:
   !   -------------
-  include "comgeom.h"
 
   !   Arguments:
Index: /LMDZ6/trunk/libf/dyn3d/advtrac.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/advtrac.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/advtrac.f90	(revision 5281)
@@ -9,4 +9,5 @@
    !            M.A Filiberti (04/2002)
    !
+   USE comgeom2_mod_h
    USE comdissip_mod_h
    USE infotrac,     ONLY: nqtot, tracers, isoCheck
@@ -22,5 +23,4 @@
 
 
-   include "comgeom2.h"
    include "description.h"
    include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3d/bilan_dyn.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/bilan_dyn.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/bilan_dyn.f90	(revision 5281)
@@ -10,4 +10,5 @@
   !             vQ..A=Cp T + L * ...
 
+  USE comgeom2_mod_h
   USE IOIPSL
   USE comconst_mod, ONLY: pi, cpp
@@ -22,5 +23,4 @@
 
 
-  include "comgeom2.h"
   include "iniprint.h"
 
Index: /LMDZ6/trunk/libf/dyn3d/caldyn.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/caldyn.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/caldyn.f90	(revision 5281)
@@ -7,4 +7,5 @@
 
 
+  USE comgeom_mod_h
   USE comvert_mod, ONLY: ap, bp
 
@@ -32,5 +33,4 @@
 
 
-  include "comgeom.h"
 
   !   Arguments:
Index: /LMDZ6/trunk/libf/dyn3d/covnat.F90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/covnat.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/covnat.F90	(revision 5281)
@@ -3,4 +3,5 @@
 !
 SUBROUTINE covnat (klevel,ucov, vcov, unat, vnat )
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
   USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -21,5 +22,4 @@
   !
   !=======================================================================
-  INCLUDE "comgeom.h"
 
   INTEGER :: klevel
Index: /LMDZ6/trunk/libf/dyn3d/dissip.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/dissip.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/dissip.f90	(revision 5281)
@@ -4,4 +4,5 @@
 SUBROUTINE dissip( vcov,ucov,teta,p, dv,du,dh )
   !
+  USE comgeom_mod_h
   USE comdissipn_mod_h
   USE comdissnew_mod_h
@@ -34,5 +35,4 @@
 
 
-  include "comgeom.h"
 
   !   Arguments:
Index: /LMDZ6/trunk/libf/dyn3d/dynetat0.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/dynetat0.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/dynetat0.f90	(revision 5281)
@@ -6,4 +6,5 @@
 ! Purpose: Initial state reading.
 !-------------------------------------------------------------------------------
+  USE comgeom2_mod_h
   USE infotrac,    ONLY: nqtot, tracers, niso, iqIsoPha, iH2O, isoName, &
                          new2oldH2O, newHNO3, oldHNO3
@@ -30,5 +31,4 @@
 
 
-  include "comgeom2.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3d/dynredem.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/dynredem.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/dynredem.f90	(revision 5281)
@@ -18,12 +18,12 @@
   USE temps_mod, ONLY: annee_ref, day_ref, itau_dyn, itaufin, start_time
   USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0
-  
+
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
-USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
+  USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
+  USE comgeom2_mod_h
 IMPLICIT NONE
 
 
-  include "comgeom2.h"
   include "description.h"
   include "iniprint.h"
@@ -31,5 +31,5 @@
 ! Arguments:
   CHARACTER(LEN=*), INTENT(IN) :: fichnom          !--- FILE NAME
-  INTEGER,          INTENT(IN) :: iday_end         !--- 
+  INTEGER,          INTENT(IN) :: iday_end         !---
   REAL,             INTENT(IN) :: phis(iip1, jjp1) !--- GROUND GEOPOTENTIAL
 !===============================================================================
@@ -72,5 +72,5 @@
   tab_cntrl(19) = preff
 
-!    .....    parameters for zoom    ......   
+!    .....    parameters for zoom    ......
   tab_cntrl(20) = clon
   tab_cntrl(21) = clat
@@ -169,13 +169,13 @@
                           err, modname, fil, msg
   USE temps_mod, ONLY: itau_dyn, itaufin
-  
+
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
-USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
+  USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
+  USE comgeom_mod_h
 IMPLICIT NONE
 
 
   include "description.h"
-  include "comgeom.h"
   include "iniprint.h"
 !===============================================================================
Index: /LMDZ6/trunk/libf/dyn3d/fluxstokenc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/fluxstokenc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/fluxstokenc.f90	(revision 5281)
@@ -6,5 +6,6 @@
   ! This routine is designed to work with ioipsl
 
-   USE IOIPSL
+   USE comgeom_mod_h
+  USE IOIPSL
   !
   ! Auteur :  F. Hourdin
@@ -20,5 +21,4 @@
 
 
-  include "comgeom.h"
   include "tracstoke.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3d/friction.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/friction.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/friction.f90	(revision 5281)
@@ -4,4 +4,5 @@
 !=======================================================================
 SUBROUTINE friction(ucov,vcov,pdt)
+  USE comgeom2_mod_h
   USE control_mod
   USE IOIPSL
@@ -25,5 +26,4 @@
 
 
-  include "comgeom2.h"
   include "iniprint.h"
   include "academic.h"
Index: /LMDZ6/trunk/libf/dyn3d/gcm.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/gcm.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/gcm.f90	(revision 5281)
@@ -5,4 +5,5 @@
 !
 PROGRAM gcm
+  USE comgeom_mod_h
   USE comdissnew_mod_h
   USE IOIPSL
@@ -57,5 +58,4 @@
   !   Declarations:
   !   -------------
-  include "comgeom.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3d/groupe.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/groupe.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/groupe.f90	(revision 5281)
@@ -4,4 +4,5 @@
 subroutine groupe(pext,pbaru,pbarv,pbarum,pbarvm,wm)
 
+  USE comgeom2_mod_h
   use comconst_mod, only: ngroup
 
@@ -24,5 +25,4 @@
 
 
-  include "comgeom2.h"
 
   ! integer ngroup
Index: /LMDZ6/trunk/libf/dyn3d/groupeun.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/groupeun.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/groupeun.f90	(revision 5281)
@@ -4,4 +4,5 @@
 SUBROUTINE groupeun(jjmax,llmax,q)
 
+  USE comgeom2_mod_h
   USE comconst_mod, ONLY: ngroup
 
@@ -13,5 +14,4 @@
 
 
-  include "comgeom2.h"
 
   INTEGER :: jjmax,llmax
@@ -141,4 +141,5 @@
 SUBROUTINE INIT_GROUPEUN(airen_tab, aires_tab)
 
+  USE comgeom2_mod_h
   USE comconst_mod, ONLY: ngroup
 
@@ -150,5 +151,4 @@
 
 
-  include "comgeom2.h"
 
   ! INTEGER ngroup
Index: /LMDZ6/trunk/libf/dyn3d/guide_mod.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/guide_mod.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/guide_mod.f90	(revision 5281)
@@ -8,6 +8,5 @@
 !            F. Codron 01/09
 !=======================================================================
-
-  USE getparam, only: ini_getparam, fin_getparam, getpar
+    USE getparam, only: ini_getparam, fin_getparam, getpar
   USE Write_Field
   use netcdf, only: nf90_nowrite, nf90_open, nf90_inq_varid, nf90_close, &
@@ -92,5 +91,5 @@
 
 ! ---------------------------------------------
-! Lecture des parametres:  
+! Lecture des parametres:
 ! ---------------------------------------------
     call ini_getparam("nudging_parameters_out.txt")
@@ -172,5 +171,5 @@
 
     call fin_getparam
-    
+
 ! ---------------------------------------------
 ! Determination du nombre de niveaux verticaux
@@ -230,5 +229,5 @@
 
 
-    endif 
+    endif
     error=nf90_inq_dimid(ncidpl,'LEVEL',rid)
     IF (error.NE.NF90_NOERR) error=nf90_inq_dimid(ncidpl,'PRESSURE',rid)
@@ -237,5 +236,5 @@
     ENDIF
     error=nf90_inquire_dimension(ncidpl,rid,len=nlevnc)
-    write(*,*)trim(modname)//' : number of vertical levels nlevnc', nlevnc 
+    write(*,*)trim(modname)//' : number of vertical levels nlevnc', nlevnc
     rcod = nf90_close(ncidpl)
 
@@ -264,5 +263,5 @@
     IF (error /= 0) CALL abort_gcm(modname,abort_message,1)
     alpha_u=0.;alpha_v=0;alpha_T=0;alpha_Q=0;alpha_P=0
-    
+
     IF (guide_u) THEN
         ALLOCATE(unat1(iip1,jjp1,nlevnc), stat = error)
@@ -288,5 +287,5 @@
         tnat1=0.;tnat2=0.;tgui1=0.;tgui2=0.
     ENDIF
-     
+
     IF (guide_Q) THEN
         ALLOCATE(qnat1(iip1,jjp1,nlevnc), stat = error)
@@ -361,5 +360,5 @@
     USE comconst_mod, ONLY: cpp, dtvr, daysec,kappa
     USE comvert_mod, ONLY: ap, bp, preff, presnivs, pressure_exner
- 
+
     USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -391,5 +390,5 @@
     REAL          :: tau,reste ! position entre 2 etats de guidage
     REAL, SAVE    :: factt ! pas de temps en fraction de jour
-    
+
     INTEGER       :: l
     CHARACTER(LEN=20) :: modname="guide_main"
@@ -402,10 +401,10 @@
     IF (first) THEN
         first=.FALSE.
-        CALL guide_init 
+        CALL guide_init
         itau_test=1001
         step_rea=1
         count_no_rea=0
 ! Calcul des constantes de rappel
-        factt=dtvr*iperiod/daysec 
+        factt=dtvr*iperiod/daysec
         call tau2alpha(3,iip1,jjm ,factt,tau_min_v,tau_max_v,alpha_v)
         call tau2alpha(2,iip1,jjp1,factt,tau_min_u,tau_max_u,alpha_u)
@@ -421,5 +420,5 @@
             enddo
         endif
-! ini_anal: etat initial egal au guidage        
+! ini_anal: etat initial egal au guidage
         IF (ini_anal) THEN
             CALL guide_interp(ps,teta)
@@ -464,5 +463,5 @@
             itau
               abort_message='stopped'
-              CALL abort_gcm(modname,abort_message,1) 
+              CALL abort_gcm(modname,abort_message,1)
           ELSE
               IF (guide_v) vnat1=vnat2
@@ -503,8 +502,8 @@
 
 !-----------------------------------------------------------------------
-!   Ajout des champs de guidage 
+!   Ajout des champs de guidage
 !-----------------------------------------------------------------------
 ! Sauvegarde du guidage?
-    f_out=((MOD(itau,iguide_sav).EQ.0).AND.guide_sav)  
+    f_out=((MOD(itau,iguide_sav).EQ.0).AND.guide_sav)
     IF (f_out) THEN
       ! compute pressures at layer interfaces
@@ -522,5 +521,5 @@
       CALL guide_out("SP",jjp1,llm,p(:,1:llm))
     ENDIF
-    
+
     if (guide_u) then
         if (guide_add) then
@@ -528,5 +527,5 @@
         else
            f_add=(1.-tau)*ugui1+tau*ugui2-ucov
-        endif 
+        endif
         if (guide_zon) CALL guide_zonave(1,jjp1,llm,f_add)
         CALL guide_addfield(ip1jmp1,llm,f_add,alpha_u)
@@ -542,5 +541,5 @@
         else
            f_add=(1.-tau)*tgui1+tau*tgui2-teta
-        endif 
+        endif
         if (guide_zon) CALL guide_zonave(2,jjp1,llm,f_add)
         CALL guide_addfield(ip1jmp1,llm,f_add,alpha_T)
@@ -554,5 +553,5 @@
         else
            f_add(1:ip1jmp1,1)=(1.-tau)*psgui1+tau*psgui2-ps
-        endif 
+        endif
         if (guide_zon) CALL guide_zonave(2,jjp1,1,f_add(1:ip1jmp1,1))
         CALL guide_addfield(ip1jmp1,1,f_add(1:ip1jmp1,1),alpha_P)
@@ -568,5 +567,5 @@
         else
            f_add=(1.-tau)*qgui1+tau*qgui2-q
-        endif 
+        endif
         if (guide_zon) CALL guide_zonave(2,jjp1,llm,f_add)
         CALL guide_addfield(ip1jmp1,llm,f_add,alpha_Q)
@@ -580,5 +579,5 @@
         else
            f_add(1:ip1jm,:)=(1.-tau)*vgui1+tau*vgui2-vcov
-        endif 
+        endif
         if (guide_zon) CALL guide_zonave(2,jjm,llm,f_add(1:ip1jm,:))
         CALL guide_addfield(ip1jm,llm,f_add(1:ip1jm,:),alpha_v)
@@ -599,5 +598,5 @@
     INTEGER,                      INTENT(IN)    :: hsize
     INTEGER,                      INTENT(IN)    :: vsize
-    REAL, DIMENSION(hsize),       INTENT(IN)    :: alpha 
+    REAL, DIMENSION(hsize),       INTENT(IN)    :: alpha
     REAL, DIMENSION(hsize,vsize), INTENT(INOUT) :: field
 
@@ -615,14 +614,14 @@
 
     USE comconst_mod, ONLY: pi
-    
+
     USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
+    USE comgeom_mod_h
 IMPLICIT NONE
 
 
 
-    INCLUDE "comgeom.h"
-    
+
     ! input/output variables
     INTEGER,                           INTENT(IN)    :: typ
@@ -665,5 +664,5 @@
                 fieldm(j,l)=fieldm(j,l)+field(ij,l)
             ENDDO
-        ENDDO 
+        ENDDO
         fieldm(:,l)=fieldm(:,l)/REAL(imax(typ)-imin(typ)+1)
     ! Compute forcing
@@ -680,5 +679,5 @@
 !=======================================================================
   SUBROUTINE guide_interp(psi,teta)
-  
+
   use exner_hyb_m, only: exner_hyb
   use exner_milieu_m, only: exner_milieu
@@ -686,12 +685,8 @@
   use comvert_mod, only: preff, pressure_exner, bp, ap
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
-USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
+  USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
+  USE comgeom2_mod_h
 IMPLICIT NONE
-
-
-
-  include "comgeom2.h"
-
   REAL, DIMENSION (iip1,jjp1),     INTENT(IN) :: psi ! Psol gcm
   REAL, DIMENSION (iip1,jjp1,llm), INTENT(IN) :: teta ! Temp. Pot. gcm
@@ -702,21 +697,21 @@
   REAL, DIMENSION (iip1,jjp1,llm)    :: plunc,plsnc !niveaux pression modele
   REAL, DIMENSION (iip1,jjm,llm)     :: plvnc       !niveaux pression modele
-  REAL, DIMENSION (iip1,jjp1,llmp1)  :: p           ! pression intercouches 
+  REAL, DIMENSION (iip1,jjp1,llmp1)  :: p           ! pression intercouches
   REAL, DIMENSION (iip1,jjp1,llm)    :: pls, pext   ! var intermediaire
-  REAL, DIMENSION (iip1,jjp1,llm)    :: pbarx 
-  REAL, DIMENSION (iip1,jjm,llm)     :: pbary 
+  REAL, DIMENSION (iip1,jjp1,llm)    :: pbarx
+  REAL, DIMENSION (iip1,jjm,llm)     :: pbary
   ! Variables pour fonction Exner (P milieu couche)
   REAL, DIMENSION (iip1,jjp1,llm)    :: pk
-  REAL, DIMENSION (iip1,jjp1)        :: pks    
+  REAL, DIMENSION (iip1,jjp1)        :: pks
   REAL                               :: prefkap,unskap
   ! Pression de vapeur saturante
   REAL, DIMENSION (ip1jmp1,llm)      :: qsat
   !Variables intermediaires interpolation
-  REAL, DIMENSION (iip1,jjp1,llm)    :: zu1,zu2 
+  REAL, DIMENSION (iip1,jjp1,llm)    :: zu1,zu2
   REAL, DIMENSION (iip1,jjm,llm)     :: zv1,zv2
-  
+
   INTEGER                            :: i,j,l,ij
   CHARACTER(LEN=20),PARAMETER :: modname="guide_interp"
-  
+
     write(*,*)trim(modname)//': interpolate nudging variables'
 ! -----------------------------------------------------------------
@@ -767,7 +762,7 @@
         endif
     endif
-    
+
 ! -----------------------------------------------------------------
-! Calcul niveaux pression modele 
+! Calcul niveaux pression modele
 ! -----------------------------------------------------------------
     CALL pression( ip1jmp1, ap, bp, psi, p )
@@ -839,6 +834,6 @@
                     ugui2(ij,l)=zu2(i,j,l)*cu(i,j)
                 enddo
-                ugui1(j*iip1,l)=ugui1((j-1)*iip1+1,l)    
-                ugui2(j*iip1,l)=ugui2((j-1)*iip1+1,l)    
+                ugui1(j*iip1,l)=ugui1((j-1)*iip1+1,l)
+                ugui2(j*iip1,l)=ugui2((j-1)*iip1+1,l)
             enddo
             do i=1,iip1
@@ -850,5 +845,5 @@
         enddo
     ENDIF
-    
+
     IF (guide_T) THEN
         CALL pres2lev(tnat1,zu1,nlevnc,llm,plnc1,plsnc,iip1,jjp1,invert_p)
@@ -869,12 +864,12 @@
                     enddo
                 ENDIF
-                tgui1(j*iip1,l)=tgui1((j-1)*iip1+1,l)    
-                tgui2(j*iip1,l)=tgui2((j-1)*iip1+1,l)    
+                tgui1(j*iip1,l)=tgui1((j-1)*iip1+1,l)
+                tgui2(j*iip1,l)=tgui2((j-1)*iip1+1,l)
             enddo
             do i=1,iip1
                 tgui1(i,l)=tgui1(1,l)
-                tgui1(ip1jm+i,l)=tgui1(ip1jm+1,l) 
+                tgui1(ip1jm+i,l)=tgui1(ip1jm+1,l)
                 tgui2(i,l)=tgui2(1,l)
-                tgui2(ip1jm+i,l)=tgui2(ip1jm+1,l) 
+                tgui2(ip1jm+i,l)=tgui2(ip1jm+1,l)
             enddo
         enddo
@@ -893,10 +888,10 @@
                     vgui2(ij,l)=zv2(i,j,l)*cv(i,j)
                 enddo
-                vgui1(j*iip1,l)=vgui1((j-1)*iip1+1,l)    
-                vgui2(j*iip1,l)=vgui2((j-1)*iip1+1,l)    
+                vgui1(j*iip1,l)=vgui1((j-1)*iip1+1,l)
+                vgui2(j*iip1,l)=vgui2((j-1)*iip1+1,l)
             enddo
         enddo
     ENDIF
-    
+
     IF (guide_Q) THEN
         ! On suppose qu'on a la bonne variable dans le fichier de guidage:
@@ -911,12 +906,12 @@
                     qgui2(ij,l)=zu2(i,j,l)
                 enddo
-                qgui1(j*iip1,l)=qgui1((j-1)*iip1+1,l)    
-                qgui2(j*iip1,l)=qgui2((j-1)*iip1+1,l)    
+                qgui1(j*iip1,l)=qgui1((j-1)*iip1+1,l)
+                qgui2(j*iip1,l)=qgui2((j-1)*iip1+1,l)
             enddo
             do i=1,iip1
                 qgui1(i,l)=qgui1(1,l)
-                qgui1(ip1jm+i,l)=qgui1(ip1jm+1,l) 
+                qgui1(ip1jm+i,l)=qgui1(ip1jm+1,l)
                 qgui2(i,l)=qgui2(1,l)
-                qgui2(ip1jm+i,l)=qgui2(ip1jm+1,l) 
+                qgui2(ip1jm+i,l)=qgui2(ip1jm+1,l)
             enddo
         enddo
@@ -924,5 +919,5 @@
             CALL q_sat(iip1*jjp1*llm,teta*pk/cpp,plsnc,qsat)
             qgui1=qgui1*qsat*0.01 !hum. rel. en %
-            qgui2=qgui2*qsat*0.01 
+            qgui2=qgui2*qsat*0.01
         ENDIF
     ENDIF
@@ -937,14 +932,10 @@
     use comconst_mod, only: pi
     use serre_mod, only: clon, clat, grossismx, grossismy
-    
+
     USE dimensions_mod, ONLY: iim, jjm, llm, ndm
-USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
+    USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
-implicit none
-
-
-
-    include "comgeom2.h"
-
+    USE comgeom2_mod_h
+    implicit none
 ! input arguments :
     INTEGER, INTENT(IN) :: typ    ! u(2),v(3), ou scalaire(1)
@@ -953,6 +944,6 @@
     REAL, INTENT(IN)    :: taumin,taumax
 ! output arguments:
-    REAL, DIMENSION(pim,pjm), INTENT(OUT) :: alpha 
-  
+    REAL, DIMENSION(pim,pjm), INTENT(OUT) :: alpha
+
 !  local variables:
     LOGICAL, SAVE               :: first=.TRUE.
@@ -1040,8 +1031,8 @@
         ENDIF
 ! Premier appel: calcul des aires min et max et de gamma.
-        IF (first) THEN 
+        IF (first) THEN
             first=.FALSE.
             ! coordonnees du centre du zoom
-            CALL coordij(clon,clat,ilon,ilat) 
+            CALL coordij(clon,clat,ilon,ilat)
             ! aire de la maille au centre du zoom
             dxdy_min=dxdys(ilon,ilat)
@@ -1067,5 +1058,5 @@
               endif
               gamma=log(0.5)/log(gamma)
-              if (gamma4) then 
+              if (gamma4) then
                 gamma=min(gamma,4.)
               endif
@@ -1168,5 +1159,5 @@
              ENDIF
              rcode = nf90_inq_varid(ncidp, 'PRES', varidp)
-             IF (rcode.NE.NF90_NOERR) THEN 
+             IF (rcode.NE.NF90_NOERR) THEN
               abort_message='Nudging: error -> no PRES variable in file P.nc'
               CALL abort_gcm(modname,abort_message,1)
@@ -1221,8 +1212,8 @@
              write(*,*) trim(modname)//' ncidv,varidv',ncidv,varidv
              if (ncidpl.eq.-99) ncidpl=ncidv
-             
+
              status=NF90_INQ_DIMID(ncidv, "LONV", dimid)
              status=NF90_INQUIRE_DIMENSION(ncidv,dimid,namedim,lendim)
-             
+
                 IF (lendim .NE. iip1) THEN
                 abort_message='dimension LONV different from iip1 in v.nc'
@@ -1237,5 +1228,5 @@
                 CALL abort_gcm(modname,abort_message,1)
              ENDIF
-        
+
          endif
 
@@ -1350,5 +1341,5 @@
      count(4)=1
 
-! Pression 
+! Pression
      if (guide_plevs.EQ.2) then
          status=NF90_GET_VAR(ncidp,varidp,pnat2,start,count)
@@ -1382,5 +1373,5 @@
            CALL invert_lat(iip1,jjp1,nlevnc,qnat2)
          ENDIF
-         
+
      endif
 
@@ -1672,5 +1663,5 @@
 
   END SUBROUTINE guide_read2D
-  
+
 !=======================================================================
   SUBROUTINE guide_out(varname,hsize,vsize,field)
@@ -1680,13 +1671,13 @@
     use netcdf95, only: nf95_def_var, nf95_put_var
     use netcdf, only: nf90_float, nf90_def_var, nf90_put_var
-    
+
     USE dimensions_mod, ONLY: iim, jjm, llm, ndm
-USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
+    USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
+    USE comgeom2_mod_h
 IMPLICIT NONE
 
 
 
-    INCLUDE "comgeom2.h"
     
     ! Variables entree
Index: /LMDZ6/trunk/libf/dyn3d/iniacademic.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/iniacademic.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/iniacademic.f90	(revision 5281)
@@ -4,4 +4,5 @@
 SUBROUTINE iniacademic(vcov,ucov,teta,q,masse,ps,phis,time_0)
 
+  USE comgeom_mod_h
   USE filtreg_mod, ONLY: inifilr
   USE infotrac,    ONLY: nqtot, niso, iqIsoPha, tracers, isoName, addPhase
@@ -36,5 +37,4 @@
 
 
-  include "comgeom.h"
   include "academic.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3d/integrd.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/integrd.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/integrd.f90	(revision 5281)
@@ -7,4 +7,5 @@
         )
 
+  USE comgeom_mod_h
   use control_mod, only : planet_type
   use comconst_mod, only: pi
@@ -36,5 +37,4 @@
 
 
-  include "comgeom.h"
   include "iniprint.h"
 
Index: /LMDZ6/trunk/libf/dyn3d/leapfrog.F90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/leapfrog.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/leapfrog.F90	(revision 5281)
@@ -6,4 +6,5 @@
 SUBROUTINE leapfrog(ucov,vcov,teta,ps,masse,phis,q,time_0)
   !IM : pour sortir les param. du modele dans un fis. netcdf 110106
+  USE comgeom_mod_h
   USE comdissnew_mod_h
   use IOIPSL
@@ -60,5 +61,4 @@
   !   Declarations:
   !   -------------
-  include "comgeom.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3d/replay3d.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/replay3d.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/replay3d.f90	(revision 5281)
@@ -5,4 +5,5 @@
 !
 PROGRAM replay3d
+  USE comgeom2_mod_h
   USE comdissnew_mod_h
   USE comvert_mod, only :  preff, pa
@@ -57,5 +58,4 @@
   !   Declarations:
   !   -------------
-  include "comgeom2.h"
 
   REAL zdtvr
Index: /LMDZ6/trunk/libf/dyn3d/sw_case_williamson91_6.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/sw_case_williamson91_6.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/sw_case_williamson91_6.f90	(revision 5281)
@@ -26,4 +26,5 @@
   !
   !=======================================================================
+  USE comgeom_mod_h
   USE comconst_mod, ONLY: cpp, omeg, rad
   USE comvert_mod, ONLY: ap, bp, preff
@@ -39,5 +40,4 @@
 
 
-  include "comgeom.h"
   include "iniprint.h"
 
Index: /LMDZ6/trunk/libf/dyn3d/top_bound.F90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/top_bound.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/top_bound.F90	(revision 5281)
@@ -4,4 +4,5 @@
 SUBROUTINE top_bound(vcov,ucov,teta,masse,dt)
 
+  USE comgeom2_mod_h
   USE comdissipn_mod_h
 USE comconst_mod, ONLY: iflag_top_bound, mode_top_bound, &
@@ -16,5 +17,4 @@
 
 
-  include "comgeom2.h"
 
 
Index: /LMDZ6/trunk/libf/dyn3d/vlsplt.F90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/vlsplt.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/vlsplt.F90	(revision 5281)
@@ -487,4 +487,5 @@
 END SUBROUTINE vlx
 RECURSIVE SUBROUTINE vly(q,pente_max,masse,masse_adv_v,iq)
+  USE comgeom_mod_h
   USE infotrac, ONLY : nqtot,tracers, & ! CRisi
         min_qParent,min_qMass,min_ratio ! MVals et CRisi
@@ -508,5 +509,4 @@
 
 
-  include "comgeom.h"
   !
   !
Index: /LMDZ6/trunk/libf/dyn3d/vlspltqs.F90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d/vlspltqs.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d/vlspltqs.F90	(revision 5281)
@@ -545,4 +545,5 @@
 END SUBROUTINE vlxqs
 SUBROUTINE vlyqs(q,pente_max,masse,masse_adv_v,qsat,iq)
+  USE comgeom_mod_h
   USE infotrac, ONLY : nqtot,tracers ! CRisi
   !
@@ -567,5 +568,4 @@
 
 
-  include "comgeom.h"
   !
   !
Index: /LMDZ6/trunk/libf/dyn3d_common/adaptdt.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/adaptdt.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/adaptdt.f90	(revision 5281)
@@ -5,4 +5,5 @@
         masse)
 
+  USE comgeom2_mod_h
   USE comdissip_mod_h
   USE comconst_mod, ONLY: dtvr
@@ -14,5 +15,4 @@
 
 
-  include "comgeom2.h"
   include "description.h"
 
Index: /LMDZ6/trunk/libf/dyn3d_common/advn.F90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/advn.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/advn.F90	(revision 5281)
@@ -15,4 +15,5 @@
   !
   !   --------------------------------------------------------------------
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -22,5 +23,4 @@
 
 
-  include "comgeom.h"
   include "iniprint.h"
 
@@ -759,4 +759,5 @@
   !
   !   --------------------------------------------------------------------
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -766,5 +767,4 @@
 
 
-  INCLUDE "comgeom.h"
   INCLUDE "iniprint.h"
   !
@@ -893,4 +893,5 @@
   !
   !   --------------------------------------------------------------------
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -900,5 +901,4 @@
 
 
-  INCLUDE "comgeom.h"
   INCLUDE "iniprint.h"
   !
Index: /LMDZ6/trunk/libf/dyn3d_common/advy.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/advy.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/advy.f90	(revision 5281)
@@ -3,4 +3,5 @@
 !
 SUBROUTINE advy(limit,dty,pbarv,sm,s0,sx,sy,sz)
+  USE comgeom2_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -31,5 +32,4 @@
 
 
-  include "comgeom2.h"
 
   !  Arguments :
Index: /LMDZ6/trunk/libf/dyn3d_common/advyp.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/advyp.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/advyp.f90	(revision 5281)
@@ -4,4 +4,5 @@
 SUBROUTINE ADVYP(LIMIT,DTY,PBARV,SM,S0,SSX,SY,SZ &
         ,SSXX,SSXY,SSXZ,SYY,SYZ,SZZ,ntra )
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -33,5 +34,4 @@
 
 
-  include "comgeom.h"
 
   !  Arguments :
Index: /LMDZ6/trunk/libf/dyn3d_common/advzp.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/advzp.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/advzp.f90	(revision 5281)
@@ -5,4 +5,5 @@
         ,SSXX,SSXY,SSXZ,SYY,SYZ,SZZ,ntra )
 
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -36,5 +37,4 @@
 
 
-  include "comgeom.h"
   !
   !  Arguments :
Index: /LMDZ6/trunk/libf/dyn3d_common/caldyn0.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/caldyn0.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/caldyn0.f90	(revision 5281)
@@ -6,5 +6,6 @@
 ! Purpose: Compute dynamic tendencies.
 !-------------------------------------------------------------------------------
-  USE control_mod, ONLY: resetvarc 
+  USE comgeom_mod_h
+  USE control_mod, ONLY: resetvarc
   USE comvert_mod, ONLY: ap, bp
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -14,5 +15,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: DZ6/trunk/libf/dyn3d_common/comgeom.h
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/comgeom.h	(revision 5280)
+++ 	(revision )
@@ -1,33 +1,0 @@
-!
-! $Header$
-!
-!CDK comgeom
-      COMMON/comgeom/                                                   &
-     & cu(ip1jmp1),cv(ip1jm),unscu2(ip1jmp1),unscv2(ip1jm),             &
-     & aire(ip1jmp1),airesurg(ip1jmp1),aireu(ip1jmp1),                  &
-     & airev(ip1jm),unsaire(ip1jmp1),apoln,apols,                       &
-     & unsairez(ip1jm),airuscv2(ip1jm),airvscu2(ip1jm),                 &
-     & aireij1(ip1jmp1),aireij2(ip1jmp1),aireij3(ip1jmp1),              &
-     & aireij4(ip1jmp1),alpha1(ip1jmp1),alpha2(ip1jmp1),                &
-     & alpha3(ip1jmp1),alpha4(ip1jmp1),alpha1p2(ip1jmp1),               &
-     & alpha1p4(ip1jmp1),alpha2p3(ip1jmp1),alpha3p4(ip1jmp1),           &
-     & fext(ip1jm),constang(ip1jmp1),rlatu(jjp1),rlatv(jjm),            &
-     & rlonu(iip1),rlonv(iip1),cuvsurcv(ip1jm),cvsurcuv(ip1jm),         &
-     & cvusurcu(ip1jmp1),cusurcvu(ip1jmp1),cuvscvgam1(ip1jm),           &
-     & cuvscvgam2(ip1jm),cvuscugam1(ip1jmp1),                           &
-     & cvuscugam2(ip1jmp1),cvscuvgam(ip1jm),cuscvugam(ip1jmp1),         &
-     & unsapolnga1,unsapolnga2,unsapolsga1,unsapolsga2,                 &
-     & unsair_gam1(ip1jmp1),unsair_gam2(ip1jmp1),unsairz_gam(ip1jm),    &
-     & aivscu2gam(ip1jm),aiuscv2gam(ip1jm),xprimu(iip1),xprimv(iip1)
-
-!
-        REAL                                                            &
-     & cu,cv,unscu2,unscv2,aire,airesurg,aireu,airev,unsaire,apoln     ,&
-     & apols,unsairez,airuscv2,airvscu2,aireij1,aireij2,aireij3,aireij4,&
-     & alpha1,alpha2,alpha3,alpha4,alpha1p2,alpha1p4,alpha2p3,alpha3p4 ,&
-     & fext,constang,rlatu,rlatv,rlonu,rlonv,cuvscvgam1,cuvscvgam2     ,&
-     & cvuscugam1,cvuscugam2,cvscuvgam,cuscvugam,unsapolnga1,unsapolnga2&
-     & ,unsapolsga1,unsapolsga2,unsair_gam1,unsair_gam2,unsairz_gam    ,&
-     & aivscu2gam ,aiuscv2gam,cuvsurcv,cvsurcuv,cvusurcu,cusurcvu,xprimu&
-     & , xprimv
-!
Index: DZ6/trunk/libf/dyn3d_common/comgeom2.h
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/comgeom2.h	(revision 5280)
+++ 	(revision )
@@ -1,33 +1,0 @@
-!
-! $Header$
-!
-!CDK comgeom2
-      COMMON/comgeom/                                                   &
-     & cu(iip1,jjp1),cv(iip1,jjm),unscu2(iip1,jjp1),unscv2(iip1,jjm)  , &
-     & aire(iip1,jjp1),airesurg(iip1,jjp1),aireu(iip1,jjp1)           , &
-     & airev(iip1,jjm),unsaire(iip1,jjp1),apoln,apols                 , &
-     & unsairez(iip1,jjm),airuscv2(iip1,jjm),airvscu2(iip1,jjm)       , &
-     & aireij1(iip1,jjp1),aireij2(iip1,jjp1),aireij3(iip1,jjp1)       , &
-     & aireij4(iip1,jjp1),alpha1(iip1,jjp1),alpha2(iip1,jjp1)         , &
-     & alpha3(iip1,jjp1),alpha4(iip1,jjp1),alpha1p2(iip1,jjp1)        , &
-     & alpha1p4(iip1,jjp1),alpha2p3(iip1,jjp1),alpha3p4(iip1,jjp1)    , &
-     & fext(iip1,jjm),constang(iip1,jjp1), rlatu(jjp1),rlatv(jjm),      &
-     & rlonu(iip1),rlonv(iip1),cuvsurcv(iip1,jjm),cvsurcuv(iip1,jjm)  , &
-     & cvusurcu(iip1,jjp1),cusurcvu(iip1,jjp1)                        , &
-     & cuvscvgam1(iip1,jjm),cuvscvgam2(iip1,jjm),cvuscugam1(iip1,jjp1), &
-     & cvuscugam2(iip1,jjp1),cvscuvgam(iip1,jjm),cuscvugam(iip1,jjp1) , &
-     & unsapolnga1,unsapolnga2,unsapolsga1,unsapolsga2                , &
-     & unsair_gam1(iip1,jjp1),unsair_gam2(iip1,jjp1)                  , &
-     & unsairz_gam(iip1,jjm),aivscu2gam(iip1,jjm),aiuscv2gam(iip1,jjm)  &
-     & , xprimu(iip1),xprimv(iip1)
-
-
-      REAL                                                               &
-     & cu,cv,unscu2,unscv2,aire,airesurg,aireu,airev,apoln,apols,unsaire &
-     & ,unsairez,airuscv2,airvscu2,aireij1,aireij2,aireij3,aireij4     , &
-     & alpha1,alpha2,alpha3,alpha4,alpha1p2,alpha1p4,alpha2p3,alpha3p4 , &
-     & fext,constang,rlatu,rlatv,rlonu,rlonv,cuvscvgam1,cuvscvgam2     , &
-     & cvuscugam1,cvuscugam2,cvscuvgam,cuscvugam,unsapolnga1           , &
-     & unsapolnga2,unsapolsga1,unsapolsga2,unsair_gam1,unsair_gam2     , &
-     & unsairz_gam,aivscu2gam,aiuscv2gam,cuvsurcv,cvsurcuv,cvusurcu    , &
-     & cusurcvu,xprimu,xprimv
Index: /LMDZ6/trunk/libf/dyn3d_common/comgeom2_mod_h.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/comgeom2_mod_h.f90	(revision 5281)
+++ /LMDZ6/trunk/libf/dyn3d_common/comgeom2_mod_h.f90	(revision 5281)
@@ -0,0 +1,49 @@
+! Replaces comgeom2.h
+
+! AB 07/24: NOTE: THIS USES A LEGACY COMMON BLOCK SHARED WITH comgeom_mod_h.f90 !
+! See comgeom_mod_h.f90 comment.
+
+MODULE comgeom2_mod_h
+  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
+  USE paramet_mod_h
+  IMPLICIT NONE; PRIVATE
+  PUBLIC                                                               &
+          cu, cv, unscu2, unscv2, aire, airesurg, aireu, airev, apoln, apols, unsaire &
+          , unsairez, airuscv2, airvscu2, aireij1, aireij2, aireij3, aireij4, &
+          alpha1, alpha2, alpha3, alpha4, alpha1p2, alpha1p4, alpha2p3, alpha3p4, &
+          fext, constang, rlatu, rlatv, rlonu, rlonv, cuvscvgam1, cuvscvgam2, &
+          cvuscugam1, cvuscugam2, cvscuvgam, cuscvugam, unsapolnga1, &
+          unsapolnga2, unsapolsga1, unsapolsga2, unsair_gam1, unsair_gam2, &
+          unsairz_gam, aivscu2gam, aiuscv2gam, cuvsurcv, cvsurcuv, cvusurcu, &
+          cusurcvu, xprimu, xprimv
+
+  COMMON/comgeom/                                                   &
+          cu(iip1, jjp1), cv(iip1, jjm), unscu2(iip1, jjp1), unscv2(iip1, jjm), &
+          aire(iip1, jjp1), airesurg(iip1, jjp1), aireu(iip1, jjp1), &
+          airev(iip1, jjm), unsaire(iip1, jjp1), apoln, apols, &
+          unsairez(iip1, jjm), airuscv2(iip1, jjm), airvscu2(iip1, jjm), &
+          aireij1(iip1, jjp1), aireij2(iip1, jjp1), aireij3(iip1, jjp1), &
+          aireij4(iip1, jjp1), alpha1(iip1, jjp1), alpha2(iip1, jjp1), &
+          alpha3(iip1, jjp1), alpha4(iip1, jjp1), alpha1p2(iip1, jjp1), &
+          alpha1p4(iip1, jjp1), alpha2p3(iip1, jjp1), alpha3p4(iip1, jjp1), &
+          fext(iip1, jjm), constang(iip1, jjp1), rlatu(jjp1), rlatv(jjm), &
+          rlonu(iip1), rlonv(iip1), cuvsurcv(iip1, jjm), cvsurcuv(iip1, jjm), &
+          cvusurcu(iip1, jjp1), cusurcvu(iip1, jjp1), &
+          cuvscvgam1(iip1, jjm), cuvscvgam2(iip1, jjm), cvuscugam1(iip1, jjp1), &
+          cvuscugam2(iip1, jjp1), cvscuvgam(iip1, jjm), cuscvugam(iip1, jjp1), &
+          unsapolnga1, unsapolnga2, unsapolsga1, unsapolsga2, &
+          unsair_gam1(iip1, jjp1), unsair_gam2(iip1, jjp1), &
+          unsairz_gam(iip1, jjm), aivscu2gam(iip1, jjm), aiuscv2gam(iip1, jjm)  &
+          , xprimu(iip1), xprimv(iip1)
+
+  REAL                                                               &
+          cu, cv, unscu2, unscv2, aire, airesurg, aireu, airev, apoln, apols, unsaire &
+          , unsairez, airuscv2, airvscu2, aireij1, aireij2, aireij3, aireij4, &
+          alpha1, alpha2, alpha3, alpha4, alpha1p2, alpha1p4, alpha2p3, alpha3p4, &
+          fext, constang, rlatu, rlatv, rlonu, rlonv, cuvscvgam1, cuvscvgam2, &
+          cvuscugam1, cvuscugam2, cvscuvgam, cuscvugam, unsapolnga1, &
+          unsapolnga2, unsapolsga1, unsapolsga2, unsair_gam1, unsair_gam2, &
+          unsairz_gam, aivscu2gam, aiuscv2gam, cuvsurcv, cvsurcuv, cvusurcu, &
+          cusurcvu, xprimu, xprimv
+
+END MODULE comgeom2_mod_h
Index: /LMDZ6/trunk/libf/dyn3d_common/comgeom_mod_h.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/comgeom_mod_h.f90	(revision 5281)
+++ /LMDZ6/trunk/libf/dyn3d_common/comgeom_mod_h.f90	(revision 5281)
@@ -0,0 +1,51 @@
+! Replaces comgeom.h
+
+! AB 07/24: NOTE: THIS USES A LEGACY COMMON BLOCK SHARED WITH comgeom2_mod_h.f90 !
+! This file defines a number of 1D arrays.
+! The COMMON block is used to transparently (but using unsafe procedures) cast those as 2D arrays, without duplication.
+! Therefore there's no trivial way to turn this into a module without COMMON.
+! In the future we'd still like to get rid of this COMMON block.
+
+MODULE comgeom_mod_h
+  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
+  USE paramet_mod_h
+  IMPLICIT NONE; PRIVATE
+  PUBLIC cu, cv, unscu2, unscv2, aire, airesurg, aireu, airev, unsaire, apoln, &
+          apols, unsairez, airuscv2, airvscu2, aireij1, aireij2, aireij3, aireij4, &
+          alpha1, alpha2, alpha3, alpha4, alpha1p2, alpha1p4, alpha2p3, alpha3p4, &
+          fext, constang, rlatu, rlatv, rlonu, rlonv, cuvscvgam1, cuvscvgam2, &
+          cvuscugam1, cvuscugam2, cvscuvgam, cuscvugam, unsapolnga1, unsapolnga2&
+          , unsapolsga1, unsapolsga2, unsair_gam1, unsair_gam2, unsairz_gam, &
+          aivscu2gam, aiuscv2gam, cuvsurcv, cvsurcuv, cvusurcu, cusurcvu, xprimu&
+          , xprimv
+
+  COMMON/comgeom/                                                   &
+          cu(ip1jmp1), cv(ip1jm), unscu2(ip1jmp1), unscv2(ip1jm), &
+          aire(ip1jmp1), airesurg(ip1jmp1), aireu(ip1jmp1), &
+          airev(ip1jm), unsaire(ip1jmp1), apoln, apols, &
+          unsairez(ip1jm), airuscv2(ip1jm), airvscu2(ip1jm), &
+          aireij1(ip1jmp1), aireij2(ip1jmp1), aireij3(ip1jmp1), &
+          aireij4(ip1jmp1), alpha1(ip1jmp1), alpha2(ip1jmp1), &
+          alpha3(ip1jmp1), alpha4(ip1jmp1), alpha1p2(ip1jmp1), &
+          alpha1p4(ip1jmp1), alpha2p3(ip1jmp1), alpha3p4(ip1jmp1), &
+          fext(ip1jm), constang(ip1jmp1), rlatu(jjp1), rlatv(jjm), &
+          rlonu(iip1), rlonv(iip1), cuvsurcv(ip1jm), cvsurcuv(ip1jm), &
+          cvusurcu(ip1jmp1), cusurcvu(ip1jmp1), cuvscvgam1(ip1jm), &
+          cuvscvgam2(ip1jm), cvuscugam1(ip1jmp1), &
+          cvuscugam2(ip1jmp1), cvscuvgam(ip1jm), cuscvugam(ip1jmp1), &
+          unsapolnga1, unsapolnga2, unsapolsga1, unsapolsga2, &
+          unsair_gam1(ip1jmp1), unsair_gam2(ip1jmp1), unsairz_gam(ip1jm), &
+          aivscu2gam(ip1jm), aiuscv2gam(ip1jm), xprimu(iip1), xprimv(iip1)
+
+  REAL                                                            &
+          cu, cv, unscu2, unscv2, aire, airesurg, aireu, airev, unsaire, apoln, &
+          apols, unsairez, airuscv2, airvscu2, aireij1, aireij2, aireij3, aireij4, &
+          alpha1, alpha2, alpha3, alpha4, alpha1p2, alpha1p4, alpha2p3, alpha3p4, &
+          fext, constang, rlatu, rlatv, rlonu, rlonv, cuvscvgam1, cuvscvgam2, &
+          cvuscugam1, cvuscugam2, cvscuvgam, cuscvugam, unsapolnga1, unsapolnga2&
+          , unsapolsga1, unsapolsga2, unsair_gam1, unsair_gam2, unsairz_gam, &
+          aivscu2gam, aiuscv2gam, cuvsurcv, cvsurcuv, cvusurcu, cusurcvu, xprimu&
+          , xprimv
+
+END MODULE comgeom_mod_h
+
Index: /LMDZ6/trunk/libf/dyn3d_common/convflu.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/convflu.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/convflu.f90	(revision 5281)
@@ -18,4 +18,5 @@
   ! nbniv   est le nombre de niveaux vert. de  xflu et de yflu .
   !
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -33,5 +34,4 @@
   !
   !
-  INCLUDE "comgeom.h"
   !
   DO l = 1,nbniv
Index: /LMDZ6/trunk/libf/dyn3d_common/convmas.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/convmas.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/convmas.f90	(revision 5281)
@@ -5,4 +5,5 @@
 !-------------------------------------------------------------------------------
 ! Purpose: Compute mass flux convergence at p levels.
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -11,5 +12,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: /LMDZ6/trunk/libf/dyn3d_common/coordij.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/coordij.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/coordij.f90	(revision 5281)
@@ -11,4 +11,5 @@
   !=======================================================================
 
+  USE comgeom_mod_h
   USE comconst_mod, ONLY: pi
 
@@ -23,5 +24,4 @@
 
 
-  include "comgeom.h"
 
   real :: zlon,zlat
Index: /LMDZ6/trunk/libf/dyn3d_common/covcont.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/covcont.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/covcont.f90	(revision 5281)
@@ -6,4 +6,5 @@
 ! Purpose: Compute contravariant components from covariant components.
 !-------------------------------------------------------------------------------
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -12,5 +13,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: /LMDZ6/trunk/libf/dyn3d_common/diagedyn.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/diagedyn.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/diagedyn.f90	(revision 5281)
@@ -53,4 +53,5 @@
   !======================================================================
 
+  USE comgeom_mod_h
   USE control_mod, ONLY : planet_type
 
@@ -62,5 +63,4 @@
 
 
-  INCLUDE "comgeom.h"
   INCLUDE "iniprint.h"
 
Index: /LMDZ6/trunk/libf/dyn3d_common/diverg.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/diverg.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/diverg.f90	(revision 5281)
@@ -11,4 +11,5 @@
   !          x et y  etant des composantes covariantes   ...
   !  *********************************************************************
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -27,5 +28,4 @@
 
 
-  INCLUDE "comgeom.h"
   !
   !    ..........          variables en arguments    ...................
Index: /LMDZ6/trunk/libf/dyn3d_common/diverg_gam.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/diverg_gam.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/diverg_gam.f90	(revision 5281)
@@ -12,4 +12,5 @@
   !          x et y  etant des composantes covariantes   ...
   !  *********************************************************************
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -28,5 +29,4 @@
 
 
-  INCLUDE "comgeom.h"
   !
   !    ..........          variables en arguments    ...................
Index: /LMDZ6/trunk/libf/dyn3d_common/divergf.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/divergf.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/divergf.f90	(revision 5281)
@@ -11,4 +11,5 @@
   !          x et y  etant des composantes covariantes   ...
   !  *********************************************************************
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -27,5 +28,4 @@
 
 
-  INCLUDE "comgeom.h"
   !
   !    ..........          variables en arguments    ...................
Index: /LMDZ6/trunk/libf/dyn3d_common/divergst.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/divergst.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/divergst.f90	(revision 5281)
@@ -3,4 +3,5 @@
 !
 SUBROUTINE divergst(klevel,x,y,div)
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
   USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -21,5 +22,4 @@
   !   -------------------------------------------------------------------
   !
-  INCLUDE "comgeom.h"
 
   INTEGER :: klevel
Index: /LMDZ6/trunk/libf/dyn3d_common/divgrad.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/divgrad.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/divgrad.f90	(revision 5281)
@@ -3,4 +3,5 @@
 !
 SUBROUTINE divgrad (klevel,h, lh, divgra )
+  USE comgeom_mod_h
   USE comdissipn_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -26,5 +27,4 @@
 
 
-  include "comgeom.h"
   !
   INTEGER :: klevel
Index: /LMDZ6/trunk/libf/dyn3d_common/divgrad2.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/divgrad2.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/divgrad2.f90	(revision 5281)
@@ -13,4 +13,5 @@
   !     divgra     est  un argument  de sortie pour le s-prg
   !
+  USE comgeom2_mod_h
   USE comdissipn_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -21,5 +22,4 @@
 
 
-  INCLUDE "comgeom2.h"
 
   !    .......    variables en arguments   .......
Index: /LMDZ6/trunk/libf/dyn3d_common/enercin.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/enercin.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/enercin.f90	(revision 5281)
@@ -5,4 +5,5 @@
 !-------------------------------------------------------------------------------
 ! Purpose: Compute kinetic energy at sigma levels.
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -11,5 +12,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: /LMDZ6/trunk/libf/dyn3d_common/exner_hyb_m.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/exner_hyb_m.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/exner_hyb_m.f90	(revision 5281)
@@ -1,5 +1,6 @@
 module exner_hyb_m
 
-  IMPLICIT NONE
+  USE comgeom_mod_h
+    IMPLICIT NONE
 
 contains
@@ -14,5 +15,5 @@
     !
     !   ************************************************************************
-    !    Calcule la fonction d'Exner pk = Cp * (p/preff) ** kappa , aux milieux des 
+    !    Calcule la fonction d'Exner pk = Cp * (p/preff) ** kappa , aux milieux des
     !    couches .   Pk(l) sera calcule aux milieux  des couches l ,entre les
     !    pressions p(l) et p(l+1) ,definis aux interfaces des llm couches .
@@ -26,7 +27,7 @@
     !    ( voir note de Fr.Hourdin )  ,
     !
-    !    on determine successivement , du haut vers le bas des couches, les 
-    !    coef. alpha(llm),beta(llm) .,.,alpha(l),beta(l),,,alpha(2),beta(2), 
-    !    puis pk(ij,1). Ensuite ,on calcule,du bas vers le haut des couches,  
+    !    on determine successivement , du haut vers le bas des couches, les
+    !    coef. alpha(llm),beta(llm) .,.,alpha(l),beta(l),,,alpha(2),beta(2),
+    !    puis pk(ij,1). Ensuite ,on calcule,du bas vers le haut des couches,
     !     pk(ij,l)  donne  par la relation (2),  pour l = 2 a l = llm .
     !
@@ -35,5 +36,5 @@
     USE comconst_mod, ONLY: jmp1, cpp, kappa, r
     USE comvert_mod, ONLY: preff
-    
+
     USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -43,5 +44,4 @@
 
 
-    include "comgeom.h"
 
     INTEGER  ngrid
Index: /LMDZ6/trunk/libf/dyn3d_common/exner_milieu_m.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/exner_milieu_m.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/exner_milieu_m.f90	(revision 5281)
@@ -1,5 +1,6 @@
 module exner_milieu_m
 
-  IMPLICIT NONE
+  USE comgeom_mod_h
+    IMPLICIT NONE
 
 contains
@@ -15,5 +16,5 @@
     !
     !   ************************************************************************
-    !    Calcule la fonction d'Exner pk = Cp * (p/preff) ** kappa , aux milieux des 
+    !    Calcule la fonction d'Exner pk = Cp * (p/preff) ** kappa , aux milieux des
     !    couches .   Pk(l) sera calcule aux milieux  des couches l ,entre les
     !    pressions p(l) et p(l+1) ,definis aux interfaces des llm couches .
@@ -23,7 +24,7 @@
     !
     !     WARNING : CECI est une version speciale de exner_hyb originale
-    !               Utilise dans la version martienne pour pouvoir 
+    !               Utilise dans la version martienne pour pouvoir
     !               tourner avec des coordonnees verticales complexe
-    !              => Il ne verifie PAS la condition la proportionalite en 
+    !              => Il ne verifie PAS la condition la proportionalite en
     !              energie totale/ interne / potentielle (F.Forget 2001)
     !    ( voir note de Fr.Hourdin )  ,
@@ -32,5 +33,5 @@
     USE comconst_mod, ONLY: jmp1, cpp, kappa, r
     USE comvert_mod, ONLY: preff
-    
+
     USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -40,5 +41,4 @@
 
 
-    include "comgeom.h"
 
     INTEGER  ngrid
Index: /LMDZ6/trunk/libf/dyn3d_common/flumass.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/flumass.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/flumass.f90	(revision 5281)
@@ -5,4 +5,5 @@
 !-------------------------------------------------------------------------------
 ! Purpose: Compute mass flux at s levels.
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -11,5 +12,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: /LMDZ6/trunk/libf/dyn3d_common/gr_u_scal.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/gr_u_scal.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/gr_u_scal.f90	(revision 5281)
@@ -25,4 +25,5 @@
   !
   !=======================================================================
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -35,5 +36,4 @@
 
 
-  INCLUDE "comgeom.h"
 
   !   Arguments:
Index: /LMDZ6/trunk/libf/dyn3d_common/gr_v_scal.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/gr_v_scal.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/gr_v_scal.f90	(revision 5281)
@@ -25,4 +25,5 @@
   !
   !=======================================================================
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -35,5 +36,4 @@
 
 
-  INCLUDE "comgeom.h"
 
   !   Arguments:
Index: /LMDZ6/trunk/libf/dyn3d_common/gradiv2.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/gradiv2.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/gradiv2.f90	(revision 5281)
@@ -16,4 +16,5 @@
   !
   !
+  USE comgeom_mod_h
   USE comdissipn_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -24,5 +25,4 @@
 
 
-  INCLUDE "comgeom.h"
   !
   ! ........    variables en arguments      ........
Index: /LMDZ6/trunk/libf/dyn3d_common/grilles_gcm_netcdf_sub.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/grilles_gcm_netcdf_sub.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/grilles_gcm_netcdf_sub.f90	(revision 5281)
@@ -12,9 +12,10 @@
 SUBROUTINE grilles_gcm_netcdf_sub(masque,phis)
 
+  USE comgeom_mod_h
   USE comconst_mod, ONLY: cpp, kappa, g, omeg, daysec, rad, pi
   USE comvert_mod, ONLY: presnivs, preff, pa
   USE netcdf, ONLY: nf90_def_var, nf90_int, nf90_float, nf90_put_var, nf90_clobber, nf90_64bit_offset, nf90_def_dim, &
           nf90_put_att, nf90_enddef, nf90_create
-  
+
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -24,5 +25,4 @@
 
 
-  INCLUDE "comgeom.h"
 
 !======================== 
Index: /LMDZ6/trunk/libf/dyn3d_common/inigeom.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/inigeom.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/inigeom.f90	(revision 5281)
@@ -16,4 +16,5 @@
   !
   !
+  USE comgeom2_mod_h
   USE comdissnew_mod_h
   use fxhyp_m, only: fxhyp
@@ -30,5 +31,4 @@
 
 
-  include "comgeom2.h"
 
   !-----------------------------------------------------------------------
Index: /LMDZ6/trunk/libf/dyn3d_common/initdynav.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/initdynav.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/initdynav.f90	(revision 5281)
@@ -3,4 +3,5 @@
 subroutine initdynav(day0,anne0,tstep,t_ops,t_wrt)
 
+  USE comgeom_mod_h
   USE IOIPSL
   USE infotrac, ONLY : nqtot
@@ -10,5 +11,5 @@
   USE comvert_mod, ONLY: presnivs
   USE temps_mod, ONLY: itau_dyn
-  
+
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -39,5 +40,4 @@
 
 
-  include "comgeom.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3d_common/initfluxsto.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/initfluxsto.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/initfluxsto.f90	(revision 5281)
@@ -6,5 +6,6 @@
         fileid,filevid,filedid)
 
-   USE IOIPSL
+   USE comgeom_mod_h
+  USE IOIPSL
   USE comconst_mod, ONLY: pi
   USE comvert_mod, ONLY: nivsigs
@@ -45,5 +46,4 @@
 
 
-  include "comgeom.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3d_common/inithist.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/inithist.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/inithist.f90	(revision 5281)
@@ -4,5 +4,6 @@
 subroutine inithist(day0,anne0,tstep,t_ops,t_wrt)
 
-   USE IOIPSL
+   USE comgeom_mod_h
+  USE IOIPSL
    USE infotrac, ONLY : nqtot
    use com_io_dyn_mod, only : histid,histvid,histuid,               &
@@ -44,5 +45,4 @@
 
 
-  include "comgeom.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3d_common/inter_barxy_m.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/inter_barxy_m.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/inter_barxy_m.f90	(revision 5281)
@@ -6,5 +6,6 @@
   ! Authors: Robert SADOURNY, Phu LE VAN, Lionel GUEZ
 
-  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
+  USE comgeom2_mod_h
+    USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
@@ -27,5 +28,4 @@
     ! (for other included files)
 
-    include "comgeom2.h"
     ! (for "aire", "apoln", "apols")
 
Index: /LMDZ6/trunk/libf/dyn3d_common/interpost.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/interpost.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/interpost.f90	(revision 5281)
@@ -4,5 +4,6 @@
   subroutine interpost(q,qppm)
 
-   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
+   USE comgeom2_mod_h
+  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
@@ -12,5 +13,4 @@
 
 
-  include "comgeom2.h"
 
   ! Arguments
Index: /LMDZ6/trunk/libf/dyn3d_common/interpre.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/interpre.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/interpre.f90	(revision 5281)
@@ -6,4 +6,5 @@
          unatppm,vnatppm,psppm)
 
+  USE comgeom2_mod_h
   USE comdissip_mod_h
   USE comconst_mod, ONLY: g
@@ -17,5 +18,4 @@
 
 
-  include "comgeom2.h"
   include "description.h"
 
Index: /LMDZ6/trunk/libf/dyn3d_common/laplacien.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/laplacien.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/laplacien.f90	(revision 5281)
@@ -12,4 +12,5 @@
   !  divgra     est  un argument  de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -19,5 +20,4 @@
 
 
-  INCLUDE "comgeom.h"
 
   !
Index: /LMDZ6/trunk/libf/dyn3d_common/laplacien_gam.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/laplacien_gam.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/laplacien_gam.f90	(revision 5281)
@@ -14,4 +14,5 @@
   !  divgra     est  un argument  de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -21,5 +22,4 @@
 
 
-  INCLUDE "comgeom.h"
 
   !
Index: /LMDZ6/trunk/libf/dyn3d_common/laplacien_rot.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/laplacien_rot.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/laplacien_rot.f90	(revision 5281)
@@ -13,4 +13,5 @@
   !  rotout           est  un argument  de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -20,5 +21,4 @@
 
 
-  INCLUDE "comgeom.h"
 
   !
Index: /LMDZ6/trunk/libf/dyn3d_common/laplacien_rotgam.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/laplacien_rotgam.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/laplacien_rotgam.f90	(revision 5281)
@@ -12,4 +12,5 @@
   !  divgra     est  un argument  de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -19,5 +20,4 @@
 
 
-  INCLUDE "comgeom.h"
 
   !
Index: /LMDZ6/trunk/libf/dyn3d_common/limx.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/limx.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/limx.f90	(revision 5281)
@@ -13,4 +13,5 @@
   !
   !   --------------------------------------------------------------------
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -20,5 +21,4 @@
 
 
-  include "comgeom.h"
   !
   !
Index: /LMDZ6/trunk/libf/dyn3d_common/limy.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/limy.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/limy.f90	(revision 5281)
@@ -14,4 +14,5 @@
   !
   !   --------------------------------------------------------------------
+  USE comgeom_mod_h
   USE comconst_mod, ONLY: pi
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -22,5 +23,4 @@
 
 
-  include "comgeom.h"
   !
   !
Index: /LMDZ6/trunk/libf/dyn3d_common/limz.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/limz.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/limz.f90	(revision 5281)
@@ -13,4 +13,5 @@
   !
   !   --------------------------------------------------------------------
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -20,5 +21,4 @@
 
 
-  include "comgeom.h"
   !
   !
Index: /LMDZ6/trunk/libf/dyn3d_common/massbar.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/massbar.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/massbar.f90	(revision 5281)
@@ -6,4 +6,5 @@
 ! Purpose: Compute air mass mean along X and Y in each cell.
 ! See iniconst for more details.
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -12,5 +13,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: /LMDZ6/trunk/libf/dyn3d_common/massbarxy.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/massbarxy.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/massbarxy.f90	(revision 5281)
@@ -6,4 +6,5 @@
 ! Purpose: Compute air mass mean along X and Y in each cell.
 ! See iniconst for more details.
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -12,5 +13,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: /LMDZ6/trunk/libf/dyn3d_common/massdair.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/massdair.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/massdair.f90	(revision 5281)
@@ -16,4 +16,5 @@
   !  ....  p est defini aux interfaces des llm couches   .....
   !
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -23,5 +24,4 @@
 
 
-  include "comgeom.h"
   !
   !  .....   arguments  ....
Index: /LMDZ6/trunk/libf/dyn3d_common/nxgrad.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/nxgrad.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/nxgrad.f90	(revision 5281)
@@ -12,4 +12,5 @@
   !   x  et y    sont des arguments de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -19,5 +20,4 @@
 
 
-  INCLUDE "comgeom.h"
   INTEGER :: klevel
   REAL :: rot( ip1jm,klevel ),x( ip1jmp1,klevel ),y(ip1jm,klevel )
Index: /LMDZ6/trunk/libf/dyn3d_common/nxgrad_gam.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/nxgrad_gam.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/nxgrad_gam.f90	(revision 5281)
@@ -12,4 +12,5 @@
   !   x  et y    sont des arguments de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -19,5 +20,4 @@
 
 
-  INCLUDE "comgeom.h"
   INTEGER :: klevel
   REAL :: rot( ip1jm,klevel ),x( ip1jmp1,klevel ),y(ip1jm,klevel )
Index: /LMDZ6/trunk/libf/dyn3d_common/nxgradst.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/nxgradst.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/nxgradst.f90	(revision 5281)
@@ -3,4 +3,5 @@
 !
 SUBROUTINE nxgradst (klevel,rot, x, y )
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
   USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -16,5 +17,4 @@
   !
 
-  INCLUDE "comgeom.h"
 
   INTEGER :: klevel
Index: /LMDZ6/trunk/libf/dyn3d_common/pbar.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/pbar.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/pbar.f90	(revision 5281)
@@ -3,4 +3,5 @@
 !
 SUBROUTINE pbar ( pext, pbarx, pbary, pbarxy )
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
   USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -77,5 +78,4 @@
   !=======================================================================
 
-  INCLUDE "comgeom.h"
 
   REAL :: pext( ip1jmp1 ),  pbarx ( ip1jmp1 )
Index: /LMDZ6/trunk/libf/dyn3d_common/pentes_ini.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/pentes_ini.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/pentes_ini.f90	(revision 5281)
@@ -3,4 +3,5 @@
 !
 SUBROUTINE pentes_ini (q,w,masse,pbaru,pbarv,mode)
+  USE comgeom2_mod_h
   USE comconst_mod, ONLY: pi, dtvr
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -24,5 +25,4 @@
   !
   !=======================================================================
-  include "comgeom2.h"
 
   !   Arguments:
Index: /LMDZ6/trunk/libf/dyn3d_common/prather.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/prather.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/prather.f90	(revision 5281)
@@ -3,4 +3,5 @@
 !
 SUBROUTINE prather (q,w,masse,pbaru,pbarv,nt,dt)
+  USE comgeom2_mod_h
   USE comconst_mod, ONLY: pi
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -21,5 +22,4 @@
   !
   !=======================================================================
-  include "comgeom2.h"
 
   !   Arguments:
Index: /LMDZ6/trunk/libf/dyn3d_common/psextbar.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/psextbar.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/psextbar.f90	(revision 5281)
@@ -3,4 +3,5 @@
 !
 SUBROUTINE psextbar ( ps, psexbarxy )
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
   USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -76,5 +77,4 @@
   !
   !=======================================================================
-  INCLUDE "comgeom.h"
 
   REAL :: ps( ip1jmp1 ), psexbarxy ( ip1jm ), pext( ip1jmp1 )
Index: /LMDZ6/trunk/libf/dyn3d_common/rotat.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/rotat.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/rotat.f90	(revision 5281)
@@ -13,4 +13,5 @@
   !        rot          est  un argument  de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -20,5 +21,4 @@
 
 
-  INCLUDE "comgeom.h"
   !
   !   .....  variables en arguments  ......
Index: /LMDZ6/trunk/libf/dyn3d_common/rotat_nfil.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/rotat_nfil.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/rotat_nfil.f90	(revision 5281)
@@ -13,4 +13,5 @@
   !        rot          est  un argument  de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -20,5 +21,4 @@
 
 
-  INCLUDE "comgeom.h"
   !
   !   .....  variables en arguments  ......
Index: /LMDZ6/trunk/libf/dyn3d_common/rotatf.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/rotatf.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/rotatf.f90	(revision 5281)
@@ -14,4 +14,5 @@
   !        rot          est  un argument  de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -21,5 +22,4 @@
 
 
-  INCLUDE "comgeom.h"
   !
   !   .....  variables en arguments  ......
Index: /LMDZ6/trunk/libf/dyn3d_common/sortvarc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/sortvarc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/sortvarc.f90	(revision 5281)
@@ -6,4 +6,5 @@
         vcov )
 
+  USE comgeom_mod_h
   USE control_mod, ONLY: resetvarc
   USE comconst_mod, ONLY: dtvr, daysec, g, rad, omeg
@@ -35,5 +36,4 @@
 
 
-  INCLUDE "comgeom.h"
   INCLUDE "iniprint.h"
 
Index: /LMDZ6/trunk/libf/dyn3d_common/tourpot.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/tourpot.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/tourpot.f90	(revision 5281)
@@ -5,4 +5,5 @@
 !-------------------------------------------------------------------------------
 ! Purpose: Compute potential vorticity.
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -11,5 +12,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: /LMDZ6/trunk/libf/dyn3d_common/traceurpole.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/traceurpole.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/traceurpole.f90	(revision 5281)
@@ -4,5 +4,6 @@
     subroutine traceurpole(q,masse)
 
-      USE comdissip_mod_h
+      USE comgeom2_mod_h
+  USE comdissip_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -12,5 +13,4 @@
 
 
-  include "comgeom2.h"
   include "description.h"
 
Index: /LMDZ6/trunk/libf/dyn3d_common/ugeostr.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/ugeostr.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/ugeostr.f90	(revision 5281)
@@ -11,6 +11,7 @@
   ! levels are pressure levels.
 
+  USE comgeom2_mod_h
   use comconst_mod, only: omeg, rad
-  
+
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -20,5 +21,4 @@
 
 
-  include "comgeom2.h"
 
   real ucov(iip1,jjp1,llm),phi(iip1,jjp1,llm)
Index: /LMDZ6/trunk/libf/dyn3d_common/writedynav.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/writedynav.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/writedynav.f90	(revision 5281)
@@ -3,4 +3,5 @@
 subroutine writedynav(time, vcov, ucov, teta, ppk, phi, q, masse, ps, phis)
 
+  USE comgeom_mod_h
   USE ioipsl
   USE infotrac, ONLY : nqtot
@@ -33,5 +34,4 @@
 
 
-  include "comgeom.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3d_common/writehist.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3d_common/writehist.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3d_common/writehist.f90	(revision 5281)
@@ -3,4 +3,5 @@
 !
 subroutine writehist(time,vcov,ucov,teta,phi,q,masse,ps,phis)
+  USE comgeom_mod_h
   USE ioipsl
   USE infotrac, ONLY : nqtot
@@ -35,5 +36,4 @@
   !   Declarations
 
-  include "comgeom.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3dmem/addfi_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/addfi_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/addfi_loc.f90	(revision 5281)
@@ -11,4 +11,5 @@
   USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
+  USE comgeom_mod_h
   IMPLICIT NONE
   !
@@ -46,8 +47,4 @@
   !-----------------------------------------------------------------------
   !
-  !    0.  Declarations :
-  !    ------------------
-  include "comgeom.h"
-  !
   !    Arguments :
   !    -----------
Index: /LMDZ6/trunk/libf/dyn3dmem/advect_new_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/advect_new_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/advect_new_loc.f90	(revision 5281)
@@ -4,4 +4,5 @@
 SUBROUTINE advect_new_loc(ucov,vcov,teta,w,massebx,masseby, &
         du,dv,dteta)
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE write_field_loc
@@ -32,5 +33,4 @@
   !   Declarations:
   !   -------------
-  include "comgeom.h"
 
   !   Arguments:
Index: /LMDZ6/trunk/libf/dyn3dmem/advtrac_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/advtrac_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/advtrac_loc.f90	(revision 5281)
@@ -8,4 +8,5 @@
    !            M.A Filiberti (04/2002)
    !
+   USE comgeom2_mod_h
    USE comdissip_mod_h
    USE infotrac,     ONLY: nqtot, tracers
@@ -29,5 +30,4 @@
 
 
-   include "comgeom2.h"
    include "description.h"
 !   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3dmem/bilan_dyn_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/bilan_dyn_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/bilan_dyn_loc.f90	(revision 5281)
@@ -10,4 +10,5 @@
   !             vQ..A=Cp T + L * ...
 
+  USE comgeom2_mod_h
   USE IOIPSL
   USE parallel_lmdz
@@ -26,5 +27,4 @@
 
 
-  include "comgeom2.h"
   include "iniprint.h"
 
Index: /LMDZ6/trunk/libf/dyn3dmem/caldyn_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/caldyn_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/caldyn_loc.f90	(revision 5281)
@@ -2,4 +2,5 @@
         (itau,ucov,vcov,teta,ps,masse,pk,pkf,phis , &
         phi,conser,du,dv,dteta,dp,w,pbaru,pbarv,time )
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE Write_Field_loc
@@ -31,5 +32,4 @@
 
 
-  include "comgeom.h"
 
   !   Arguments:
Index: /LMDZ6/trunk/libf/dyn3dmem/call_dissip_mod.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/call_dissip_mod.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/call_dissip_mod.f90	(revision 5281)
@@ -47,10 +47,10 @@
     CALL allocate_u(dtetadis,llm,d)
     CALL allocate_u(dtetaecdt,llm,d)
-    
-    
+
+
     CALL dissip_allocate
-    
+
   END SUBROUTINE call_dissip_allocate
-  
+
   SUBROUTINE call_dissip_switch_dissip(dist)
   USE allocate_field_mod
@@ -76,9 +76,9 @@
 
     CALL dissip_switch_dissip(dist)
-    
-  END SUBROUTINE call_dissip_switch_dissip  
-  
-
-  
+
+  END SUBROUTINE call_dissip_switch_dissip
+
+
+
   SUBROUTINE call_dissip(ucov_dyn,vcov_dyn,teta_dyn,p_dyn,pk_dyn,ps_dyn)
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -92,6 +92,6 @@
   USE write_field_loc
   USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_DEBUGIO
+  USE comgeom_mod_h
   IMPLICIT NONE
-    INCLUDE 'comgeom.h'
     REAL,INTENT(INOUT) :: ucov_dyn(ijb_u:ije_u,llm) ! covariant zonal wind
     REAL,INTENT(INOUT) :: vcov_dyn(ijb_v:ije_v,llm) ! covariant meridional wind
Index: /LMDZ6/trunk/libf/dyn3dmem/convflu_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/convflu_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/convflu_loc.f90	(revision 5281)
@@ -15,4 +15,5 @@
   ! nbniv   est le nombre de niveaux vert. de  xflu et de yflu .
   !
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -33,5 +34,4 @@
   !
   !
-  INCLUDE "comgeom.h"
   !
 
Index: /LMDZ6/trunk/libf/dyn3dmem/convmas1_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/convmas1_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/convmas1_loc.f90	(revision 5281)
@@ -6,4 +6,5 @@
 ! Purpose: Compute mass flux convergence at p levels.
 !          Equivalent to convmas_loc if convmas2_loc is called after.
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE mod_filtreg_p
@@ -14,5 +15,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: /LMDZ6/trunk/libf/dyn3dmem/convmas2_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/convmas2_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/convmas2_loc.f90	(revision 5281)
@@ -6,4 +6,5 @@
 ! Purpose: Compute mass flux convergence at p levels.
 !          Equivalent to convmas_loc if convmas1_loc is called before.
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -13,5 +14,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: /LMDZ6/trunk/libf/dyn3dmem/convmas_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/convmas_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/convmas_loc.f90	(revision 5281)
@@ -5,4 +5,5 @@
 !-------------------------------------------------------------------------------
 ! Purpose: Compute mass flux convergence at p levels.
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE mod_filtreg_p
@@ -13,5 +14,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: /LMDZ6/trunk/libf/dyn3dmem/covcont_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/covcont_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/covcont_loc.f90	(revision 5281)
@@ -1,3 +1,4 @@
 SUBROUTINE covcont_loc (klevel,ucov, vcov, ucont, vcont )
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -19,5 +20,4 @@
   !
   !=======================================================================
-  INCLUDE "comgeom.h"
 
   INTEGER :: klevel
Index: /LMDZ6/trunk/libf/dyn3dmem/covnat_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/covnat_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/covnat_loc.f90	(revision 5281)
@@ -3,4 +3,5 @@
 !
 SUBROUTINE covnat_loc(klevel,ucov, vcov, unat, vnat )
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -22,5 +23,4 @@
   !
   !=======================================================================
-  INCLUDE "comgeom.h"
 
   INTEGER :: klevel
Index: /LMDZ6/trunk/libf/dyn3dmem/dissip_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/dissip_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/dissip_loc.f90	(revision 5281)
@@ -4,4 +4,5 @@
 SUBROUTINE dissip_loc( vcov,ucov,teta,p, dv,du,dh )
   !
+  USE comgeom_mod_h
   USE comdissipn_mod_h
   USE comdissnew_mod_h
@@ -37,5 +38,4 @@
 
 
-  include "comgeom.h"
 
   !   Arguments:
Index: /LMDZ6/trunk/libf/dyn3dmem/diverg_gam_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/diverg_gam_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/diverg_gam_loc.f90	(revision 5281)
@@ -9,4 +9,5 @@
   !          x et y  etant des composantes covariantes   ...
   !  *********************************************************************
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -26,5 +27,4 @@
 
 
-  INCLUDE "comgeom.h"
   !
   !    ..........          variables en arguments    ...................
Index: /LMDZ6/trunk/libf/dyn3dmem/diverg_p.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/diverg_p.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/diverg_p.f90	(revision 5281)
@@ -8,4 +8,5 @@
   !          x et y  etant des composantes covariantes   ...
   !  *********************************************************************
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -25,5 +26,4 @@
 
 
-  INCLUDE "comgeom.h"
   !
   !    ..........          variables en arguments    ...................
Index: /LMDZ6/trunk/libf/dyn3dmem/divergf_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/divergf_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/divergf_loc.f90	(revision 5281)
@@ -8,4 +8,5 @@
   !          x et y  etant des composantes covariantes   ...
   !  *********************************************************************
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE mod_filtreg_p
@@ -26,5 +27,4 @@
 
 
-  INCLUDE "comgeom.h"
   !
   !    ..........          variables en arguments    ...................
Index: /LMDZ6/trunk/libf/dyn3dmem/divgrad2_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/divgrad2_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/divgrad2_loc.f90	(revision 5281)
@@ -10,4 +10,5 @@
   !     divgra     est  un argument  de sortie pour le s-prg
   !
+  USE comgeom2_mod_h
   USE comdissipn_mod_h
   USE parallel_lmdz
@@ -22,5 +23,4 @@
 
 
-  INCLUDE "comgeom2.h"
 
   !    .......    variables en arguments   .......
Index: /LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.f90	(revision 5281)
@@ -6,4 +6,5 @@
 ! Purpose: Initial state reading.
 !-------------------------------------------------------------------------------
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE infotrac,    ONLY: nqtot, tracers, niso, iqIsoPha, iH2O, isoName, &
@@ -30,5 +31,4 @@
 
 
-  include "comgeom.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3dmem/dynredem_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/dynredem_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/dynredem_loc.f90	(revision 5281)
@@ -4,4 +4,5 @@
 ! Write the NetCDF restart file (initialization).
 !-------------------------------------------------------------------------------
+  USE comgeom_mod_h
   USE IOIPSL
   USE parallel_lmdz
@@ -23,10 +24,9 @@
 
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
-USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
+  USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
-IMPLICIT NONE
-
-
-  include "comgeom.h"
+  USE comgeom_mod_h
+  IMPLICIT NONE
+
   include "description.h"
   include "iniprint.h"
@@ -34,5 +34,5 @@
 ! Arguments:
   CHARACTER(LEN=*), INTENT(IN) :: fichnom          !--- FILE NAME
-  INTEGER,          INTENT(IN) :: iday_end         !--- 
+  INTEGER,          INTENT(IN) :: iday_end         !---
   REAL,             INTENT(IN) :: phis(ijb_u:ije_u)!--- GROUND GEOPOTENTIAL
 !===============================================================================
@@ -79,5 +79,5 @@
   tab_cntrl(19) = preff
 
-!    .....    parameters for zoom    ......   
+!    .....    parameters for zoom    ......
   tab_cntrl(20) = clon
   tab_cntrl(21) = clat
@@ -178,5 +178,5 @@
                           err, modname, fil, msg
   USE temps_mod, ONLY: itau_dyn, itaufin
-  
+
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -186,5 +186,4 @@
 
   include "description.h"
-  include "comgeom.h"
   include "iniprint.h"
 !===============================================================================
Index: /LMDZ6/trunk/libf/dyn3dmem/enercin_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/enercin_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/enercin_loc.f90	(revision 5281)
@@ -5,4 +5,5 @@
 !-------------------------------------------------------------------------------
 ! Purpose: Compute kinetic energy at sigma levels.
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -12,5 +13,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: /LMDZ6/trunk/libf/dyn3dmem/exner_hyb_loc_m.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/exner_hyb_loc_m.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/exner_hyb_loc_m.f90	(revision 5281)
@@ -1,5 +1,6 @@
 module exner_hyb_loc_m
 
-  IMPLICIT NONE
+  USE comgeom_mod_h
+    IMPLICIT NONE
 
 contains
@@ -14,5 +15,5 @@
     !
     !   ************************************************************************
-    !    Calcule la fonction d'Exner pk = Cp * (p/preff) ** kappa , aux milieux des 
+    !    Calcule la fonction d'Exner pk = Cp * (p/preff) ** kappa , aux milieux des
     !    couches .   Pk(l) sera calcule aux milieux  des couches l ,entre les
     !    pressions p(l) et p(l+1) ,definis aux interfaces des llm couches .
@@ -26,7 +27,7 @@
     !    ( voir note de Fr.Hourdin )  ,
     !
-    !    on determine successivement , du haut vers le bas des couches, les 
-    !    coef. alpha(llm),beta(llm) .,.,alpha(l),beta(l),,,alpha(2),beta(2), 
-    !    puis pk(ij,1). Ensuite ,on calcule,du bas vers le haut des couches,  
+    !    on determine successivement , du haut vers le bas des couches, les
+    !    coef. alpha(llm),beta(llm) .,.,alpha(l),beta(l),,,alpha(2),beta(2),
+    !    puis pk(ij,1). Ensuite ,on calcule,du bas vers le haut des couches,
     !     pk(ij,l)  donne  par la relation (2),  pour l = 2 a l = llm .
     !
@@ -45,5 +46,4 @@
 
 
-    include "comgeom.h"
 
     INTEGER  ngrid
Index: /LMDZ6/trunk/libf/dyn3dmem/exner_milieu_loc_m.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/exner_milieu_loc_m.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/exner_milieu_loc_m.f90	(revision 5281)
@@ -1,5 +1,6 @@
 module exner_milieu_loc_m
 
-  IMPLICIT NONE
+  USE comgeom_mod_h
+    IMPLICIT NONE
 
 contains
@@ -15,5 +16,5 @@
     !
     !   ************************************************************************
-    !    Calcule la fonction d'Exner pk = Cp * (p/preff) ** kappa , aux milieux des 
+    !    Calcule la fonction d'Exner pk = Cp * (p/preff) ** kappa , aux milieux des
     !    couches .   Pk(l) sera calcule aux milieux  des couches l ,entre les
     !    pressions p(l) et p(l+1) ,definis aux interfaces des llm couches .
@@ -23,7 +24,7 @@
     !
     !     WARNING : CECI est une version speciale de exner_hyb originale
-    !               Utilise dans la version martienne pour pouvoir 
+    !               Utilise dans la version martienne pour pouvoir
     !               tourner avec des coordonnees verticales complexe
-    !              => Il ne verifie PAS la condition la proportionalite en 
+    !              => Il ne verifie PAS la condition la proportionalite en
     !              energie totale/ interne / potentielle (F.Forget 2001)
     !    ( voir note de Fr.Hourdin )  ,
@@ -33,5 +34,5 @@
     USE comconst_mod, ONLY: cpp, kappa, r, jmp1
     USE comvert_mod, ONLY: preff
-    
+
     USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -41,5 +42,4 @@
 
 
-    include "comgeom.h"
 
     INTEGER  ngrid
Index: /LMDZ6/trunk/libf/dyn3dmem/flumass_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/flumass_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/flumass_loc.f90	(revision 5281)
@@ -5,4 +5,5 @@
 !-------------------------------------------------------------------------------
 ! Purpose: Compute mass flux at s levels.
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -12,5 +13,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: /LMDZ6/trunk/libf/dyn3dmem/friction_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/friction_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/friction_loc.f90	(revision 5281)
@@ -4,4 +4,5 @@
 !=======================================================================
 SUBROUTINE friction_loc(ucov,vcov,pdt)
+  USE comgeom2_mod_h
   USE parallel_lmdz
   USE control_mod
@@ -26,5 +27,4 @@
 
 
-  include "comgeom2.h"
   include "iniprint.h"
   include "academic.h"
Index: /LMDZ6/trunk/libf/dyn3dmem/gcm.F90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/gcm.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/gcm.F90	(revision 5281)
@@ -13,4 +13,5 @@
 !  USE mod_interface_dyn_phys, ONLY: init_interface_dyn_phys
 !#endif
+  USE comgeom_mod_h
   USE comdissnew_mod_h
   USE mod_hallo
@@ -61,5 +62,4 @@
   !   Declarations:
   !   -------------
-  include "comgeom.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3dmem/gr_u_scal_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/gr_u_scal_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/gr_u_scal_loc.f90	(revision 5281)
@@ -25,4 +25,5 @@
   !
   !=======================================================================
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -36,5 +37,4 @@
 
 
-  INCLUDE "comgeom.h"
 
   !   Arguments:
Index: /LMDZ6/trunk/libf/dyn3dmem/gr_v_scal_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/gr_v_scal_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/gr_v_scal_loc.f90	(revision 5281)
@@ -25,4 +25,5 @@
   !
   !=======================================================================
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -36,5 +37,4 @@
 
 
-  INCLUDE "comgeom.h"
 
   !   Arguments:
Index: /LMDZ6/trunk/libf/dyn3dmem/gradiv2_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/gradiv2_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/gradiv2_loc.f90	(revision 5281)
@@ -13,4 +13,5 @@
   !
   !
+  USE comgeom_mod_h
   USE comdissipn_mod_h
   USE parallel_lmdz
@@ -27,5 +28,4 @@
 
 
-  INCLUDE "comgeom.h"
   !
   ! ........    variables en arguments      ........
Index: /LMDZ6/trunk/libf/dyn3dmem/groupe_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/groupe_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/groupe_loc.f90	(revision 5281)
@@ -1,3 +1,4 @@
 subroutine groupe_loc(pext,pbaru,pbarv,pbarum,pbarvm,wm)
+  USE comgeom2_mod_h
   USE parallel_lmdz
   USE Write_field_loc
@@ -23,5 +24,4 @@
 
 
-  include "comgeom2.h"
 
   ! integer ngroup
Index: /LMDZ6/trunk/libf/dyn3dmem/groupeun_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/groupeun_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/groupeun_loc.f90	(revision 5281)
@@ -1,3 +1,4 @@
 SUBROUTINE groupeun_loc(jjmax,llmax,sb,se,jjb,jje,q)
+  USE comgeom2_mod_h
   USE parallel_lmdz
   USE Write_Field_p
@@ -10,5 +11,4 @@
 
 
-  include "comgeom2.h"
 
   INTEGER :: jjmax,llmax,sb,se,jjb,jje
@@ -138,4 +138,5 @@
 SUBROUTINE init_groupeun_loc(airen_tab, aires_tab)
 
+  USE comgeom2_mod_h
   USE parallel_lmdz
   USE comconst_mod, ONLY: ngroup
@@ -147,5 +148,4 @@
 
 
-  include "comgeom2.h"
 
   ! INTEGER ngroup
Index: /LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/guide_loc_mod.f90	(revision 5281)
@@ -787,8 +787,9 @@
   SUBROUTINE guide_zonave_u(typ,vsize,field)
 
+    USE comgeom_mod_h
     USE comconst_mod, ONLY: pi
 
     USE dimensions_mod, ONLY: iim, jjm, llm, ndm
-USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
+    USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
 IMPLICIT NONE
@@ -796,5 +797,4 @@
 
 
-    INCLUDE "comgeom.h"
 
     ! input/output variables
@@ -862,4 +862,5 @@
   SUBROUTINE guide_zonave_v(typ,hsize,vsize,field)
 
+    USE comgeom_mod_h
     USE comconst_mod, ONLY: pi
 
@@ -871,5 +872,4 @@
 
 
-    INCLUDE "comgeom.h"
 
     ! input/output variables
@@ -932,5 +932,6 @@
 !=======================================================================
   SUBROUTINE guide_interp(psi,teta)
-    use exner_hyb_loc_m, only: exner_hyb_loc
+    USE comgeom2_mod_h
+  use exner_hyb_loc_m, only: exner_hyb_loc
     use exner_milieu_loc_m, only: exner_milieu_loc
   USE parallel_lmdz
@@ -946,5 +947,4 @@
 
 
-  include "comgeom2.h"
 
   REAL, DIMENSION (iip1,jjb_u:jje_u),     INTENT(IN) :: psi ! Psol gcm
@@ -1431,11 +1431,9 @@
 
     USE dimensions_mod, ONLY: iim, jjm, llm, ndm
-USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
+    USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
-implicit none
-
-
-
-    include "comgeom2.h"
+    USE comgeom2_mod_h
+    implicit none
+
 
 ! input arguments :
@@ -2202,11 +2200,11 @@
 
     USE dimensions_mod, ONLY: iim, jjm, llm, ndm
-USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
+    USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
+    USE comgeom2_mod_h
 IMPLICIT NONE
 
 
 
-    INCLUDE "comgeom2.h"
 
     ! Variables entree
@@ -2424,5 +2422,5 @@
 
 subroutine dump2du(var,varname)
-use parallel_lmdz
+     use parallel_lmdz
 use mod_hallo
 USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -2460,10 +2458,10 @@
 subroutine dumpall
      USE dimensions_mod, ONLY: iim, jjm, llm, ndm
-USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
+     USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
+     USE comgeom_mod_h
 implicit none
 
 
-     include "comgeom.h"
      call barrier
      call dump2du(alpha_u(ijb_u:ije_u),'  alpha_u couche 1')
Index: /LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.f90	(revision 5281)
@@ -4,4 +4,5 @@
 SUBROUTINE iniacademic_loc(vcov,ucov,teta,q,masse,ps,phis,time_0)
 
+  USE comgeom_mod_h
   USE filtreg_mod, ONLY: inifilr
   USE infotrac,    ONLY: nqtot, niso, iqIsoPha, tracers, isoName, addPhase
@@ -37,5 +38,4 @@
 
 
-  include "comgeom.h"
   include "academic.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3dmem/initdynav_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/initdynav_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/initdynav_loc.f90	(revision 5281)
@@ -5,5 +5,6 @@
 
   ! This routine needs IOIPSL
-   USE IOIPSL
+   USE comgeom_mod_h
+  USE IOIPSL
 
    USE parallel_lmdz
@@ -49,5 +50,4 @@
 
 
-  include "comgeom.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3dmem/initfluxsto_p.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/initfluxsto_p.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/initfluxsto_p.f90	(revision 5281)
@@ -7,5 +7,6 @@
 
   ! This routine needs IOIPSL
-   USE IOIPSL
+   USE comgeom_mod_h
+  USE IOIPSL
 
    USE parallel_lmdz
@@ -50,5 +51,4 @@
 
 
-  include "comgeom.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3dmem/inithist_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/inithist_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/inithist_loc.f90	(revision 5281)
@@ -5,5 +5,6 @@
 
   ! This routine needs IOIPSL
-   USE IOIPSL
+   USE comgeom_mod_h
+  USE IOIPSL
 
    USE parallel_lmdz
@@ -47,5 +48,4 @@
 
 
-  include "comgeom.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3dmem/integrd_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/integrd_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/integrd_loc.f90	(revision 5281)
@@ -5,4 +5,5 @@
         (  nq,vcovm1,ucovm1,tetam1,psm1,massem1, &
         dv,du,dteta,dq,dp,vcov,ucov,teta,q,ps0,masse,phis) !,finvmaold)
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE control_mod
@@ -39,5 +40,4 @@
 
 
-  include "comgeom.h"
   include "iniprint.h"
 
Index: /LMDZ6/trunk/libf/dyn3dmem/laplacien_gam_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/laplacien_gam_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/laplacien_gam_loc.f90	(revision 5281)
@@ -11,4 +11,5 @@
   !  divgra     est  un argument  de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -19,5 +20,4 @@
 
 
-  INCLUDE "comgeom.h"
 
   !
Index: /LMDZ6/trunk/libf/dyn3dmem/laplacien_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/laplacien_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/laplacien_loc.f90	(revision 5281)
@@ -9,4 +9,5 @@
   !  divgra     est  un argument  de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE mod_filtreg_p
@@ -18,5 +19,4 @@
 
 
-  INCLUDE "comgeom.h"
 
   !
Index: /LMDZ6/trunk/libf/dyn3dmem/laplacien_rot_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/laplacien_rot_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/laplacien_rot_loc.f90	(revision 5281)
@@ -10,4 +10,5 @@
   !  rotout           est  un argument  de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE mod_filtreg_p
@@ -19,5 +20,4 @@
 
 
-  INCLUDE "comgeom.h"
 
   !
Index: /LMDZ6/trunk/libf/dyn3dmem/laplacien_rotgam_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/laplacien_rotgam_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/laplacien_rotgam_loc.f90	(revision 5281)
@@ -9,4 +9,5 @@
   !  divgra     est  un argument  de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -17,5 +18,4 @@
 
 
-  INCLUDE "comgeom.h"
 
   !
Index: /LMDZ6/trunk/libf/dyn3dmem/leapfrog_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/leapfrog_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/leapfrog_loc.f90	(revision 5281)
@@ -2,5 +2,6 @@
         masse0,phis0,q0,time_0)
 
-   USE comdissnew_mod_h
+   USE comgeom_mod_h
+  USE comdissnew_mod_h
    USE misc_mod
    USE parallel_lmdz
@@ -76,5 +77,4 @@
   !   Declarations:
   !   -------------
-  include "comgeom.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3dmem/massbar_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/massbar_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/massbar_loc.f90	(revision 5281)
@@ -6,4 +6,5 @@
 ! Purpose: Compute air mass mean along X and Y in each cell.
 ! See iniconst for more details.
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -13,5 +14,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: /LMDZ6/trunk/libf/dyn3dmem/massbarxy_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/massbarxy_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/massbarxy_loc.f90	(revision 5281)
@@ -6,4 +6,5 @@
 ! Purpose: Compute air mass mean along X and Y in each cell.
 ! See iniconst for more details.
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -13,5 +14,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: /LMDZ6/trunk/libf/dyn3dmem/massdair_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/massdair_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/massdair_loc.f90	(revision 5281)
@@ -1,3 +1,4 @@
 SUBROUTINE massdair_loc( p, masse )
+  USE comgeom_mod_h
   USE parallel_lmdz
   !
@@ -21,5 +22,4 @@
 
 
-  include "comgeom.h"
   !
   !  .....   arguments  ....
Index: /LMDZ6/trunk/libf/dyn3dmem/mod_xios_dyn3dmem.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/mod_xios_dyn3dmem.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/mod_xios_dyn3dmem.f90	(revision 5281)
@@ -37,4 +37,5 @@
 
 
+     USE comgeom_mod_h
      USE comvert_mod, ONLY: presnivs
      USE parallel_lmdz
@@ -47,5 +48,4 @@
 
 
-     INCLUDE 'comgeom.h'
 
      TYPE(xios_duration) :: tstep_xios
Index: /LMDZ6/trunk/libf/dyn3dmem/nxgrad_gam_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/nxgrad_gam_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/nxgrad_gam_loc.f90	(revision 5281)
@@ -9,4 +9,5 @@
   !   x  et y    sont des arguments de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE parallel_lmdz
 
@@ -18,5 +19,4 @@
 
 
-  INCLUDE "comgeom.h"
   INTEGER :: klevel
   REAL :: rot( ijb_v:ije_v,klevel )
Index: /LMDZ6/trunk/libf/dyn3dmem/nxgrad_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/nxgrad_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/nxgrad_loc.f90	(revision 5281)
@@ -9,4 +9,5 @@
   !   x  et y    sont des arguments de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -17,5 +18,4 @@
 
 
-  INCLUDE "comgeom.h"
   INTEGER :: klevel
   REAL :: rot( ijb_v:ije_v,klevel ),x( ijb_u:ije_u,klevel )
Index: /LMDZ6/trunk/libf/dyn3dmem/rotat_nfil_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/rotat_nfil_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/rotat_nfil_loc.f90	(revision 5281)
@@ -10,4 +10,5 @@
   !        rot          est  un argument  de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -18,5 +19,4 @@
 
 
-  INCLUDE "comgeom.h"
   !
   !   .....  variables en arguments  ......
Index: /LMDZ6/trunk/libf/dyn3dmem/rotat_p.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/rotat_p.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/rotat_p.f90	(revision 5281)
@@ -10,4 +10,5 @@
   !        rot          est  un argument  de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
@@ -18,5 +19,4 @@
 
 
-  INCLUDE "comgeom.h"
   !
   !   .....  variables en arguments  ......
Index: /LMDZ6/trunk/libf/dyn3dmem/rotatf_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/rotatf_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/rotatf_loc.f90	(revision 5281)
@@ -10,4 +10,5 @@
   !        rot          est  un argument  de sortie pour le s-prog
   !
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE mod_filtreg_p
@@ -19,5 +20,4 @@
 
 
-  INCLUDE "comgeom.h"
   !
   !   .....  variables en arguments  ......
Index: /LMDZ6/trunk/libf/dyn3dmem/sw_case_williamson91_6_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/sw_case_williamson91_6_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/sw_case_williamson91_6_loc.f90	(revision 5281)
@@ -26,4 +26,5 @@
   !
   !=======================================================================
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE comconst_mod, ONLY: cpp, omeg, rad
@@ -40,5 +41,4 @@
 
 
-  include "comgeom.h"
   include "iniprint.h"
 
Index: /LMDZ6/trunk/libf/dyn3dmem/top_bound_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/top_bound_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/top_bound_loc.f90	(revision 5281)
@@ -3,4 +3,5 @@
 !
 SUBROUTINE top_bound_loc(vcov,ucov,teta,masse,dt)
+  USE comgeom2_mod_h
   USE comdissipn_mod_h
   USE parallel_lmdz
@@ -16,5 +17,4 @@
 
 
-  include "comgeom2.h"
 
 
Index: /LMDZ6/trunk/libf/dyn3dmem/tourpot_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/tourpot_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/tourpot_loc.f90	(revision 5281)
@@ -5,4 +5,5 @@
 !-------------------------------------------------------------------------------
 ! Purpose: Compute potential vorticity.
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE mod_filtreg_p
@@ -13,5 +14,4 @@
 
 
-  include "comgeom.h"
 !===============================================================================
 ! Arguments:
Index: /LMDZ6/trunk/libf/dyn3dmem/vlsplt_loc.F90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/vlsplt_loc.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/vlsplt_loc.F90	(revision 5281)
@@ -417,4 +417,5 @@
   !
   !   --------------------------------------------------------------------
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE infotrac, ONLY : nqtot,tracers, & ! CRisi                 &
@@ -428,5 +429,4 @@
 
 
-  include "comgeom.h"
   !
   !
Index: /LMDZ6/trunk/libf/dyn3dmem/vlspltqs_loc.F90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/vlspltqs_loc.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/vlspltqs_loc.F90	(revision 5281)
@@ -422,4 +422,5 @@
   !
   !   --------------------------------------------------------------------
+  USE comgeom_mod_h
   USE parallel_lmdz
   USE infotrac, ONLY : nqtot,tracers, & ! CRisi                 &
@@ -433,5 +434,4 @@
 
 
-  include "comgeom.h"
   include "iniprint.h"
   !
Index: /LMDZ6/trunk/libf/dyn3dmem/writedyn_xios.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/writedyn_xios.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/writedyn_xios.f90	(revision 5281)
@@ -5,4 +5,5 @@
      &                           masse,ps,phis)
 
+      USE comgeom_mod_h
       USE lmdz_xios
       USE parallel_lmdz
@@ -13,5 +14,5 @@
       USE temps_mod, ONLY: itau_dyn
       USE mod_xios_dyn3dmem, ONLY : writefield_dyn_u, writefield_dyn_v
-      
+
       USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -39,5 +40,4 @@
 
 
-      include "comgeom.h"
       include "description.h"
       include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3dmem/writedynav_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/writedynav_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/writedynav_loc.f90	(revision 5281)
@@ -6,4 +6,5 @@
 
   ! This routine needs IOIPSL
+  USE comgeom_mod_h
   USE ioipsl
 
@@ -48,5 +49,4 @@
 
 
-  include "comgeom.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dyn3dmem/writehist_loc.f90
===================================================================
--- /LMDZ6/trunk/libf/dyn3dmem/writehist_loc.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dyn3dmem/writehist_loc.f90	(revision 5281)
@@ -6,4 +6,5 @@
 
   ! This routine needs IOIPSL
+  USE comgeom_mod_h
   USE ioipsl
 
@@ -48,5 +49,4 @@
 
 
-  include "comgeom.h"
   include "description.h"
   include "iniprint.h"
Index: /LMDZ6/trunk/libf/dynphy_lonlat/calfis.f90
===================================================================
--- /LMDZ6/trunk/libf/dynphy_lonlat/calfis.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dynphy_lonlat/calfis.f90	(revision 5281)
@@ -29,4 +29,5 @@
   !    Auteur :  P. Le Van, F. Hourdin
   !   .........
+  USE comgeom2_mod_h
   USE infotrac, ONLY: nqtot, tracers
   USE control_mod, ONLY: planet_type, nsplit_phys
@@ -94,5 +95,4 @@
   PARAMETER( ngridmx = 2+(jjm-1)*iim - 1/jjm   )
 
-  include "comgeom2.h"
   include "iniprint.h"
 
Index: /LMDZ6/trunk/libf/dynphy_lonlat/calfis_loc.F90
===================================================================
--- /LMDZ6/trunk/libf/dynphy_lonlat/calfis_loc.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/dynphy_lonlat/calfis_loc.F90	(revision 5281)
@@ -51,4 +51,5 @@
   USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_PHYS
 #ifdef CPP_PARA
+  USE comgeom2_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
   USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -111,5 +112,4 @@
   PARAMETER( ngridmx = 2+(jjm-1)*iim - 1/jjm   )
 
-  include "comgeom2.h"
   include "iniprint.h"
   !    Arguments :
Index: /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/ce0l.F90
===================================================================
--- /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/ce0l.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/ce0l.F90	(revision 5281)
@@ -39,4 +39,5 @@
 #endif
 
+  USE comgeom2_mod_h
   USE comconst_mod, ONLY: cpp, daysec, dtphys, dtvr, g, kappa, omeg, r, rad, &
                           pi, jmp1
@@ -55,5 +56,4 @@
 
 
-  include "comgeom2.h"
   include "iniprint.h"
   
Index: /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0dyn_netcdf.F90
===================================================================
--- /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0dyn_netcdf.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0dyn_netcdf.F90	(revision 5281)
@@ -45,4 +45,5 @@
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
   USE comdissnew_mod_h
+  USE comgeom2_mod_h
 IMPLICIT NONE
 
@@ -53,5 +54,4 @@
 
 
-  include "comgeom2.h"
   REAL, SAVE :: deg2rad
   INTEGER,            SAVE      :: iml_dyn, jml_dyn, llm_dyn, ttm_dyn, fid_dyn
Index: /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.f90
===================================================================
--- /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.f90	(revision 5281)
@@ -57,9 +57,9 @@
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
   USE dimsoil_mod_h, ONLY: nsoilmx
+  USE comgeom2_mod_h
   PRIVATE
   PUBLIC :: etat0phys_netcdf
 
   include "iniprint.h"
-  include "comgeom2.h"
   include "clesphys.h"
   REAL, SAVE :: deg2rad
Index: /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/iniphysiq_mod.F90
===================================================================
--- /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/iniphysiq_mod.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/iniphysiq_mod.F90	(revision 5281)
@@ -38,5 +38,6 @@
   USE bands, ONLY : distrib_phys
 #endif
-  USE mod_phys_lmdz_omp_data, ONLY: klon_omp 
+  USE comgeom_mod_h
+  USE mod_phys_lmdz_omp_data, ONLY: klon_omp
   USE ioipsl_getin_p_mod, ONLY: getin_p
   USE slab_heat_transp_mod, ONLY: ini_slab_transp_geom
@@ -56,5 +57,4 @@
   include "iniprint.h"
   include "tracstoke.h"
-  include "comgeom.h"
 
   REAL, INTENT (IN) :: prad ! radius of the planet (m)
Index: /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/init_ssrf_m.f90
===================================================================
--- /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/init_ssrf_m.f90	(revision 5280)
+++ /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/init_ssrf_m.f90	(revision 5281)
@@ -3,4 +3,5 @@
 !*******************************************************************************
 
+  USE comgeom2_mod_h
   USE indice_sol_mod,     ONLY: is_ter, is_oce, is_oce, is_lic, epsfra
   USE dimphy,             ONLY: klon, zmasq
@@ -20,5 +21,4 @@
   PUBLIC :: start_init_subsurf
   include "iniprint.h"
-  include "comgeom2.h"
 
 CONTAINS
Index: /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/limit_netcdf.F90
===================================================================
--- /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/limit_netcdf.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/dynphy_lonlat/phylmd/limit_netcdf.F90	(revision 5281)
@@ -70,4 +70,5 @@
 !-------------------------------------------------------------------------------
 #ifndef CPP_1D
+  USE comgeom2_mod_h
   USE indice_sol_mod
   USE netcdf,             ONLY: NF90_OPEN,    NF90_CREATE,  NF90_CLOSE,        &
@@ -82,7 +83,8 @@
   USE lmdz_cppkeys_wrapper, ONLY: nf90_format
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
-USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
+  USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
           ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
-IMPLICIT NONE
+  USE comgeom2_mod_h
+  IMPLICIT NONE
 !-------------------------------------------------------------------------------
 ! Arguments:
@@ -93,7 +95,4 @@
   REAL, DIMENSION(iip1,jjp1), INTENT(INOUT) :: phis   ! ground geopotential
   LOGICAL,                    INTENT(IN)    :: extrap ! SST extrapolation flag
-!-------------------------------------------------------------------------------
-! Local variables:
-  include "comgeom2.h"
 
 !--- INPUT NETCDF FILES AND VARIABLES NAMES ------------------------------------
@@ -339,5 +338,4 @@
 
 
-  include "comgeom2.h"
 !-----------------------------------------------------------------------------
 ! Arguments:
Index: /LMDZ6/trunk/libf/filtrez/filtreg_mod.F90
===================================================================
--- /LMDZ6/trunk/libf/filtrez/filtreg_mod.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/filtrez/filtreg_mod.F90	(revision 5281)
@@ -14,5 +14,6 @@
   USE mod_filtre_fft_loc, ONLY : Init_filtre_fft_loc=>Init_filtre_fft    !
 #endif
-  USE serre_mod, ONLY: alphax
+  USE comgeom_mod_h
+    USE serre_mod, ONLY: alphax
   USE logic_mod, ONLY: fxyhypb, ysinus
   USE comconst_mod, ONLY: maxlatfilter
@@ -32,5 +33,4 @@
 
     !  -------------------------------------------------------------------
-    include "comgeom.h"
     include "coefils.h"
 
Index: /LMDZ6/trunk/libf/filtrez/inifgn.F90
===================================================================
--- /LMDZ6/trunk/libf/filtrez/inifgn.F90	(revision 5280)
+++ /LMDZ6/trunk/libf/filtrez/inifgn.F90	(revision 5281)
@@ -6,4 +6,5 @@
   !    ...  H.Upadyaya , O.Sharma  ...
   !
+  USE comgeom_mod_h
   USE dimensions_mod, ONLY: iim, jjm, llm, ndm
 USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
@@ -13,5 +14,4 @@
 
 
-  include "comgeom.h"
 
   !
Index: DZ6/trunk/libf/phylmd/dyn1d/comgeom.h
===================================================================
--- /LMDZ6/trunk/libf/phylmd/dyn1d/comgeom.h	(revision 5280)
+++ 	(revision )
@@ -1,1 +1,0 @@
-link ../../dyn3d_common/comgeom.h
Index: DZ6/trunk/libf/phylmd/dyn1d/comgeom2.h
===================================================================
--- /LMDZ6/trunk/libf/phylmd/dyn1d/comgeom2.h	(revision 5280)
+++ 	(revision )
@@ -1,1 +1,0 @@
-link ../../dyn3d_common/comgeom2.h
Index: /LMDZ6/trunk/libf/phylmd/dyn1d/comgeom2_mod_h.f90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/dyn1d/comgeom2_mod_h.f90	(revision 5281)
+++ /LMDZ6/trunk/libf/phylmd/dyn1d/comgeom2_mod_h.f90	(revision 5281)
@@ -0,0 +1,1 @@
+link ../../dyn3d_common/comgeom2_mod_h.f90
Index: /LMDZ6/trunk/libf/phylmd/dyn1d/comgeom_mod_h.f90
===================================================================
--- /LMDZ6/trunk/libf/phylmd/dyn1d/comgeom_mod_h.f90	(revision 5281)
+++ /LMDZ6/trunk/libf/phylmd/dyn1d/comgeom_mod_h.f90	(revision 5281)
@@ -0,0 +1,1 @@
+link ../../dyn3d_common/comgeom_mod_h.f90
