/[lmdze]/trunk/Sources/phylmd/CV30_routines/cv30_unsat.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/CV30_routines/cv30_unsat.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 193 by guez, Thu May 12 13:22:19 2016 UTC revision 194 by guez, Thu May 12 14:35:35 2016 UTC
# Line 41  contains Line 41  contains
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
# Line 86  contains Line 86  contains
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)
Legend:
Removed from v.193  
changed lines
  Added in v.194
  ViewVC Help
Powered by ViewVC 1.1.21