anchor.diff:    type(lj_params_t) :: colloid_lj
     type(lj_params_t) :: walls_colloid_lj
   
  -  type(profile_t) :: tz
  -  type(histogram_t) :: rhoz
     type(profile_t) :: vx
   
     integer :: rho
  
positive.code:   integer(HID_T) :: box_group
  
    type(PTo) :: config
    integer(c_int64_t) :: seed
    integer :: i, L(3),  n_threads
    integer :: j, k, m
    
  
negative.code: 
    T = PTread_d(config, 'T')
    d = PTread_d(config, 'd')
    order = PTread_l(config, 'order')
  
    wall_v = 0
    wall_t = [T, T]
  