SUBROUTINE GetLz(nI, NElec, Lz) INTEGER :: NElec INTEGER :: nI(NElec), Lz, i Lz = 0 do i = 1, NElec Lz = Lz + G1(nI(i))%Ml end do END SUBROUTINE GetLz