Bugfixes: eigenso -> broke SOC, strgn -> variable not initialized

eigen_soc/eigenso.f90

@@ -264,7 +264,7 @@ CONTAINS
                eig=eig_so(:nsz))
 
        ELSE
-          zmat%nbasfcn=size(eig_so,1)
+          zmat%nbasfcn=size(zso,1)
           allocate(zmat%z_c(zmat%nbasfcn,nsz))
           zmat%l_real=.false.
           zmat%nbands=nsz        

init/strgn.f90

@@ -317,10 +317,11 @@ CONTAINS
     !
     !--->    store number of star with respect to z-index in igz
     !
-    ngz = izmax - izmin + 1
-    DO k = 1,stars%ng3
-       igz(k) = stars%kv3(3,k) - izmin + 1
-    ENDDO
+    !ngz = izmax - izmin + 1
+    !DO k = 1,stars%ng3
+    !   igz(k) = stars%kv3(3,k) - izmin + 1
+    !ENDDO
+    igz=0;ngz=0;izmin=0;izmax=0
 
     !--->    generate all star members
     !+gu
@@ -767,11 +768,11 @@ CONTAINS
     !
     !--->    store number of star with respect to z-index in igz
     !
-    ngz = izmax - izmin + 1
-    DO k = 1,stars%ng3
-       igz(k) = stars%kv3(3,k) - izmin + 1
-    ENDDO
-
+    !ngz = izmax - izmin + 1
+    !DO k = 1,stars%ng3
+    !   igz(k) = stars%kv3(3,k) - izmin + 1
+    !ENDDO
+    igz=0;ngz=0;izmax=0;izmin=0
     !--->    generate all star members
     !+gu
     kidx=0

io/eig66_mem.F90

@@ -66,10 +66,8 @@ CONTAINS
     ALLOCATE(d%eig_eig(neig,jspins*nkpts))
     !d%eig_vec
     if (l_real.and..not.l_soc) THEN
-       print *, "Allocate real in eig66_mem"
        ALLOCATE(d%eig_vecr(nmat*neig,length*nkpts))
     else
-       print *, "Allocate complex in eig66_mem"
        ALLOCATE(d%eig_vecc(nmat*neig,length*nkpts))
     endif
     length=length*nkpts
@@ -270,9 +268,6 @@ CONTAINS
     IF (PRESENT(eig)) THEN
        eig=0.0
        eig=d%eig_eig(:SIZE(eig),nrec)
-       !print *,"R-eig:",nrec,shape(eig)
-       !print*,"R-eig(data):",shape(d%eig_eig)
-       !print*,"R:",eig
     ENDIF
     !data from d%eig_vec
 
@@ -338,8 +333,6 @@ CONTAINS
     !data from d%eig_eig
     IF (PRESENT(eig)) THEN
        d%eig_eig(:SIZE(eig),nrec)=eig
-       !print*,"W-eig:",nrec,shape(eig)
-       !print*,"W:",eig
     ENDIF
     !data from d%eig_vec
     IF (PRESENT(zmat)) THEN