Index: LMDZ6/trunk/libf/dyn3d_common/infotrac.f90
===================================================================
--- LMDZ6/trunk/libf/dyn3d_common/infotrac.f90	(revision 5486)
+++ LMDZ6/trunk/libf/dyn3d_common/infotrac.f90	(revision 5487)
@@ -144,7 +144,4 @@
 ! Local variables
    INTEGER, ALLOCATABLE :: hadv(:), vadv(:)                          !--- Horizontal/vertical transport scheme number
-   INTEGER, ALLOCATABLE :: had (:), hadv_inca(:), conv_flg_inca(:), &!--- Variables specific to INCA
-                           vad (:), vadv_inca(:),  pbl_flg_inca(:)
-   CHARACTER(LEN=8), ALLOCATABLE :: solsym_inca(:)                   !--- Tracers names for INCA
    INTEGER :: nqINCA
    CHARACTER(LEN=2)      ::   suff(9)                                !--- Suffixes for schemes of order 3 or 4 (Prather)
Index: LMDZ6/trunk/libf/phylmd/cv3_routines.f90
===================================================================
--- LMDZ6/trunk/libf/phylmd/cv3_routines.f90	(revision 5486)
+++ LMDZ6/trunk/libf/phylmd/cv3_routines.f90	(revision 5487)
@@ -4963,5 +4963,5 @@
   LOGICAL,SAVE   ::  first=.TRUE.
   LOGICAL,SAVE   ::  keep_bug_indices_cv3_tracer
-  !$OMP_THREADPRIVATE (first, keep_bug_indices_cv3_tracer)
+!$OMP THREADPRIVATE(first, keep_bug_indices_cv3_tracer)
 
 ! variables d'Emanuel : du second indice au troisieme
Index: LMDZ6/trunk/libf/phylmd/infotrac_phy.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/infotrac_phy.F90	(revision 5486)
+++ LMDZ6/trunk/libf/phylmd/infotrac_phy.F90	(revision 5487)
@@ -158,7 +158,4 @@
 ! Local variables
    INTEGER, ALLOCATABLE :: hadv(:), vadv(:)                          !--- Horizontal/vertical transport scheme number
-   INTEGER, ALLOCATABLE :: had (:), hadv_inca(:), conv_flg_inca(:), &!--- Variables specific to INCA
-                           vad (:), vadv_inca(:),  pbl_flg_inca(:)
-   CHARACTER(LEN=8), ALLOCATABLE :: solsym_inca(:)                   !--- Tracers names for INCA
    INTEGER :: nqINCA
    CHARACTER(LEN=maxlen), ALLOCATABLE :: tnames(:)
@@ -257,8 +254,4 @@
    IF(CPPKEY_REPROBUS) CALL Init_chem_rep_trac(nbtr, nqo, tracers(:)%name)     !--- Transfert the number of tracers to Reprobus
 
-!##############################################################################################################################
-   IF(lInit) THEN                                                    !=== SKIPED IF ALREADY DONE IN dyn3d_common/infotrac  ####
-!##############################################################################################################################
-
 !==============================================================================================================================
 ! 2) Calculate nqtot, number of tracers needed (greater if advection schemes 20 or 30 have been chosen).
@@ -281,4 +274,8 @@
       CALL msg('The total number of tracers needed is '//TRIM(int2str(nqtot)))
    END IF
+
+!##############################################################################################################################
+   IF(lInit) THEN                                                    !=== SKIPED IF ALREADY DONE IN dyn3d_common/infotrac  ####
+!##############################################################################################################################
 
 !==============================================================================================================================
@@ -321,16 +318,4 @@
    IF(indexUpdate(tracers)) CALL abort_physic(modname, 'problem with tracers indices update', 1)
 
-!##############################################################################################################################
-   END IF
-!##############################################################################################################################
-
-!##############################################################################################################################
-   IF(.NOT.lInit) THEN
-!##############################################################################################################################
-     nqtot = SIZE(tracers)
-!##############################################################################################################################
-   ELSE
-!##############################################################################################################################
-
    !=== READ PHYSICAL PARAMETERS FOR ISOTOPES
    niso = 0; nzone = 0; nphas = nqo; ntiso = 0; isoCheck = .FALSE.
@@ -338,6 +323,16 @@
 
 !##############################################################################################################################
-   END IF
-!##############################################################################################################################
+   ELSE
+!##############################################################################################################################
+   DO iq = 1, nqtrue
+      t1 => tracers(iq)
+      IF(hadv(iq)     ==    vadv(iq)    ) iad = hadv(iq)
+      IF(hadv(iq)==10 .AND. vadv(iq)==16) iad = 11
+      tracers(iq)%isInPhysics= iad >= 0 .AND. (delPhase(t1%gen0Name) /= 'H2O' .OR. t1%component /= 'lmdz')
+   END DO
+!##############################################################################################################################
+   END IF
+!##############################################################################################################################
+
    !--- Convection / boundary layer activation for all tracers
    IF(.NOT.ALLOCATED(conv_flg)) ALLOCATE(conv_flg(nbtr)); conv_flg(1:nbtr) = 1