Index: LMDZ6/trunk/libf/dyn3d_common/infotrac.F90
===================================================================
--- LMDZ6/trunk/libf/dyn3d_common/infotrac.F90	(revision 4078)
+++ LMDZ6/trunk/libf/dyn3d_common/infotrac.F90	(revision 4082)
@@ -308,4 +308,5 @@
       CALL init_transport(hadv_inca, vadv_inca, conv_flg_inca, pbl_flg_inca, solsym_inca)
       ! DC passive CO2 tracer is at position 1: H2O was removed ; nqCO2/=0 in "inco" case only
+      ALLOCATE(conv_flg(nbtr),pbl_flg(nbtr),solsym(nbtr))
       conv_flg = [(  1,        ic=1, nqCO2),conv_flg_inca]
        pbl_flg = [(  1,        ic=1, nqCO2), pbl_flg_inca]
@@ -363,7 +364,4 @@
 
 #ifndef INCA
-   conv_flg = [(1, ic=1, nbtr)]                                      !--- Convection activated for all tracers
-    pbl_flg = [(1, ic=1, nbtr)]                                      !--- Boundary layer activated for all tracers
-   ALLOCATE(solsym(nbtr))
    CALL msg('Valeur de traceur.def :', modname)
    CALL msg('nombre total de traceurs '//TRIM(int2str(nqtrue)), modname)
@@ -374,4 +372,7 @@
    IF(planet_type == 'earth') nqo = COUNT(delPhase(tnom_0) == 'H2O') !--- for all planets except for Earth
    nbtr = nqtrue - nqo               
+   ALLOCATE(conv_flg(nbtr),pbl_flg(nbtr),solsym(nbtr))
+   conv_flg(1:nbtr) = 1                                     !--- Convection activated for all tracers
+   pbl_flg(1:nbtr) = 1                                     !--- Boundary layer activated for all tracers
 #endif
 
@@ -417,4 +418,5 @@
       CALL init_transport(hadv_inca, vadv_inca, conv_flg_inca, pbl_flg_inca, solsym_inca)
       ! DC passive CO2 tracer is at position 1: H2O was removed ; nqCO2/=0 in "inco" case only
+      
       conv_flg = [(  1       , k=1, nqCO2), conv_flg_inca]
       pbl_flg  = [(  1       , k=1, nqCO2), pbl_flg_inca]
@@ -440,4 +442,7 @@
       lerr = getKey('hadv', hadv, ky=tracers(:)%keys)
       lerr = getKey('vadv', vadv, ky=tracers(:)%keys)
+      ALLOCATE(solsym(nbtr))
+      conv_flg(1:nbtr)=1  !--- Convection activated for all tracers
+      pbl_flg(1:nbtr)=1   !--- Boundary layer activated for all tracers
 #endif
    !---------------------------------------------------------------------------------------------------------------------------
@@ -450,6 +455,6 @@
       lerr = getKey('vadv', vadv, ky=tracers(:)%keys)
       ALLOCATE(solsym(nbtr))
-      conv_flg = [(1, it=1, nbtr)]
-       pbl_flg = [(1, it=1, nbtr)]
+      conv_flg(1:nbtr)=1  !--- Convection activated for all tracers
+       pbl_flg(1:nbtr)=1  !--- Boundary layer activated for all tracers
    !---------------------------------------------------------------------------------------------------------------------------
    END IF