anchor.diff:  
           use json_out
           use calc, only: doing_dmqmc_calc, dmqmc_energy, dmqmc_energy_squared, dmqmc_correlation, &
  -                        dmqmc_staggered_magnetisation, dmqmc_rdm_r2, dmqmc_full_r2
  +                        dmqmc_staggered_magnetisation, dmqmc_rdm_r2, dmqmc_full_r2, dmqmc_kinetic_energy, &
  +                        dmqmc_potential_energy, dmqmc_H0_energy, dmqmc_HI_energy
   
           type(json_out_t), intent(inout) :: js
           logical, intent(in), optional :: terminal
  
positive.code:         call json_object_init(js, 'operators')
          call json_write_key(js, 'energy', doing_dmqmc_calc(dmqmc_energy))
          call json_write_key(js, 'energy_squared', doing_dmqmc_calc(dmqmc_energy_squared))
          call json_write_key(js, 'correlation_fn', doing_dmqmc_calc(dmqmc_correlation))
          call json_write_key(js, 'staggered_mad_ind', doing_dmqmc_calc(dmqmc_staggered_magnetisation))
          call json_write_key(js, 'rdm_r2', doing_dmqmc_calc(dmqmc_rdm_r2))
  
negative.code:         call json_write_key(js, 'calc_excit_dist', dmqmc%calc_excit_dist)
          call json_write_key(js, 'all_sym_sectors', dmqmc%all_sym_sectors)
          call json_write_key(js, 'all_spin_sectors', dmqmc%all_spin_sectors)
          call json_write_key(js, 'initiator_level', dmqmc%initiator_level)
          if (allocated(dmqmc%sampling_probs)) then
              call json_write_key(js, 'sampling_probs', dmqmc%sampling_probs)
          else
  