anchor.diff:                   L  = p%IDY((K3(1)-1)*6+1)! starting index for node K3(1) within yout ! TODO different DOF order
                    L2 = p%IDY((K3(2)-1)*6+1)! starting index for node K3(2) within yout
                    CALL CALC_NODE_FORCES(DIRCOS,p%MOutLst(I)%Me(:,:,J,p%NavgEls),p%MOutLst(I)%Ke(:,:,J,p%NavgEls),(/uddout( L : L+5  ),uddout( L2 : L2+5 )/), &
  -                                 (/yout( L : L+5 ),       yout( L2 : L2+5 )/), p%MoutLst(I)%Fg(:,J,p%NavgEls), K2,FM_elm,FK_elm ) 
  +                                 (/yout( L : L+5 ),       yout( L2 : L2+5 )/), p%MOutLst(I)%Fg(:,J,p%NavgEls), K2,FM_elm,FK_elm ) 
                                      
                    FM_elm2=0.5*( FM_elm2 + sgn*FM_elm ) !Now Average
                    FK_elm2=0.5*( FK_elm2 + sgn*FK_elm) !Now Average
  
positive.code:    DO I=1,9
            !I know el # and whether it is 1st node or second node
        if (I <= p%NMOutputs) then
           INDX=p%MoutLst(I)%NOutCnt+1
        else
           INDX = 1
        end if
  
negative.code:           !I need to get at most 2 elements that belong to the same I Member
            !make use of MemberNodes and NodesConnE
            
            NconEls=Init%NodesConnE(p%MOutLst2(I)%NodeIDs(J),1) !Number of elements connecting to the 1st or last node of the member
            
            K2= J/(Init%Ndiv+1)+1  !store this quantity used later, basically 1 or 2 depending on J
           
            DO K=1, NconEls 
                L=Init%NodesConnE(p%MOutLst2(I)%NodeIDs(J),k+1)  !k-th Element Number in the set of elements attached to the selected node 
                M=p%Elems(L,2:3) !1st and 2nd node of the k-th element
               !Select only the other node, not the one where elements connect to
               Junk=M(1)
               IF (M(1) .EQ. p%MOutLst2(I)%NodeIDs(J)) Junk=M(2)
               
               IF (ANY(Init%MemberNodes(p%MOutLst2(I)%MemberID,:) .EQ. Junk)) THEN  !This means we are in the selected member
                    p%MOutLst2(I)%ElmID2s(K2)=L     !This array has for each node requested NODEID(J), for each member I, the element to get results for 
                    p%MOutLst2(I)%ElmNd2s(K2)=1                        !store whether first or second node of element  
                    IF (M(2) .EQ. p%MOutLst2(I)%NodeIDs(J) ) p%MOutLst2(I)%ElmNd2s(K2)=2 !store whether first or second node of element  
                    ! --- Element Me, Ke, Fg, Fce
                    CALL ElemM(p%ElemProps(L),         p%MOutLst2(I)%Me2(:,:,K2))
                    CALL ElemK(p%ElemProps(L),         p%MOutLst2(I)%Ke2(:,:,K2))
                    CALL ElemF(p%ElemProps(L), Init%g, p%MOutLst2(I)%Fg2(:,K2), FCe)
                    p%MOutLst2(I)%Fg2(:,K2) = p%MOutLst2(I)%Fg2(:,K2) + FCe(1:12) ! Adding cable element force 
                     
                    EXIT   !We found the element for that node, exit loop on elements
                ENDIF
  