anchor.diff:      
      ! --- Eigenvalue  analysis
      ! note: SGGEV seems to have memory issues in certain cases. The eigenvalues seem to be okay, but the eigenvectors vary wildly with different compiling options.
  -   !       DGGEV seems to work better, so I'm making these variables FEKi (which is set to R8Ki for now)   - bjj 4/25/2014
  +   !       DGGEV seems to work better, so I'm making these variables LaKi (which is set to R8Ki for now)   - bjj 4/25/2014
      ! bjj: This comes from the LAPACK documentation:
      !   Note: the quotients AlphaR(j)/BETA(j) and AlphaI(j)/BETA(j) may easily over- or underflow, and BETA(j) may even be zero.
      !   Thus, the user should avoid naively computing the ratio Alpha/beta.  However, AlphaR and AlphaI will be always less
  
positive.code:             call WrScr('[WARN] Complex eigenvalue found with large imaginary value)')
              Omega2(I) = MAX_EIGENVALUE
           endif
           call Fatal('Complex eigenvalue found, system may be ill-conditioned'); return
        endif
        ! Capping to avoid overflow
        if (Omega2(I)> MAX_EIGENVALUE) then
  
negative.code:    call BreakSysMtrx(MM, KK, IDR, IDL, nR, nL, MRR, MLL, MRL, KRR, KLL, KRL, FG=FG, FGR=FGR, FGL=FGL, CC=CC, CRR=CRR, CLL=CLL, CRL=CRL)
     ! --- CB reduction
     call CraigBamptonReduction_FromPartition( MRR, MLL, MRL, KRR, KLL, KRL, nR, nL, nM, nM_Out,& !< Inputs 
                                MBB, MBM, KBB, PhiL, PhiR, OmegaL, ErrStat2, ErrMsg2, & !< Outputs
                                CRR=CRR, CLL=CLL, CRL=CRL,& !< Optional inputs
                                CBB=CBB, CBM=CBM, CMM=CMM)  !< Optional Outputs
     if(Failed()) return
  
  
     ! --- Reduction of force if provided
     if (present(FG).and.present(FGR).and.present(FGL)) then
  