do i=1,#energy
emult(i)=energy(i)
enddo
which redefines the outer loop variable i, which is referenced in the geometry. This causes the hamiltonian evaluator to get confused. You can simply change i to j and it will work - or, more simply,
emult = energy
although emult doesn't get used, so maybe you intend something else.
Peter
External email to Cardiff University - Take care when replying/opening attachments or links. |
Nid ebost mewnol o Brifysgol Caerdydd yw hwn -
Cymerwch ofal wrth ateb/agor atodiadau neu ddolenni. |