anchor.diff:          real, dimension(nz),   target :: q1z,lz
           real, dimension(nz-1), target :: fz, U2z
           logical :: flux_is_w
  -
  +        
           real :: qmax_i,qmin_i,qmax_i2,qmin_i2
           real :: beta_in_i, beta_out_i, beta_in_i2, beta_out_i2
           real :: fin_i, fout_i, fin_i2, fout_i2
  -
  -        ! NOTE: before inclusion of FCT_core the following variables must be setup:
  +        
  +        ! NOTE: before inclusion of FCT_core the following variables must be setup: 
           ! q1 and l (l=q0)
           !$omp parallel shared(q2,q,u,v,w) firstprivate(nx,ny,nz) default(private)
           !$omp do schedule(static)
  
positive.code:         real,dimension(1:nx,1:nz,1:ny-1),intent(inout) :: v
          real,dimension(1:nx,1:nz,1:ny),  intent(inout) :: w
          integer, intent(in) :: nx,nz,ny
  
          integer :: i,j,k,n
          real, dimension(:), pointer :: q1, U2, l, f
          ! q1 = q after applying previous iteration advection
  
negative.code:                 w(k,1:n-1,j)=U2
                  w(k,n,j)=0
              end do
              
          end do
          !$omp end do
          
          flux_is_w=.False.
          n=ny
          q1=>q1y
  