anchor.diff:       end if
   
        call varia%tic()
  -     call compute_rho_xy
  +     if (store_rho_xy) then
  +        call compute_rho_xy
  +        call rho_xy_el%append(rho_xy)
  +     end if
        call varia%tac()
  -     call rho_xy_el%append(rho_xy)
   
        temperature = compute_temperature(solvent, solvent_cells)
        if (dimer) then
  
positive.code:   prob = PTread_d(config,'probability')
  
    number_of_angles = PTread_i(config, 'number_of_angles')
    dimer = PTread_l(config, 'dimer')
    if (dimer) then
       N_colloids = 2
  
negative.code: 
    double precision :: g(3) !gravity
    logical :: fixed, on_track, stopped,order,dimer 
    logical :: store_rho_xy
    integer :: bufferlength
    double precision :: max_speed, z, Lz
    integer :: steps_fixed
  