anchor.diff:         DO j=1,p%node_elem
              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_CrvMatrixR(temp_CRV,TmpDCM)
  +           TmpDCM(:,:) = TRANSPOSE(TmpDCM)
              temp_id = (i-1)*(p%node_elem-1)+j
              CALL MeshPositionNode ( Mesh    = u%PointLoad  &
                                     ,INode   = temp_id      &
                                     ,Pos     = TmpPos       &
                                     ,ErrStat = ErrStat2      &
                                     ,ErrMess = ErrMsg2       &
  -                                  , Orient = ??? ) !bjj: add orient=DCM ! Orientation (direction cosine matrix) of node; identity by default
  +                                  , Orient = TmpDCM ) !bjj: add orient=DCM ! Orientation (direction cosine matrix) of node; identity by default
               CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
          ENDDO
      ENDDO
      DO i=1,p%elem_total
          DO j=1,p%node_elem
              temp_id = (j-1)*p%dof_node
  -!           TmpPos(1:3) = p%uuN0(temp_id+1:temp_id+3,i)
              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_CrvMatrixR(temp_CRV,TmpDCM)
  +           TmpDCM(:,:) = TRANSPOSE(TmpDCM)
              temp_id = (i-1)*p%node_elem+j
              CALL MeshPositionNode ( Mesh    = y%BldMotion  &
                                     ,INode   = temp_id      &
                                     ,Pos     = TmpPos       &
                                     ,ErrStat = ErrStat2      &
                                     ,ErrMess = ErrMsg2       &
  -                                  , Orient = ??? ) !bjj: add orient=DCM ! Orientation (direction cosine matrix) of node; identity by default
  +                                  , Orient = TmpDCM ) !bjj: add orient=DCM ! Orientation (direction cosine matrix) of node; identity by default
               CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
              
          ENDDO
  
positive.code:                    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)
                     EXIT
                 ENDIF
             ENDDO
  
negative.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)
         TmpDCM(:,:) = TRANSPOSE(TmpDCM)
         CALL MeshPositionNode ( Mesh    = u%DistrLoad  &
                                ,INode   = i            &
                                ,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 )
         
     ENDDO
  