--- trunk/Sources/phylmd/phytrac.f	2017/10/16 12:35:41	225
+++ trunk/phylmd/phytrac.f	2018/04/19 17:54:55	266
@@ -8,8 +8,8 @@
 contains
 
   SUBROUTINE phytrac(julien, gmtime, firstcal, lafin, pdtphys, t_seri, paprs, &
-       pplay, pmfu, pmfd, pde_u, pen_d, coefh, fm_therm, entr_therm, yu1, &
-       yv1, ftsol, pctsrf, frac_impa, frac_nucl, da, phi, mp, upwd, dnwd, &
+       pplay, pmfu, pmfd, pde_u, pen_d, coefh, cdragh, fm_therm, entr_therm, &
+       yu1, yv1, ftsol, pctsrf, frac_impa, frac_nucl, da, phi, mp, upwd, dnwd, &
        tr_seri, zmasse, ncid_startphy)
 
     ! From phylmd/phytrac.F, version 1.15, 2006/02/21 08:08:30 (SVN
@@ -35,7 +35,7 @@
     USE conf_gcm_m, ONLY: lmt_pas
     use ctherm, only: iflag_thermals
     use cvltr_m, only: cvltr
-    use dimens_m, only: llm, nqmx
+    use dimensions, only: llm, nqmx
     use dimphy, only: klon
     use histwrite_phy_m, only: histwrite_phy
     use indicesol, only: nbsrf
@@ -77,9 +77,10 @@
 
     REAL pde_u(klon, llm) ! flux detraine dans le panache montant
     REAL pen_d(klon, llm) ! flux entraine dans le panache descendant
-    REAL coefh(klon, llm) ! coeff melange couche limite
+    REAL coefh(:, 2:) ! (klon, 2:llm) coeff melange couche limite
+    real cdragh(:) ! (klon)
     real fm_therm(klon, llm+1), entr_therm(klon, llm) ! thermiques
-    REAL, intent(in):: yu1(klon), yv1(klon) ! vent au premier niveau
+    REAL, intent(in):: yu1(:), yv1(:) ! (klon) vent au premier niveau
 
     ! Arguments n\'ecessaires pour les sources et puits de traceur :
     real, intent(in):: ftsol(:, :) ! (klon, nbsrf) surface temperature (K)
@@ -118,21 +119,20 @@
 
     REAL, save:: trs(klon, nqmx - 2) ! Concentration de traceur dans le sol
 
-    REAL masktr(klon, nqmx - 2) ! Masque reservoir de sol traceur
+    REAL, save:: masktr(klon, nqmx - 2) ! Masque reservoir de sol traceur
     ! Masque de l'echange avec la surface
     ! (1 = reservoir) ou (possible => 1)
-    SAVE masktr
-    REAL fshtr(klon, nqmx - 2) ! Flux surfacique dans le reservoir de sol
-    SAVE fshtr
-    REAL hsoltr(nqmx - 2) ! Epaisseur equivalente du reservoir de sol
-    SAVE hsoltr
-    REAL tautr(nqmx - 2) ! Constante de decroissance radioactive
-    SAVE tautr
-    REAL vdeptr(nqmx - 2) ! Vitesse de depot sec dans la couche Brownienne
-    SAVE vdeptr
-    REAL scavtr(nqmx - 2) ! Coefficient de lessivage
-    SAVE scavtr
+    
+    REAL, save:: fshtr(klon, nqmx - 2)
+    ! Flux surfacique dans le reservoir de sol
+    
+    REAL, save:: hsoltr(nqmx - 2) ! Epaisseur equivalente du reservoir de sol
+    REAL, save:: tautr(nqmx - 2) ! constante de d\'ecroissance radioactive
 
+    REAL, save:: vdeptr(nqmx - 2)
+    ! Vitesse de depot sec dans la couche Brownienne
+
+    REAL, save:: scavtr(nqmx - 2) ! Coefficient de lessivage
     CHARACTER itn
 
     logical, save:: aerosol(nqmx - 2) ! Nature du traceur
@@ -279,7 +279,7 @@
     DO it=1, nqmx - 2
        if (clsol(it)) then 
           ! couche limite avec quantite dans le sol calculee
-          CALL cltracrn(it, pdtphys, yu1, yv1, coefh, t_seri, ftsol, &
+          CALL cltracrn(it, pdtphys, yu1, yv1, coefh, cdragh, t_seri, ftsol, &
                pctsrf, tr_seri(:, :, it), trs(:, it), paprs, pplay, delp, &
                masktr(1, it), fshtr(1, it), hsoltr(it), tautr(it), &
                vdeptr(it), rlat, d_tr_cl(1, 1, it), d_trs)