This subroutine scans along a prescribed plane in wavevector space to map out in this space. It is required that n_scan=2, and is invoked with option =2 KGK: This line causes the solution to (occasionally) smoothly transition to unphysical values. Suppressing until we understand the error.
KGK: This line causes the solution to (occasionally) smoothly transition to unphysical values. Suppressing until we understand the error.