Incomplete: part way implementation of boltzmann weight probability calculations