anchor.diff:                  + qc(:) * CVW(I_QC) &
                   + qr(:) * CVW(I_QR)
   
  -          fq(ij,kmin:kmax,I_QV) = ( qv(:) - q(ij,kmin:kmax,I_QV) ) / dt
  -          fq(ij,kmin:kmax,I_QC) = ( qc(:) - q(ij,kmin:kmax,I_QC) ) / dt
  -          fq(ij,kmin:kmax,I_QR) = ( qr(:) - q(ij,kmin:kmax,I_QR) ) / dt
  -          fe(ij,kmin:kmax)      = ( cv(:) * theta(:) * pk(:) - ein(ij,kmin:kmax) ) / dt
  +          fq(ij,kmin:kmax,I_QV) = fq(ij,kmin:kmax,I_QV) + ( qv(:) - q(ij,kmin:kmax,I_QV) ) / dt
  +          fq(ij,kmin:kmax,I_QC) = fq(ij,kmin:kmax,I_QC) + ( qc(:) - q(ij,kmin:kmax,I_QC) ) / dt
  +          fq(ij,kmin:kmax,I_QR) = fq(ij,kmin:kmax,I_QR) + ( qr(:) - q(ij,kmin:kmax,I_QR) ) / dt
  +          fe(ij,kmin:kmax)      = fe(ij,kmin:kmax)      + ( cv(:) * theta(:) * pk(:) - ein(ij,kmin:kmax) ) / dt
          enddo
       endif
   
  
positive.code:        do k = 1, vlayer
            kk = kmax - k + 1 ! reverse order
  
            fvx(:,kk) = ( u(:,k) * ix(:) + v(:,k) * jx(:) - vx(:,kk) ) / dt
            fvy(:,kk) = ( u(:,k) * iy(:) + v(:,k) * jy(:) - vy(:,kk) ) / dt
            fvz(:,kk) = ( u(:,k) * iz(:) + v(:,k) * jz(:) - vz(:,kk) ) / dt
         enddo
  
         do k = 1, vlayer
  
negative.code: 
            call tendency_Terminator( lat_deg, lon_deg, cl, cl2, dt, cl_f, cl2_f )
  
            fq(ij,k,NCHEM_STR) = fq(ij,k,NCHEM_STR) + cl_f
            fq(ij,k,NCHEM_END) = fq(ij,k,NCHEM_END) + cl2_f
         enddo
         enddo
      endif
  