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(:, :) ! (klon, llm) coeff melange couche limite |
81 |
real fm_therm(klon, llm+1), entr_therm(klon, llm) ! thermiques |
real fm_therm(klon, llm+1), entr_therm(klon, llm) ! thermiques |
82 |
REAL, intent(in):: yu1(klon), yv1(klon) ! vent au premier niveau |
REAL, intent(in):: yu1(:), yv1(:) ! (klon) vent au premier niveau |
83 |
|
|
84 |
! Arguments n\'ecessaires pour les sources et puits de traceur : |
! Arguments n\'ecessaires pour les sources et puits de traceur : |
85 |
real, intent(in):: ftsol(:, :) ! (klon, nbsrf) surface temperature (K) |
real, intent(in):: ftsol(:, :) ! (klon, nbsrf) surface temperature (K) |
279 |
DO it=1, nqmx - 2 |
DO it=1, nqmx - 2 |
280 |
if (clsol(it)) then |
if (clsol(it)) then |
281 |
! couche limite avec quantite dans le sol calculee |
! couche limite avec quantite dans le sol calculee |
282 |
CALL cltracrn(it, pdtphys, yu1, yv1, coefh, t_seri, ftsol, & |
CALL cltracrn(it, pdtphys, yu1, yv1, coefh(:, 2:llm), coefh(:, 1), & |
283 |
pctsrf, tr_seri(:, :, it), trs(:, it), paprs, pplay, delp, & |
t_seri, ftsol, pctsrf, tr_seri(:, :, it), trs(:, it), paprs, & |
284 |
masktr(1, it), fshtr(1, it), hsoltr(it), tautr(it), & |
pplay, delp, masktr(1, it), fshtr(1, it), hsoltr(it), & |
285 |
vdeptr(it), rlat, d_tr_cl(1, 1, it), d_trs) |
tautr(it), vdeptr(it), rlat, d_tr_cl(1, 1, it), d_trs) |
286 |
DO k = 1, llm |
DO k = 1, llm |
287 |
DO i = 1, klon |
DO i = 1, klon |
288 |
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) |
297 |
source(i) = 0. ! pas de source, pour l'instant |
source(i) = 0. ! pas de source, pour l'instant |
298 |
ENDDO |
ENDDO |
299 |
|
|
300 |
CALL cltrac(pdtphys, coefh, t_seri, tr_seri(:, :, it), source, & |
CALL cltrac(pdtphys, coefh(:, 2:llm), t_seri, tr_seri(:, :, it), source, & |
301 |
paprs, pplay, delp, d_tr_cl(1, 1, it)) |
paprs, pplay, delp, d_tr_cl(1, 1, it)) |
302 |
DO k = 1, llm |
DO k = 1, llm |
303 |
DO i = 1, klon |
DO i = 1, klon |