da_accre.inc

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