anchor.diff:    do subiter = 1, ITERATIONS_PER_STEP
       ! Compute lamnda
       call compute_lambda
  -  print*,'after lambda', maxval(cvx), maxval(cvy), maxval(cvz)
  +!  print*,'after lambda', maxval(cvx), maxval(cvy), maxval(cvz)
       ! Jiggle the particles
       call particle_jiggle
  -  print*,'after jiggle', maxval(cvx), maxval(cvy), maxval(cvz)
  -  end do
  +!  print*,'after jiggle', maxval(cvx), maxval(cvy), maxval(cvz)
  +  end do
  +  call update_candidate_vel
   
     ! Confine vorticity, update vel, and calculate force
     if(USE_VORTICITY_CONFINEMENT) then
  
positive.code: end subroutine deinit
  
  
  ! ===================================== !
  ! One step of SPH                       !
  ! ===================================== !
  
negative.code:       dvz = vz(nb) - vz(i)
  
        k = viscositykernel(i,nb)
        r = (cpx(i)-cpx(nb))**2+(cpy(i)-cpy(nb))**2+(cpz(i)-cpz(nb))**2
        r = sqrt(r)
        print*,r-KERNEL_SIZE,k, i, nb, k*dvx, dvx, vx(nb), vx(i)
        sx = sx + k * dvx
        sy = sy + k * dvy
        sz = sz + k * dvz
  