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