anchor.diff:    return
   end subroutine
   
  -subroutine vel_update_vort
  +subroutine compute_omega
     use state
     use posvel
     use forces
  
positive.code:   end do
  
    return
  end subroutine vel_update_vort
  
  subroutine vorticity_force
    use state
  
negative.code: 
    real(WP), parameter :: ARTIFICIAL_VISCOSITY = 1.0e-2_WP
    real(WP), parameter :: VORTICITY_COEFFICIENT = 1.0e-4_WP
    logical,  parameter :: USE_VORTICITY_CONFINEMENT = .true.
  
    real(WP) :: PRESSURE_RADIUS_FACTOR
  
    real(WP),dimension(:), allocatable :: lm ! Lambda
  
  end module constants
  
  ! ===================================== !
  