anchor.diff:            KLJASTR = KLIASTR
             KLJBSTR = KLIBSTR
           ELSE
  -          CALL MEMMAN(KLJASTR,NPRVAR*NAEL,'ADDL  ',1,'LJASTR')
  -          CALL MEMMAN(KLJBSTR,NPRVAR*NBEL,'ADDL  ',1,'LJBSTR')
  +          CALL MEMMAN(KLJASTR,NPRVAR*NAEL,'ADDL  ',1,'LJASTR') !done
  +          CALL MEMMAN(KLJBSTR,NPRVAR*NBEL,'ADDL  ',1,'LJBSTR') !done
             CALL ICOPVE(int_mb(KLIASTR),int_mb(KLJASTR),NPRVAR*NAEL)
             CALL ICOPVE(int_mb(KLIBSTR),int_mb(KLJBSTR),NPRVAR*NBEL)
           END IF
  
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: COLD &           NGAS,NORB,NACOB,NINOB)
        IF(NOCSF.EQ.1) THEN
  *. Obtain alpha and beta-strings for the selected determinants
          CALL MEMMAN(KLIASTR,NPRVAR*NAEL,'ADDL  ',1,'IASTR ') !done
          CALL MEMMAN(KLIBSTR,NPRVAR*NBEL,'ADDL  ',1,'IBSTR ') !done
  C       GET_IAIB_FOR_SEL_DETS(ISM,ISPC,ISEL,NSEL,IA,IB)
          CALL GET_IAIB_FOR_SEL_DETS(ISM,ISPC,ISEL,NPRVAR,
       &        int_mb(KLIASTR),int_mb(KLIBSTR))
  