41 |
|
|
42 |
! Local: |
! Local: |
43 |
integer ncum |
integer ncum |
44 |
integer i, j, il, imax |
integer i, il, imax |
45 |
real tinv, delti |
real tinv, delti |
46 |
real awat, afac, afac1, afac2, bfac |
real afac, afac1, afac2, bfac |
47 |
real pr1, pr2, sigt, b6, c6, revap, tevap, delth |
real pr1, pr2, sigt, b6, c6, revap, tevap, delth |
48 |
real amfac, amp2, xf, tf, fac2, ur, sru, fac, d, af, bf |
real amfac, amp2, xf, tf, fac2, ur, sru, fac, d, af, bf |
49 |
real ampmax |
real ampmax |
86 |
! and condensed water flux |
! and condensed water flux |
87 |
|
|
88 |
! Calculate detrained precipitation |
! Calculate detrained precipitation |
89 |
|
forall (il = 1:ncum, inb(il) >= i .and. lwork(il)) wdtrain(il) = grav & |
90 |
do il = 1, ncum |
* (ep(il, i) * m(il, i) * clw(il, i) & |
91 |
if (i <= inb(il) .and. lwork(il)) then |
+ sum(max(elij(il, :i - 1, i) - (1. - ep(il, i)) * clw(il, i), 0.) & |
92 |
wdtrain(il) = grav * ep(il, i) * m(il, i) * clw(il, i) |
* ment(il, :i - 1, i))) |
|
endif |
|
|
enddo |
|
|
|
|
|
if (i > 1) then |
|
|
do j = 1, i - 1 |
|
|
do il = 1, ncum |
|
|
if (i <= inb(il) .and. lwork(il)) then |
|
|
awat = elij(il, j, i) - (1. - ep(il, i)) * clw(il, i) |
|
|
awat = max(awat, 0.) |
|
|
wdtrain(il) = wdtrain(il) + grav * awat * ment(il, j, i) |
|
|
endif |
|
|
enddo |
|
|
end do |
|
|
endif |
|
93 |
|
|
94 |
! Find rain water and evaporation using provisional |
! Find rain water and evaporation using provisional |
95 |
! estimates of qp(i) and qp(i - 1) |
! estimates of qp(i) and qp(i - 1) |