anchor.diff:             temp_id = (j-1) * p%dof_node
              TmpPos(1:3) = p%GlbPos(1:3) + MATMUL(p%GlbRot,p%uuN0(temp_id+1:temp_id+3,i))
              temp_CRV(:) = MATMUL(p%GlbRot,p%uuN0(temp_id+4:temp_id+6,i))
  +           CALL BD_CrvCompose(temp_CRV,temp_CRV,temp_glb,0,ErrStat2,ErrMsg2)
  +              CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
              CALL BD_CrvMatrixR(temp_CRV,TmpDCM,ErrStat2,ErrMsg2)
  +              CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
              TmpDCM(:,:) = TRANSPOSE(TmpDCM)
              temp_id = (i-1)*(p%node_elem-1)+j
              CALL MeshPositionNode ( Mesh    = u%PointLoad  &
  
positive.code:    DO i=1,p%ngp*p%elem_total+2
         TmpPos(1:3) = p%GlbPos(1:3) + MATMUL(p%GlbRot,temp_L2(1:3,i))
         temp_CRV(:) = MATMUL(p%GlbRot,temp_L2(4:6,i))
         CALL BD_CrvMatrixR(temp_CRV,TmpDCM,ErrStat2,ErrMsg2)
         TmpDCM(:,:) = TRANSPOSE(TmpDCM)
         CALL MeshPositionNode ( Mesh    = u%DistrLoad  &
                                ,INode   = i            &
  
negative.code:            temp_id = (j-1)*p%dof_node
             TmpPos(1:3) = p%GlbPos(1:3) + MATMUL(p%GlbRot,p%uuN0(temp_id+1:temp_id+3,i))
             temp_CRV(:) = MATMUL(p%GlbRot,p%uuN0(temp_id+4:temp_id+6,i))
             CALL BD_CrvCompose(temp_CRV,temp_CRV,temp_glb,0,ErrStat2,ErrMsg2)
                CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
             CALL BD_CrvMatrixR(temp_CRV,TmpDCM,ErrStat2,ErrMsg2)
                CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
             TmpDCM(:,:) = TRANSPOSE(TmpDCM)
             temp_id = (i-1)*p%node_elem+j
             CALL MeshPositionNode ( Mesh    = y%BldMotion  &
  