anchor.diff:          !$omp default(shared)
           !$omp do schedule(static)
           do j=jds,jde
  -            RAINCV(:,j)=0
  -            p8(:,:kde-1,j)=(domain%p(:,kds:kde-1,j)+domain%p(:,kds+1:kde,j))/2
  -            p8(:,kde,j)=p8(:,kde-1,j)+(p8(:,kde-1,j)-p8(:,kde-2,j))/2
  +            RAINCV(:,j)     = 0
  +            RQVCUTEN(:,:,j) = 0
  +            RTHCUTEN(:,:,j) = 0
  +            RQCCUTEN(:,:,j) = 0
  +            RQICUTEN(:,:,j) = 0
  +            p8(:,:kde-1,j)  = (domain%p(:,kds:kde-1,j) + domain%p(:,kds+1:kde,j))/2
  +            p8(:,kde,j)     =  domain%p(:,kde,j) + (domain%p(:,kde,j)-p8(:,kde-1,j))
               
               V3D(:,:,j)=(domain%v(:,:,j+1)+domain%v(:,:,j))/2
               U3D(:,:,j)=(domain%u(ids:ide,:,j)+domain%u(ids+1:ide+1,:,j))/2
  
positive.code: 
          call CU_TIEDTKE(                                          &
                   dt_in,itimestep,STEPCU                           &
                  ,RAINCV,PRATEC,domain%latent_heat/LH_vaporization,domain%sensible_heat,ZNU(kds:kde) &
                  ,U3D,V3D,domain%w*domain%dz/domain%dx,T3D,domain%qv,domain%cloud,domain%ice,domain%pii,domain%rho &
                  ,domain%qv_adv_tendency,domain%qv_pbl_tendency    &
                  ,domain%dz,p8,domain%p,XLAND,CU_ACT_FLAG          &
                  ,ids+1,ide-1, jds+1,jde-1, kds,kde-1              &
  
negative.code:     
  subroutine convect(domain,options,dt_in)
      implicit none
      type(domain_type),  intent(inout) :: domain
      type(options_type), intent(in)    :: options
      real, intent(in) :: dt_in
      
      integer :: ids,ide,jds,jde,kds,kde,j,itimestep,STEPCU
      
      itimestep=1
      STEPCU=1
  