anchor.diff:  module convection
       use data_structures
       use module_cu_tiedtke,  only: tiedtkeinit, CU_TIEDTKE
  -    use string,             only: str
  -    use io_routines,        only: io_write3d, io_write2d
       
       implicit none
  -    integer :: output_number
  -    real,allocatable,dimension(:,:,:)::p8,U3D,V3D,T3D, & 
  -                                       RTHCUTEN,RQVCUTEN,RQCCUTEN,RQICUTEN,    &
  -                                       RUCUTEN,RVCUTEN
  -    logical,allocatable,dimension(:,:)::CU_ACT_FLAG
  -    real,allocatable,dimension(:,:) ::XLAND, RAINCV, PRATEC
  -    real,allocatable,dimension(:) :: ZNU ! = (/0.9975, 0.9915, 0.984, 0.975, 0.965, 0.9525, 0.9375, &
  -                                        !    0.92, 0.9, 0.88, 0.86, 0.835, 0.805, 0.775, 0.745,   &
  -                                        !    0.71, 0.67, 0.63, 0.59, 0.55, 0.51, 0.47, 0.43, 0.39,&
  -                                        !    0.355, 0.325, 0.295, 0.27, 0.25, 0.23, 0.21, 0.19,   &
  -                                        !    0.17, 0.15, 0.13, 0.11, 0.09100001, 0.074, 0.059,    &
  -                                        !    0.046, 0.035, 0.025, 0.015, 0.005/)
  +    real,   allocatable, dimension(:,:,:) :: p8,U3D,V3D,T3D, & 
  +                                             RTHCUTEN,RQVCUTEN,RQCCUTEN,RQICUTEN,    &
  +                                             RUCUTEN,RVCUTEN
  +    logical,allocatable, dimension(:,:) :: CU_ACT_FLAG
  +    real,   allocatable, dimension(:,:) :: XLAND, RAINCV, PRATEC
  +    real,   allocatable, dimension(:)   :: ZNU
                                           
   contains
       subroutine znu_init(domain)
           type(domain_type), intent(in) :: domain
  -        integer::n_levels,i,xpt,ypt
  -        real::ptop,psfc
  +        integer :: n_levels,i,xpt,ypt
  +        real    :: ptop,psfc
           
           n_levels=size(domain%p,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:         do i=1,n_levels
              znu(i)=(domain%p(xpt,i,ypt)-ptop)/(psfc-ptop)
          enddo
      end subroutine znu_init
  
      subroutine init_convection(domain,options)
          implicit none
          type(domain_type),  intent(in) :: domain
          type(options_type), intent(in) :: options
          integer :: ids,ide,jds,jde,kds,kde
          
          write(*,*) "Initializing Cumulus Scheme"
          ids=1
  