anchor.diff:                        + qr(k) * CVW(I_QR)
                endif
   
  -             fq(ij,kk,I_QV) = ( qv(k) - q(ij,kk,I_QV) ) / dt
  -             fe(ij,kk)      = ( cv(k) * t(ij,k) - ein(ij,kk) ) / dt
  +             fq(ij,kk,I_QV) = fq(ij,kk,I_QV) + ( qv(k) - q(ij,kk,I_QV) ) / dt
  +             fe(ij,kk)      = fe(ij,kk)      + ( cv(k) * t(ij,k) - ein(ij,kk) ) / dt
             enddo
          enddo
   
  
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
  