da_accre_lin.inc

1 subroutine da_accre_lin(DT,SCR2,SCR3,PRA,SCR29,SCR39,PRA9,kts,kte,kms,kme)
2 
3 
4    !-----------------------------------------------------------------------
5    ! Purpose: Accretion (Collection)
6    !-----------------------------------------------------------------------
7 
8    implicit none
9 
10    integer, intent(in)                   :: kts, kte, kms, kme
11    real, dimension(kms:kme), intent(in)  :: DT, SCR2, SCR3, SCR29, SCR39
12    real, dimension(kms:kme), intent(out) :: PRA, PRA9
13 
14    integer   :: k
15    real      :: gamma, QRTH
16 
17    gamma = 0.002
18    QRTH = 1.e-6
19 
20    do K = kts, kte
21 
22       if ( DT(k) <= 0. ) cycle
23 
24       if (SCR2(k) > 0. .and. SCR3(k) > QRTH ) then
25          PRA9(K) = gamma * 0.875 * SCR2(k) * (SCR3(K)*1.e3)**(-0.125) * 1.e3 * SCR39(K)  &
26                       + gamma * SCR29(k) * (SCR3(K)*1.e3)**0.875
27          PRA(k) = gamma * SCR2(k) * (SCR3(k)*1.e3)**0.875
28       else if (SCR2(k) > 0. .and. SCR3(k) <= QRTH ) then
29          PRA9(K) = gamma * SCR29(k) * (QRTH*1.e3)**0.875
30          PRA(k) = gamma * SCR2(k) * (QRTH*1.e3)**0.875
31       else
32          PRA9(K) = 0.
33          PRA(k) = 0.
34       end if
35 
36    end do
37 
38 end subroutine da_accre_lin