Skip to content
Snippets Groups Projects
Commit 9a1545a1 authored by Philipp Rüssmann's avatar Philipp Rüssmann
Browse files

Fix orbital momentum writeout

This fixes issue #130 and it was found by Juba Bouaziz
parent 112784da
No related branches found
No related tags found
No related merge requests found
......@@ -760,7 +760,7 @@ contains
!----------------------------------------------------------------------------
! Orbital magnetic moments
!----------------------------------------------------------------------------
if (krel==1) then
if ((krel+korbit)==1) then
do ispin = 1, 3
do l = 0, lmax + 1
muorb(lmaxd1+1, ispin, i1) = muorb(lmaxd1+1, ispin, i1) + muorb(l, ispin, i1)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment