anchor.diff:    type(thermo_t) :: thermo_data
     double precision :: temperature, kin_e
     double precision :: v_com(3), x(3)
  +  double precision :: wall_v(3,2)
     double precision :: tmp_mass
     type(particle_system_io_t) :: janus_io
     type(particle_system_io_t) :: solvent_io
  
positive.code:         co_max = colloids% maximum_displacement()
  
          if ( (co_max >= skin*0.1d0) .or. (so_max >= skin*0.9d0) ) then
             if ((do_ywall) .and. (solvent_cells%bc(2)/=PERIODIC_BC)) then
                call cell_md_pos_ywall(solvent_cells, solvent, (N_MD_steps*i+j - loop_i_last_sort)*dt, md_flag=.false.)
             else
                call cell_md_pos(solvent_cells, solvent, (N_MD_steps*i+j - loop_i_last_sort)*dt, md_flag=.false.)
             end if
  
negative.code: 
    call polar%init(N_species, 64, sigma, polar_r_max, 64)
    call timer_list%append(polar%time_polar_update)
    call planar%init(N_species, 64, -12.d0, 12.d0, 64, -12.d0, 12.d0, 1.d0)
    call timer_list%append(planar%timer)
    call neigh% init(colloids% Nmax, int(500*max(sigma,1.15d0)**3))
  
  