/[lmdze]/trunk/phylmd/phytrac.f
ViewVC logotype

Diff of /trunk/phylmd/phytrac.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 134 by guez, Wed Apr 29 15:47:56 2015 UTC revision 156 by guez, Thu Jul 16 17:39:10 2015 UTC
# Line 116  contains Line 116  contains
116      !      !
117      ! Pour la source de radon et son reservoir de sol      ! Pour la source de radon et son reservoir de sol
118    
119      REAL, save:: trs(klon, nqmx - 2) ! Concentration de radon dans le sol      REAL, save:: trs(klon, nqmx - 2) ! Concentration de traceur dans le sol
120    
121      REAL masktr(klon, nqmx - 2) ! Masque reservoir de sol traceur      REAL masktr(klon, nqmx - 2) ! Masque reservoir de sol traceur
122      ! Masque de l'echange avec la surface      ! Masque de l'echange avec la surface
# Line 155  contains Line 155  contains
155    
156      REAL d_tr(klon, llm), d_trs(klon) ! tendances de traceurs      REAL d_tr(klon, llm), d_trs(klon) ! tendances de traceurs
157      REAL d_tr_cl(klon, llm, nqmx - 2) ! tendance de traceurs couche limite      REAL d_tr_cl(klon, llm, nqmx - 2) ! tendance de traceurs couche limite
158      REAL d_tr_cv(klon, llm, nqmx - 2) ! tendance de traceurs conv pour chq traceur  
159        REAL d_tr_cv(klon, llm, nqmx - 2)
160        ! tendance de traceurs conv pour chq traceur
161    
162      REAL d_tr_th(klon, llm, nqmx - 2) ! la tendance des thermiques      REAL d_tr_th(klon, llm, nqmx - 2) ! la tendance des thermiques
163      REAL d_tr_dec(klon, llm, 2) ! la tendance de la decroissance      REAL d_tr_dec(klon, llm, 2) ! la tendance de la decroissance
164      ! ! radioactive du rn - > pb      ! ! radioactive du rn - > pb
# Line 296  contains Line 299  contains
299         DO it=1, nqmx - 2         DO it=1, nqmx - 2
300            if (clsol(it)) then            if (clsol(it)) then
301               ! couche limite avec quantite dans le sol calculee               ! couche limite avec quantite dans le sol calculee
302               CALL cltracrn(it, pdtphys, yu1, yv1, &               CALL cltracrn(it, pdtphys, yu1, yv1, coefh, t_seri, ftsol, &
303                    coefh, t_seri, ftsol, pctsrf, &                    pctsrf, tr_seri(:, :, it), trs(:, it), paprs, pplay, delp, &
304                    tr_seri(:, :, it), trs(1, it), &                    masktr(1, it), fshtr(1, it), hsoltr(it), tautr(it), &
305                    paprs, pplay, delp, &                    vdeptr(it), rlat, d_tr_cl(1, 1, it), d_trs)
                   masktr(1, it), fshtr(1, it), hsoltr(it), &  
                   tautr(it), vdeptr(it), &  
                   rlat, &  
                   d_tr_cl(1, 1, it), d_trs)  
306               DO k = 1, llm               DO k = 1, llm
307                  DO i = 1, klon                  DO i = 1, klon
308                     tr_seri(i, k, it) = tr_seri(i, k, it) + d_tr_cl(i, k, it)                     tr_seri(i, k, it) = tr_seri(i, k, it) + d_tr_cl(i, k, it)
309                  ENDDO                  ENDDO
310               ENDDO               ENDDO
311    
312               ! Traceur ds sol               trs(:, it) = trs(:, it) + d_trs
313              else
314               DO i = 1, klon               ! couche limite avec flux prescrit
                 trs(i, it) = trs(i, it) + d_trs(i)  
              END DO  
           else ! couche limite avec flux prescrit  
315               !MAF provisoire source / traceur a creer               !MAF provisoire source / traceur a creer
316               DO i=1, klon               DO i=1, klon
317                  source(i) = 0.0 ! pas de source, pour l'instant                  source(i) = 0. ! pas de source, pour l'instant
318               ENDDO               ENDDO
319    
320               CALL cltrac(pdtphys, coefh, t_seri, tr_seri(:, :, it), source, &               CALL cltrac(pdtphys, coefh, t_seri, tr_seri(:, :, it), source, &
Legend:
Removed from v.134  
changed lines
  Added in v.156
  ViewVC Help
Powered by ViewVC 1.1.21