anchor.diff:     CALL AllocAry( CB%PhiR,   nL, nR,    'CB%PhiR',   ErrStat2, ErrMsg2 ); CALL SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
      CALL AllocAry( CB%OmegaL, nM_Out,    'CB%OmegaL', ErrStat2, ErrMsg2 ); if(Failed()) return
   
  -   CALL CraigBamptonReduction(Init%M, Init%K, IDR, nR, p%ID__L, nL, nM, nM_Out, CB%MBB, CB%MBM, CB%KBB, CB%PhiL, CB%PhiR, CB%OmegaL, ErrStat2, ErrMsg2,&
  +   CALL CraigBamptonReduction(Init%M, Init%K, p%ID__Rb, nR, p%ID__L, nL, nM, nM_Out, CB%MBB, CB%MBM, CB%KBB, CB%PhiL, CB%PhiR, CB%OmegaL, ErrStat2, ErrMsg2,&
                                 Init%FG, FGR, FGL, FGB, FGM)
      if(Failed()) return
   
  
positive.code:    type(SD_InputType),     intent(in   )  :: u ! Inputs
     type(SD_ParameterType), intent(in   )  :: p ! Parameters
     type(SD_MiscVarType),   intent(inout)  :: m ! Misc, for storage optimization of Fext and Fext_red
     real(ReKi)          ,   intent(out)    :: UFL(p%nDOFL_L)
     integer :: iMeshNode, iSDNode ! indices of u-mesh nodes and SD nodes
     integer :: nMembers
     real(ReKi), parameter :: myNaN = -9999998.989_ReKi 
  
negative.code:    ! write CB system matrices
     !-------------------------------------------------------------------------------------------------------------
     WRITE(UnSum, '(A)') SectionDivide
     WRITE(UnSum, '(A)') 'CB Matrices (PhiM,PhiR) (constraint applied)'
     
     WRITE(UnSum, '(A)') SubSectionDivide
     IF (p%nDOFM > 0) THEN
  