55
55
if (debug) print *, "HS read: nh, nsp, nnao: ", nh, nspin, nnao
56
56
if (nnao.ne.nao) STOP " nnao .ne. nao in HS"
57
58
if (wfs_x.and.(nspin.ne.nsp)) STOP " nspin .ne. nsp in HS"
59
61
allocate (numh(nao), listhptr(nao), listh(nh))
61
63
allocate (hamilt(nh,nspin))