anchor.diff:     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
  
positive.code:    DEALLOCATE(temp_w)
     ! temp_L2: the DistrLoad mesh node location
     ! temp_L2: physical coordinates of Gauss points and two end points
     CALL AllocAry(temp_L2,3,p%ngp*p%elem_total+2,'temp_L2',ErrStat2,ErrMsg2)
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
     ! Temporary Gauss point intrinsic coordinates array
     CALL AllocAry(temp_GL,p%ngp,'temp_GL',ErrStat2,ErrMsg2)
  
negative.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))
     temp_CRV(:) = MATMUL(p%GlbRot,p%uuN0(4:6,1))
     CALL BD_CrvMatrixR(temp_CRV,TmpDCM)
     TmpDCM(:,:) = TRANSPOSE(TmpDCM)
     CALL MeshPositionNode ( Mesh = u%RootMotion      &
                           , INode = 1                &
                           , Pos = TmpPos             &
                           , ErrStat   = ErrStat2     &
                           , ErrMess   = ErrMsg2      &
                           , Orient = TmpDCM ) !bjj: add orient=DCM ! Orientation (direction cosine matrix) of node; identity by default
        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
  
  
  