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:        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_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
              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_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
              CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
             
         ENDDO
  
negative.code:                    CALL BD_ComputeIniNodalPosition(temp_Coef,eta,temp_POS,temp_e1,temp_twist,&
                                                     ErrStat2, ErrMsg2)
                        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
                     CALL BD_ComputeIniNodalCrv(temp_e1,temp_twist,temp_CRV,ErrStat2,ErrMsg2)
                        CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
                     temp_id2 = (i-1)*p%ngp+j+1
                     temp_L2(1:3,temp_id2) = temp_POS(1:3)
                     temp_L2(4:6,temp_id2) = temp_CRV(1:3)
                     EXIT
                 ENDIF
             ENDDO
         ENDDO
     ENDDO
     temp_L2(1:3,1) = p%uuN0(1:3,1)
     temp_L2(4:6,1) = p%uuN0(4:6,1)
     temp_L2(1:3,p%ngp*p%elem_total+2) = p%uuN0(temp_int-5:temp_int-3,p%elem_total)
     temp_L2(4:6,p%ngp*p%elem_total+2) = p%uuN0(temp_int-2:temp_int,p%elem_total)
     DEALLOCATE(temp_GLL)
     DEALLOCATE(SP_Coef)
  
  