anchor.diff:        CALL TRIPAK(EIGVEC,H0,2,NPRVAR,NPRVAR)
   C          TRIPAK(AUTPAK,APAK,IWAY,MATDIM,NDIM)
   *. and diagonalize
  -      CALL MEMMAN(KLSCRVEC,NSEL,'ADDL  ',2,'SCRVEC')
  +      CALL MEMMAN(KLSCRVEC,NSEL,'ADDL  ',2,'SCRVEC') !done
   C          DIAG_SYMMAT_EISPACK(A,EIGVAL,SCRVEC,NDIM,IRETURN)
         CALL DIAG_SYMMAT_EISPACK(EIGVEC,EIGVAL,dbl_mb(KLSCRVEC),NPRVAR,
        &     IRETURN)
  
positive.code: *. Obtain subspace from lowest elements of CI diagonal
  *
  *. Local scratch
          CALL MEMMAN(KL1,3*(MXDM+1),'ADDL  ',1,'KL1   ')
          CALL MEMMAN(KL2,2*(MXDM+1),'ADDL  ',2,'KL2   ')
          CALL MEMMAN(KL3,2*(MXDM+1),'ADDL  ',2,'KL3   ')
          WRITE(6,*) ' MXDM, KL1, KL2, KL3 = ', KL1, KL2, KL3
          LBLK = -1
  *. And determine total subspace space 
  
negative.code:      &     int_mb(KCNOCCLS_BAT),int_mb(KCIBOCCLS_BAT),ILTEST)
  *. Space for strings
        IF(NOCSF.EQ.1) THEN
          CALL MEMMAN(KLASTR,MXNSTR*NAEL,'ADDL  ',1,'KLASTR') !done
          CALL MEMMAN(KLBSTR,MXNSTR*NBEL,'ADDL  ',1,'KLBSTR') !done
  *. 
  C       GET_MINMAX_ADR_IN_CISPACE_SD(
  C    &           IADR,NDET_UT,MINAC,MAXAC,MINMAX_ORB,NSSOA,NSSOB,NOCTPA,NOCTPB,
  