subroutine assignPSingles()
use constants, only: stdout
! set the initial values for pSingles/pDoubles
pSingles = 0.9_dp
write(stdout, '(A,F14.6)') " pDoubles set to: ", 1.0_dp - pSingles
write(stdout, '(A,F14.6)') " pSingles set to: ", pSingles
end subroutine assignPSingles