anchor.diff:                     CALL BD_ComputeIniNodalCrv(temp_e1,temp_twist,temp_CRV,ErrStat2,ErrMsg2)
                         CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
                      temp_id2 = (j-1)*p%dof_node
  -                   p%uuN0(temp_id2+1,i) = temp_POS(1)
  -                   p%uuN0(temp_id2+2,i) = temp_POS(2)
  -                   p%uuN0(temp_id2+3,i) = temp_POS(3)
  -                   p%uuN0(temp_id2+4,i) = temp_CRV(1)
  -                   p%uuN0(temp_id2+5,i) = temp_CRV(2)
  -                   p%uuN0(temp_id2+6,i) = temp_CRV(3)
  +                   p%uuN0(temp_id2+1:temp_id2+3,i) = temp_POS(1:3)
  +                   p%uuN0(temp_id2+4:temp_id2+6,i) = temp_CRV(1:3)
                      EXIT
                  ENDIF
              ENDDO
  
positive.code:    ENDDO
     ! place single node at origin; position affects mapping/coupling with other modules
     TmpPos(:) = p%GlbPos(1:3) + MATMUL(p%GlbRot,p%uuN0(1:3,1))
     CALL MeshPositionNode ( Mesh = u%RootMotion      &
                           , INode = 1                &
                           , Pos = TmpPos             &
                           , ErrStat   = ErrStat2     &
                           , ErrMess   = ErrMsg2      &
                           , Orient = ??? ) !bjj: add orient=DCM ! Orientation (direction cosine matrix) of node; identity by default
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  
  
  
negative.code:    REAL(ReKi),ALLOCATABLE  :: temp_L2(:,:)
     REAL(ReKi),ALLOCATABLE  :: SP_Coef(:,:,:)
     REAL(ReKi)              :: TmpPos(3)
     REAL(ReKi)              :: TmpDCM(3,3)
  
     INTEGER(IntKi)          :: ErrStat2                     ! Temporary Error status
     CHARACTER(ErrMsgLen)    :: ErrMsg2                      ! Temporary Error message
  