rezero_double_occ_stats Subroutine

    subroutine rezero_double_occ_stats
        ! at the end of each cycle i should rezero the non-summed and
        ! non-averaged quantities to 0
        character(*), parameter :: this_routine = "rezero_double_occ_stats"

        ! there is more stuff to come..
        inst_double_occ = 0.0_dp

        if (t_spin_measurements) then
            inst_spatial_doub_occ = 0.0_dp
        end if

    end subroutine rezero_double_occ_stats