16 |
! Vertical profile of buoyancy computed here (use of buoybase). |
! Vertical profile of buoyancy computed here (use of buoybase). |
17 |
|
|
18 |
use conema3_m, only: epmax |
use conema3_m, only: epmax |
19 |
use cv30_param_m, only: dtovsh, minorig, nl, nlp, pbcrit, ptcrit, spfac |
use cv30_param_m, only: dtovsh, minorig, nl, pbcrit, ptcrit, spfac |
20 |
use cvthermo, only: cl, clmcpv, cpd, cpv, eps, lv0, rrv |
use cvthermo, only: cl, clmcpv, cpd, cpv, eps, lv0, rrv |
21 |
|
|
22 |
! inputs: |
! inputs: |
137 |
! tvp est calcule en une seule fois, et sans retirer |
! tvp est calcule en une seule fois, et sans retirer |
138 |
! l'eau condensee (~> reversible CAPE) |
! l'eau condensee (~> reversible CAPE) |
139 |
do i = 1, ncum |
do i = 1, ncum |
140 |
tp(i, nlp) = tp(i, nl) |
tp(i, nl + 1) = tp(i, nl) |
141 |
end do |
end do |
142 |
|
|
143 |
! EFFECTIVE VERTICAL PROFILE OF BUOYANCY: |
! EFFECTIVE VERTICAL PROFILE OF BUOYANCY: |
177 |
|
|
178 |
! CALCULATE LIQUID WATER STATIC ENERGY OF LIFTED PARCEL |
! CALCULATE LIQUID WATER STATIC ENERGY OF LIFTED PARCEL |
179 |
|
|
180 |
do k = 1, nlp |
do k = 1, nl + 1 |
181 |
do i = 1, ncum |
do i = 1, ncum |
182 |
hp(i, k) = h(i, k) |
hp(i, k) = h(i, k) |
183 |
enddo |
enddo |