/[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 232 by guez, Tue Nov 7 10:23:25 2017 UTC revision 241 by guez, Mon Nov 13 11:51:04 2017 UTC
# Line 8  module phytrac_m Line 8  module phytrac_m
8  contains  contains
9    
10    SUBROUTINE phytrac(julien, gmtime, firstcal, lafin, pdtphys, t_seri, paprs, &    SUBROUTINE phytrac(julien, gmtime, firstcal, lafin, pdtphys, t_seri, paprs, &
11         pplay, pmfu, pmfd, pde_u, pen_d, coefh, fm_therm, entr_therm, yu1, &         pplay, pmfu, pmfd, pde_u, pen_d, coefh, cdragh, fm_therm, entr_therm, yu1, &
12         yv1, ftsol, pctsrf, frac_impa, frac_nucl, da, phi, mp, upwd, dnwd, &         yv1, ftsol, pctsrf, frac_impa, frac_nucl, da, phi, mp, upwd, dnwd, &
13         tr_seri, zmasse, ncid_startphy)         tr_seri, zmasse, ncid_startphy)
14    
# Line 77  contains Line 77  contains
77    
78      REAL pde_u(klon, llm) ! flux detraine dans le panache montant      REAL pde_u(klon, llm) ! flux detraine dans le panache montant
79      REAL pen_d(klon, llm) ! flux entraine dans le panache descendant      REAL pen_d(klon, llm) ! flux entraine dans le panache descendant
80      REAL coefh(:, :) ! (klon, llm) coeff melange couche limite      REAL coefh(:, 2:) ! (klon, 2:llm) coeff melange couche limite
81        real cdragh(:) ! (klon)
82      real fm_therm(klon, llm+1), entr_therm(klon, llm) ! thermiques      real fm_therm(klon, llm+1), entr_therm(klon, llm) ! thermiques
83      REAL, intent(in):: yu1(:), yv1(:) ! (klon) vent au premier niveau      REAL, intent(in):: yu1(:), yv1(:) ! (klon) vent au premier niveau
84    
# Line 279  contains Line 280  contains
280      DO it=1, nqmx - 2      DO it=1, nqmx - 2
281         if (clsol(it)) then         if (clsol(it)) then
282            ! couche limite avec quantite dans le sol calculee            ! couche limite avec quantite dans le sol calculee
283            CALL cltracrn(it, pdtphys, yu1, yv1, coefh(:, 2:llm), coefh(:, 1), &            CALL cltracrn(it, pdtphys, yu1, yv1, coefh, cdragh, t_seri, ftsol, &
284                 t_seri, ftsol, pctsrf, tr_seri(:, :, it), trs(:, it), paprs, &                 pctsrf, tr_seri(:, :, it), trs(:, it), paprs, pplay, delp, &
285                 pplay, delp, masktr(1, it), fshtr(1, it), hsoltr(it), &                 masktr(1, it), fshtr(1, it), hsoltr(it), tautr(it), &
286                 tautr(it), vdeptr(it), rlat, d_tr_cl(1, 1, it), d_trs)                 vdeptr(it), rlat, d_tr_cl(1, 1, it), d_trs)
287            DO k = 1, llm            DO k = 1, llm
288               DO i = 1, klon               DO i = 1, klon
289                  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)
# Line 297  contains Line 298  contains
298               source(i) = 0. ! pas de source, pour l'instant               source(i) = 0. ! pas de source, pour l'instant
299            ENDDO            ENDDO
300    
301            CALL cltrac(pdtphys, coefh(:, 2:llm), t_seri, tr_seri(:, :, it), source, &            CALL cltrac(pdtphys, coefh, t_seri, tr_seri(:, :, it), source, &
302                 paprs, pplay, delp, d_tr_cl(1, 1, it))                 paprs, pplay, delp, d_tr_cl(1, 1, it))
303            DO k = 1, llm            DO k = 1, llm
304               DO i = 1, klon               DO i = 1, klon
Legend:
Removed from v.232  
changed lines
  Added in v.241
  ViewVC Help
Powered by ViewVC 1.1.21