da_autoc.inc

1 subroutine da_autoc(DT,SCR2,PRC,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
11    real, dimension(kms:kme), intent(out) :: PRC
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          PRC(k) = alpha * ( SCR2(k) - QCTH )
25       else
26          PRC(k) = 0.
27       end if
28    end do
29 
30 end subroutine da_autoc
31