da_autoc_lin.inc

1 subroutine da_autoc_lin(DT,SCR2,PRC,SCR29,PRC9,kts,kte,kms,kme)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: Auto-conversion
5    !-----------------------------------------------------------------------
6 
7    implicit none
8 
9    integer, intent(in)                   :: kts, kte, kms, kme
10    real, dimension(kms:kme), intent(in)  :: DT, SCR2, SCR29
11    real, dimension(kms:kme), intent(out) :: PRC, PRC9
12 
13    integer      :: k
14    real         :: QCTH, alpha
15 
16    QCTH = 0.5E-3
17    alpha = 0.001
18 
19    do k = kts, kte
20 
21       if ( DT(k) <= 0. ) cycle
22 
23       if ( SCR2(k) >= QCTH ) then
24          PRC9(K) = alpha * SCR29(K)
25          PRC(k) = alpha * ( SCR2(k) - QCTH )
26       else
27          PRC9(K) = 0.
28          PRC(k) = 0.
29       end if
30    end do
31 
32 end subroutine da_autoc_lin
33 
34