subroutine molpro_ci_coeff
integer :: iCI, idxAlphaBetaOrbs(nbasis)
class(CI_coefficients_t), allocatable :: CI_coeff(:)
idxAlphaBetaOrbs = findAlphaBetaOrbs(symmax, nbasis)
do iCI = 1, n_store_ci_level
call read_ci_coeff(iCI, idxAlphaBetaOrbs, CI_coeff)
call dyn_sort_ci_coeff(CI_coeff)
call dyn_write_ci_coeff(CI_coeff)
end do
end subroutine molpro_ci_coeff