Commit e6c0f501 authored by Daniel Wortmann's avatar Daniel Wortmann

Added STOP if q-vector for self-consistency is not the first in force-theorem...

Added STOP if q-vector for self-consistency is not the first in force-theorem loop. Fixes part of #246
parent 3615dabb
......@@ -71,6 +71,7 @@ CONTAINS
q_vectors(:,1)=noco%qss
END SELECT
if (any(abs(noco%qss-q_vectors(:,1))>1E-4)) CALL judft_warn("q-vector for self-consistency should be first in list for force-theorem")
noco%qss=q_vectors(:,1) ! Usually does not do anything, but ensures that in
!force theorem CASE noco%qss is first q-vector in list
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment