anchor.diff:  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 ')
  -        CALL MEMMAN(KLIBSTR,NPRVAR*NBEL,'ADDL  ',1,'IBSTR ')
  +        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))
  
positive.code:       CALL TRIPAK(EIGVEC,H0,2,NPRVAR,NPRVAR)
  C          TRIPAK(AUTPAK,APAK,IWAY,MATDIM,NDIM)
  *. and diagonalize
        CALL MEMMAN(KLSCRVEC,NSEL,'ADDL  ',2,'SCRVEC')
  C          DIAG_SYMMAT_EISPACK(A,EIGVAL,SCRVEC,NDIM,IRETURN)
        CALL DIAG_SYMMAT_EISPACK(EIGVEC,EIGVAL,dbl_mb(KLSCRVEC),NPRVAR,
       &     IRETURN)
  
negative.code:           KLJASTR = KLIASTR
            KLJBSTR = KLIBSTR
          ELSE
            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
  