subroutine update_first_reference() use FciMCData, only: tSinglePartPhase implicit none call setup_reference_space(all(.not. tSinglePartPhase)) end subroutine update_first_reference