17 |
|
|
18 |
use conema3_m, only: epmax |
use conema3_m, only: epmax |
19 |
use cv30_param_m, only: dtovsh, minorig, nl, 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 cv_thermo_m, only: cl, clmcpv, cpd, cpv, eps, lv0, rrv |
21 |
USE dimphy, ONLY: klon, klev |
USE dimphy, ONLY: klon, klev |
22 |
|
|
23 |
integer, intent(in):: ncum |
integer, intent(in):: ncum |
35 |
! outputs: |
! outputs: |
36 |
integer, intent(out):: inb(:) ! (ncum) |
integer, intent(out):: inb(:) ! (ncum) |
37 |
! first model level above the level of neutral buoyancy of the |
! first model level above the level of neutral buoyancy of the |
38 |
! parcel (<= nl - 1) |
! parcel (1 <= inb <= nl - 1) |
39 |
|
|
40 |
real tp(klon, klev), tvp(klon, klev), clw(klon, klev) |
real tp(klon, klev), tvp(klon, klev), clw(klon, klev) |
41 |
! condensed water not removed from tvp |
! condensed water not removed from tvp |
131 |
do i = 1, ncum |
do i = 1, ncum |
132 |
pden = ptcrit - pbcrit |
pden = ptcrit - pbcrit |
133 |
ep(i, k) = (plcl(i) - p(i, k) - pbcrit) / pden * epmax |
ep(i, k) = (plcl(i) - p(i, k) - pbcrit) / pden * epmax |
134 |
ep(i, k) = amax1(ep(i, k), 0.0) |
ep(i, k) = max(ep(i, k), 0.0) |
135 |
ep(i, k) = amin1(ep(i, k), epmax) |
ep(i, k) = min(ep(i, k), epmax) |
136 |
sigp(i, k) = spfac |
sigp(i, k) = spfac |
137 |
end do |
end do |
138 |
end do |
end do |