anchor.diff:    use precision
     implicit none
   
  -  real(WP) :: neigh_t, dens_t, pres_t, visc_t, ef_t, ps_t, tot_t
  +  real(WP) :: neigh_t, dens_t, force_t, ef_t, ps_t, tot_t
   
   end module timers
   
  
positive.code:   use posvel
    implicit none
  
    ! Only gravity for now
    fy = fy + GRAVITY
  
  
negative.code:   call cpu_time(t2)
    call compute_density
    call cpu_time(t3)
    call compute_forces
    call cpu_time(t5)
    call ext_forces
    call cpu_time(t6)
    call posvel_update
  