anchor.diff:    call h5md_write_dataset(timers_group, solvent%time_sort%name, solvent%time_sort%total)
     call h5md_write_dataset(timers_group, solvent%time_ct%name, solvent%time_ct%total)
     call h5md_write_dataset(timers_group, solvent%time_max_disp%name, solvent%time_max_disp%total)
  +  call h5md_write_dataset(timers_group, solvent%time_apply_pbc%name, solvent%time_apply_pbc%total)
     call h5md_write_dataset(timers_group, neigh%time_update%name, neigh%time_update%total)
     call h5md_write_dataset(timers_group, varia%name, varia%total)
     call h5md_write_dataset(timers_group, neigh%time_force%name, neigh%time_force%total)
  
positive.code: 
    end function compute_force_harmonic_trap
  
  end program setup_sphere_thermo_trap
  
negative.code:   colloids% force_old = colloids% force
  
    i = 0
    wall_v = 0
  
    write(*,*) 'Running for', N_loop, 'loops'
  