da_autoc_adj.inc

1 subroutine da_autoc_adj(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
11    real, dimension(kms:kme), intent(in)    :: PRC, PRC9
12    real, dimension(kms:kme), intent(inout) :: SCR29
13 
14    integer      :: k
15    real         :: QCTH, alpha
16 
17    QCTH = 0.5E-3
18    alpha = 0.001
19 
20    do k = kts, kte
21       if ( DT(k) <= 0. ) cycle
22 
23       if ( SCR2(k) >= QCTH ) then
24          SCR29(K) = SCR29(K) + alpha * PRC9(K)
25       end if
26    end do
27 
28 end subroutine da_autoc_adj