anchor.diff:  subroutine computePressureRadiusFactor
     use constants
     implicit none
  -  real(WP) :: r
  +  real(WP) :: r, c
   
     r = ARTIFICIAL_PRESSURE_RADIUS
  +  c = 315.0_WP / (64.0_WP * PI)
   
     PRESSURE_RADIUS_FACTOR = 1.0_WP / &
       (c * (KERNEL_SIZE*KERNEL_SIZE - r*r)**3 &
  
positive.code:       use posvel
        implicit none
        integer, intent(in) :: i,j
        real(WP) :: dx, dy, dz, r
  
        dx = cpx(i) - cpx(j)
        dy = cpy(i) - cpy(j)
        dz = cpz(i) - cpz(j)
        r = sqrt(dx*dx + dy*dy + dz*dz)
  
        if(r .gt. KERNEL_SIZE) then
          poly6kernel = 0.0_WP
  
negative.code: ! ===================================== !
  module neighbors
    use precision
    use constants
    implicit none
    integer,dimension(:),allocatable :: nc
    integer,dimension(:,:),allocatable :: nbs
    integer, parameter :: num_of_bins = ceiling(2*BOX_SIZE/KERNEL_SIZE) ! Bins in each direction
  end module neighbors
  
  ! ===================================== !
  