<div dir="ltr">I recently ran the code with array boundary checking and it crashed in vloc_psi.f90. In particular the array igk_k has components that are 0 and then accessing nls of these elements is not within the bounds of nls (which is 1:ngms).<br><br>Martin<br><br><font size="1">Martin Schlipf <br>Department of Materials<br>University of Oxford<br>Parks Road<br>Oxford OX1 3PH, UK</font></div>