anchor.diff:  
     nbs = 0
     nc = 0
  +
     do i = 1, n
       binx = NINT((px(i))/(BOX_SIZE)*(real(nbinx,WP)-1.0_WP)) + 1
       biny = NINT((py(i))/(BOX_SIZE)*(real(nbiny,WP)-1.0_WP)) + 1
  
positive.code:   call deinit
    call cpu_time(deinit_stop)
  
    write(*,"(A)")                 "Timing Information:                Total Time              Time Per Step"
    write(*,"(A)")                 "---------------------             ------------            ---------------"
    write(*,"(A,ES25.11)")         "Initialization:        ",init_stop-init_start
  
negative.code:   real(WP), parameter :: ARTIFICIAL_PRESSURE_STRENGTH = 1.0e-5_WP
    real(WP)  :: ARTIFICIAL_PRESSURE_RADIUS, PRESSURE_RADIUS_FACTOR
    logical, parameter :: USE_SURFACE_TENSION = .false.
    real(WP), parameter :: STC = 0.0078
  
  
  end module constants
  