74 |
firstcal = .FALSE. |
firstcal = .FALSE. |
75 |
ENDIF |
ENDIF |
76 |
|
|
77 |
DO i = 1, klon |
ldcum = .FALSE. |
78 |
ldcum(i) = .FALSE. |
dt_con = 0. |
79 |
ENDDO |
dq_con = 0. |
|
DO k = 1, klev |
|
|
DO i = 1, klon |
|
|
dt_con(i, k) = 0.0 |
|
|
dq_con(i, k) = 0.0 |
|
|
ENDDO |
|
|
ENDDO |
|
|
|
|
|
! initialiser les variables et faire l'interpolation verticale |
|
80 |
|
|
81 |
|
! Initialiser les variables et faire l'interpolation verticale : |
82 |
CALL flxini(ten, qen, qsen, pgeo, paph, zgeoh, ztenh, zqenh, zqsenh, & |
CALL flxini(ten, qen, qsen, pgeo, paph, zgeoh, ztenh, zqenh, zqsenh, & |
83 |
ptu, pqu, ptd, pqd, mfd, zmfds, zmfdq, zdmfdp, mfu, mfus, mfuq, & |
ptu, pqu, ptd, pqd, mfd, zmfds, zmfdq, zdmfdp, mfu, mfus, mfuq, & |
84 |
zdmfup, zdpmel, plu, plude, ilab, pen_u, pde_u, pen_d, pde_d) |
zdmfup, zdpmel, plu, plude, ilab, pen_u, pde_u, pen_d, pde_d) |
85 |
|
|
86 |
! determiner les valeurs au niveau de base de la tour convective |
! Déterminer les valeurs au niveau de base de la tour convective : |
87 |
|
CALL flxbase(ztenh, zqenh, zgeoh, paph, ptu, pqu, plu, ldcum, kcbot, ilab) |
|
CALL flxbase(ztenh, zqenh, zgeoh, paph, & |
|
|
ptu, pqu, plu, ldcum, kcbot, ilab) |
|
|
|
|
|
! calculer la convergence totale de l'humidite et celle en provenance |
|
|
! de la couche limite, plus precisement, la convergence integree entre |
|
|
! le sol et la base de la convection. Cette derniere convergence est |
|
|
! comparee avec l'evaporation obtenue dans la couche limite pour |
|
|
! determiner le type de la convection |
|
88 |
|
|
89 |
k=1 |
! Calculer la convergence totale de l'humidité et celle en |
90 |
DO i = 1, klon |
! provenance de la couche limite, plus précisément, la convergence |
91 |
zdqcv(i) = pqte(i, k)*(paph(i, k+1)-paph(i, k)) |
! intégrée entre le sol et la base de la convection. Cette |
92 |
zdhpbl(i) = 0.0 |
! dernière convergence est comparée avec l'&vaporation obtenue |
93 |
zdqpbl(i) = 0.0 |
! dans la couche limite pour déterminer le type de la convection. |
94 |
ENDDO |
|
95 |
|
zdqcv = pqte(:, 1) * (paph(:, 2) - paph(:, 1)) |
96 |
|
zdhpbl = 0. |
97 |
|
zdqpbl = 0. |
98 |
|
|
99 |
DO k=2, klev |
DO k=2, klev |
100 |
DO i = 1, klon |
DO i = 1, klon |
182 |
IF(ktype(i) == 2) zentr(i)=ENTRSCV |
IF(ktype(i) == 2) zentr(i)=ENTRSCV |
183 |
ENDDO |
ENDDO |
184 |
|
|
185 |
IF (lmfdd) THEN ! si l'on considere le panache descendant |
downdraft: IF (lmfdd) THEN |
186 |
|
! si l'on considere le panache descendant |
187 |
! calculer la precipitation issue du panache ascendant pour |
! calculer la precipitation issue du panache ascendant pour |
188 |
! determiner l'existence du panache descendant dans la convection |
! determiner l'existence du panache descendant dans la convection |
189 |
DO i = 1, klon |
DO i = 1, klon |
255 |
DO i = 1, klon |
DO i = 1, klon |
256 |
IF (lddraf(i)) mfub(i)=mfub1(i) |
IF (lddraf(i)) mfub(i)=mfub1(i) |
257 |
ENDDO |
ENDDO |
258 |
ENDIF ! fin de test sur lmfdd |
ENDIF downdraft |
259 |
|
|
260 |
! calculer de nouveau le panache ascendant |
! calculer de nouveau le panache ascendant |
261 |
|
|