anchor.diff:  !         write(*,*) MAXVAL(RQICUTEN(2:ide-1,:,2:jde-1)), MINVAL(RQICUTEN(2:ide-1,:,2:jde-1))
   !         write(*,*) "rain range"
   !         write(*,*) MAXVAL(RAINCV(2:ide-1,2:jde-1)), MINVAL(RAINCV(2:ide-1,2:jde-1))
  -        domain%qv=domain%qv+RQVCUTEN*dt_in
  -        domain%cloud=domain%cloud+RQCCUTEN*dt_in
  -        domain%th=domain%th+RTHCUTEN*dt_in
  -        domain%ice=domain%ice+RQICUTEN*dt_in
  +
  +! COMMENTED OUT BECAUSE IT WAS REMOVING MORE WATER THAN IT WAS FORMING PRECIPITATION
  +!         domain%qv=domain%qv+RQVCUTEN*dt_in
  +! COMMENTED OUT BECAUSE IT WAS ALWAYS 0 (AT LEAST IN A SHORT RUN)
  +!         domain%cloud=domain%cloud+RQCCUTEN*dt_in
  +! COMMENTED OUT TO BE CONSISTENT WITH REMOVING QVTEN
  +!         domain%th=domain%th+RTHCUTEN*dt_in
  +! COMMENTED OUT BECAUSE IT WAS ALWAYS 0 (AT LEAST IN A SHORT RUN)
  +!         domain%ice=domain%ice+RQICUTEN*dt_in
           domain%rain=domain%rain+RAINCV
           domain%crain=domain%crain+RAINCV
  +        
  +!         call io_write3d("convtest/cu_rho_"//trim(str(output_number))//".nc","data",domain%rho)
  +!         call io_write3d("convtest/cuten_qv_"//trim(str(output_number))//".nc","data",RQVCUTEN*dt_in)
  +!         call io_write3d("convtest/cuten_qc_"//trim(str(output_number))//".nc","data",RQCCUTEN*dt_in)
  +!         call io_write3d("convtest/cuten_th_"//trim(str(output_number))//".nc","data",RTHCUTEN*dt_in)
  +!         call io_write3d("convtest/cuten_qi_"//trim(str(output_number))//".nc","data",RQICUTEN*dt_in)
  +!         call io_write2d("convtest/cuten_cr_"//trim(str(output_number))//".nc","data",RAINCV)
  +!         output_number=output_number+1
  +        
   !       write(*,*) MAXVAL(domain%qv(2:ide-1,:,2:jde-1)), MINVAL(domain%qv(2:ide-1,:,2:jde-1))
       endif
       
  
positive.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)
  
negative.code: !! ----------------------------------------------------------------------------
  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
  