anchor.diff:                                                             'chem_pot: chemical potential not supplied.')
   
           opts = aot_table_top(lua_state)
  -        call read_kinetic_args(lua_state, opts, fermi_temperature, beta, nattempts, ncycles, seed, have_seed)
  +        call read_kinetic_args(lua_state, opts, fermi_temperature, beta, nattempts, ncycles, rng_seed, have_seed)
           call aot_table_close(lua_state, opts)
   
           if (have_seed) then
  -            call estimate_kinetic_energy(sys, fermi_temperature, beta, nattempts, ncycles, seed)
  +            call estimate_kinetic_energy(sys, fermi_temperature, beta, nattempts, ncycles, rng_seed)
           else
               call estimate_kinetic_energy(sys, fermi_temperature, beta, nattempts, ncycles)
           end if
  
positive.code:         use const, only: p
  
          use flu_binding, only: flu_State, flu_copyptr
          use aot_top_module, only: aot_top_get_val
          use aot_table_module, only: aot_table_top, aot_get_val, aot_exists, aot_table_close
          use aot_vector_module, only: aot_get_val
  
          use calc, only: ms_in
          use errors, only: stop_all
          use lua_hande_system, only: get_sys_t
          use lua_hande_utils, only: warn_unused_args
          use system, only: sys_t
  
          use canonical_kinetic_energy, only: estimate_kinetic_energy
  
negative.code: 
          type(c_ptr) :: sys_ptr
          type(sys_t), pointer :: sys
          integer :: opts, err, rng_seed, ncycles, nattempts
          logical :: fermi_temperature, have_seed
          real(p) :: beta
  
  