This subroutine reads K_PATH and Q_PATH cards from standard input if present.
IMPORTANT: some cards may be missing, but the present ones must be ordered as described below.
K_PATH
label(1) k_x(1) k_y(1) k_z(1)
label(2) k_x(2) k_y(2) k_z(2)
...
label(nkspecial) k_x(nkspecial) k_y(nkspecial) k_z(nkspecial)
Q_PATH
nqpath nqspecial
label(1) q_x(1) q_y(1) q_z(1)
label(2) q_x(2) q_y(2) q_z(2)
...
label(nqspecial) q_x(nqspecial) q_y(nqspecial) q_z(nqspecial)
where nkpath indicate the desired total number of k-points along the path,
and nkspecial specifies the number of intermediate special points,
which are given below in crystal coordinates.
MBR 03/05/2024