anchor.diff:               Tmp_y2(1: 6)      = yout( L : L+5 )
                Tmp_y2(7:12)      = yout( L2 : L2+5 )
                CALL CALC_NODE_FORCES( DIRCOS,p%MOutLst(I)%Me(:,:,J,1),p%MOutLst(I)%Ke(:,:,J,1),Tmp_Udotdot, &
  -                              Tmp_y2,p%MoutLst(I)%Fg(:,J,1), K2,FM_elm,FK_elm) 
  +                              Tmp_y2,p%MOutLst(I)%Fg(:,J,1), K2,FM_elm,FK_elm) 
                
                FM_elm2=sgn*FM_elm
                FK_elm2=sgn*FK_elm
  
positive.code:                 L2=p%IDY((K3(2)-1)*6+1)! starting index for node K3(2) within yout
                  DIRCOS=transpose(p%elemprops(K)%DirCos)! global to local
                  CALL CALC_NODE_FORCES( DIRCOS,p%MOutLst2(I)%Me2(:,:,J),p%MOutLst2(I)%Ke2(:,:,J),(/uddout( L : L+5  ),uddout( L2 : L2+5 )/), &
                                   (/yout( L : L+5 ),       yout( L2 : L2+5 )/), p%MoutLst2(I)%Fg2(:,J),  K2,FM_elm,FK_elm) 
                   ! Store in All Outs
                   L=MaxOutPts+(I-1)*24+(J-1)*12+1!start index
                   L2=L+11
  
negative.code:                 DIRCOS=transpose(p%elemprops(K)%DirCos)! global to local
                  
                  CALL CALC_NODE_FORCES( DIRCOS,p%MOutLst3(I)%Me(:,:,1,J),p%MOutLst3(I)%Ke(:,:,1,J),(/uddout( L : L+5  ),uddout( L2 : L2+5 )/), &
                                   (/yout( L : L+5 ),       yout( L2 : L2+5 )/),  p%MOutLst3(I)%Fg(:,1,J),   K2,FM_elm,FK_elm) 
                  !transform back to global, need to do 3 at a time since cosine matrix is 3x3
                  DO L=1,2  
                     FM_elm2((L-1)*3+1:L*3) = FM_elm2((L-1)*3+1:L*3) + matmul(p%elemprops(K)%DirCos,FM_elm((L-1)*3+1:L*3))  !sum forces at joint in GLOBAL REF
  