18 |
|
|
19 |
! inputs: |
! inputs: |
20 |
integer, intent(in):: ncum, nd, na, nloc |
integer, intent(in):: ncum, nd, na, nloc |
21 |
integer icb(nloc), inb(nloc), nk(nloc) |
integer, intent(in):: icb(nloc), inb(nloc), nk(nloc) |
22 |
real sig(nloc, nd) |
real sig(nloc, nd) |
23 |
real t(nloc, nd), rr(nloc, nd), rs(nloc, nd) |
real t(nloc, nd), rr(nloc, nd), rs(nloc, nd) |
24 |
real u(nloc, nd), v(nloc, nd) |
real u(nloc, nd), v(nloc, nd) |
123 |
! *** if no air can entrain at level i assume that updraft detrains *** |
! *** if no air can entrain at level i assume that updraft detrains *** |
124 |
! *** at that level and calculate detrained air flux and properties *** |
! *** at that level and calculate detrained air flux and properties *** |
125 |
|
|
|
!@ do 170 i=icb(il), inb(il) |
|
|
|
|
126 |
do il=1, ncum |
do il=1, ncum |
127 |
if ((i >= icb(il)).and.(i <= inb(il)).and.(nent(il, i) == 0)) then |
if ((i >= icb(il)).and.(i <= inb(il)).and.(nent(il, i) == 0)) then |
128 |
!@ if(nent(il, i) == 0)then |
!@ if(nent(il, i) == 0)then |
140 |
! --- NORMALIZE ENTRAINED AIR MASS FLUXES |
! --- NORMALIZE ENTRAINED AIR MASS FLUXES |
141 |
! --- TO REPRESENT EQUAL PROBABILITIES OF MIXING |
! --- TO REPRESENT EQUAL PROBABILITIES OF MIXING |
142 |
|
|
143 |
call zilch(asum, nloc*nd) |
asum = 0. |
144 |
call zilch(csum, nloc*nd) |
csum = 0. |
|
call zilch(csum, nloc*nd) |
|
145 |
|
|
146 |
do il=1, ncum |
do il=1, ncum |
147 |
lwork(il) = .FALSE. |
lwork(il) = .FALSE. |