anchor.diff:    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
  
positive.code:   write(*,"(A,ES25.11,ES25.11)") "Ext. Force Calculation:",ef_t,ef_t/iter
    write(*,"(A,ES25.11,ES25.11)") "Posvel Calculation:    ",ps_t,ps_t/iter
    write(*,"(A,ES25.11)")         "Deinitialization Time: ",deinit_stop-deinit_start
    total_time = init_stop-init_start+neigh_t+dens_t &
                 +pres_t+visc_t+ef_t+ps_t+deinit_stop-deinit_start
    write(*,"(A,ES25.11,ES25.11)") "Total Time:            ",total_time,total_time/iter
  
  end program sph_run
  
negative.code: 
    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
  