4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
SUBROUTINE conflx(dtime, pres_h, pres_f, t, q, con_t, con_q, qhfl, omega, & |
SUBROUTINE conflx(pres_h, pres_f, t, q, con_t, con_q, qhfl, omega, d_t, d_q, & |
8 |
d_t, d_q, rain, snow, mfu, mfd, pen_u, pde_u, pen_d, pde_d, kcbot, & |
rain, snow, mfu, mfd, pen_u, pde_u, pen_d, pde_d, kcbot, kctop, kdtop, & |
9 |
kctop, kdtop, pmflxr, pmflxs) |
pmflxr, pmflxs) |
10 |
|
|
11 |
! From LMDZ4/libf/phylmd/conflx.F, version 1.1.1.1 2004/05/19 12:53:08 |
! From LMDZ4/libf/phylmd/conflx.F, version 1.1.1.1 2004/05/19 12:53:08 |
12 |
|
|
20 |
! par Olivier Boucher et Alexandre Armengaud pour le mélange et le |
! par Olivier Boucher et Alexandre Armengaud pour le mélange et le |
21 |
! lessivage des traceurs passifs. |
! lessivage des traceurs passifs. |
22 |
|
|
23 |
use flxmain_m, only: flxmain |
use comconst, only: dtphys |
24 |
USE dimphy, ONLY: klev, klon |
USE dimphy, ONLY: klev, klon |
25 |
|
USE fcttre, ONLY: foeew |
26 |
|
use flxmain_m, only: flxmain |
27 |
USE suphec_m, ONLY: rd, retv, rtt |
USE suphec_m, ONLY: rd, retv, rtt |
28 |
USE yoethf_m, ONLY: r2es |
USE yoethf_m, ONLY: r2es |
|
USE fcttre, ONLY: foeew |
|
29 |
|
|
|
REAL, intent(in):: dtime ! pas d'integration (s) |
|
30 |
REAL, intent(in):: pres_h(:, :) ! (klon, klev + 1) pression half-level (Pa) |
REAL, intent(in):: pres_h(:, :) ! (klon, klev + 1) pression half-level (Pa) |
31 |
REAL, intent(in):: pres_f(:, :) ! (klon, klev) pression full-level (Pa) |
REAL, intent(in):: pres_f(:, :) ! (klon, klev) pression full-level (Pa) |
32 |
REAL, intent(in):: t(:, :) ! (klon, klev) temperature (K) |
REAL, intent(in):: t(:, :) ! (klon, klev) temperature (K) |
137 |
ENDDO |
ENDDO |
138 |
|
|
139 |
! Appeler la routine principale : |
! Appeler la routine principale : |
140 |
CALL flxmain(dtime, t, q, qsen, qhfl, paprsf, paprs, zgeom, land, & |
CALL flxmain(dtphys, t, q, qsen, qhfl, paprsf, paprs, zgeom, land, & |
141 |
zcvgt, zcvgq, pvervel, rain, snow, kcbot, kctop, kdtop, mfu, mfd, & |
zcvgt, zcvgq, pvervel, rain, snow, kcbot, kctop, kdtop, mfu, mfd, & |
142 |
zen_u, zde_u, zen_d, zde_d, d_t_bis, d_q_bis, zmflxr, zmflxs) |
zen_u, zde_u, zen_d, zde_d, d_t_bis, d_q_bis, zmflxr, zmflxs) |
143 |
|
|
146 |
! la convection des traceurs. |
! la convection des traceurs. |
147 |
DO k = 1, klev |
DO k = 1, klev |
148 |
DO i = 1, klon |
DO i = 1, klon |
149 |
d_q(i, klev + 1-k) = dtime*d_q_bis(i, k) |
d_q(i, klev + 1-k) = dtphys*d_q_bis(i, k) |
150 |
d_t(i, klev + 1-k) = dtime*d_t_bis(i, k) |
d_t(i, klev + 1-k) = dtphys*d_t_bis(i, k) |
151 |
ENDDO |
ENDDO |
152 |
ENDDO |
ENDDO |
153 |
|
|