anchor.diff:        p%MMB = TRANSPOSE( p%MBM )                          != MMBt
         p%CMB = TRANSPOSE( p%CBM )                          != CMBt
   
  -      p%PhiM  = PhiL(:,1:p%nDOFM)
  +      p%PhiM = real( PhiL(:,1:p%nDOFM), ReKi)
         
         ! A_21=-Kmm (diagonal), A_22=-Cmm (approximated as diagonal) 
         p%KMMDiag=             OmegaL(1:p%nDOFM) * OmegaL(1:p%nDOFM)          ! OmegaM is a one-dimensional array
         p%CMMDiag = 2.0_ReKi * OmegaL(1:p%nDOFM) * Init%JDampings(1:p%nDOFM)  ! Init%JDampings is also a one-dimensional array
   
  -   
  -      ! B_23, B_24
  -      !p%PhiM_T =  TRANSPOSE( p%PhiM  )
  -   
         ! FX = matmul( transpose(PhiM), FGL ) (output of CraigBamptonReduction)
         p%FX = FGM
      
  
positive.code:       ENDDO   
        
        ! D2_53, D2_63, D2_64 
        p%D2_63 = MATMUL( p%PhiM, p%MMB )
        p%D2_63 = p%PhiRb_TI - p%D2_63
  
        !p%D2_64 = MATMUL( p%PhiM, p%PhiM_T )
  
negative.code:       else
           CALL WrScr('   Using static improvement method for gravity only')
        endif
        ! Allocations - NOTE: type conversion belows from FEKi to ReKi
        CALL AllocAry( p%PhiL_T,        p%nDOF__L, p%nDOF__L, 'p%PhiL_T',        ErrStat2, ErrMsg2 ); if(Failed())return
        CALL AllocAry( p%PhiLInvOmgL2,  p%nDOF__L, p%nDOF__L, 'p%PhiLInvOmgL2',  ErrStat2, ErrMsg2 ); if(Failed())return
        CALL AllocAry( p%KLLm1       ,  p%nDOF__L, p%nDOF__L, 'p%KLLm1',         ErrStat2, ErrMsg2 ); if(Failed())return
  