anchor.diff:         ELSE
            LEN_MAT = LEN_FP
          END IF
  -       CALL VECSUM(FI,H,WORK(KLBLF),ONE,ONE,LEN_MAT)
  +       CALL VECSUM(FI,H,dbl_mb(KLBLF),ONE,ONE,LEN_MAT)
   *
          IF(NTEST.GE.100) THEN
            WRITE(6,*) ' Output from FI_FROM_INIINT_G: '
  
positive.code:       DO IFRZ_T = 1, NFRZ_ORB
          IFRZ_S = IREOTS(IFRZ_ORB(IFRZ_T))
  C            ICOPVE3(IIN,IOFFIN,IOUT,IOFFOUT,NDIM)
          CALL ICOPVE3(IFRZ_S,1,WORK(KLFRZS),IFRZ_T,1)
        END DO
  *
        IF(NTEST.GE.10) THEN
          WRITE(6,*) ' Frozen orbitals in symmetry numbering '
          CALL IWRTMA(WORK(KLFRZS),1,NFRZ_ORB,1,NFRZ_ORB)
        END IF
  *. And orthonormalize
  C     ORT_MOS_TO_SELECTED_MOS(CMOAO,NORT,IORT)
        CALL ORT_MOS_TO_SELECTED_MOS(CMOAO,NFRZ_ORB,WORK(KLFRZS))
  *
        IF(NTEST.GE.100) THEN 
          WRITE(6,*) ' Orbitals orthogonalized to frozen '
  
negative.code:       ELSE
          VFINAL = ERROR(MAXIT+1)
        END IF
        IF(NTEST .GT. 0 ) THEN
  *
        IF(CONVER) THEN
  