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 |