anchor.diff:  
      WRITE(UnSum, '()') 
      WRITE(UnSum, '(A,I6)')  'Number of elements (NElems):',Init%NElem
  -   WRITE(UnSum, '(A10,5(A10))')  'Elem No.',    'Node_I',     'Node_J',   'Type',    'Prop_I',      'Prop_J'
  +   WRITE(UnSum, '(A10,5(A10))') 'Elem No.','Node_I','Node_J','Prop_I','Prop_J','Type'
      WRITE(UnSum, '(6(I10))') ((p%Elems(i, j), j = 1, MembersCol), i = 1, Init%NElem)
      
      WRITE(UnSum, '()') 
  
positive.code: #endif   
     
     CALL SDOut_CloseSum( UnSum, ErrStat, ErrMsg )  
  
  END SUBROUTINE OutSummary
  
  
negative.code:    WRITE(UnSum, '()') 
     WRITE(UnSum, '(A,I6)')  'Number of members',p%NMembers
     WRITE(UnSum, '(A,I6)')  'Number of nodes per member:', Init%Ndiv+1
     WRITE(UnSum, '(A9,A10,A10,A10,A10,A15,A15,A16)')  'Member ID', 'Joint1_ID', 'Joint2_ID','Prop_I','Prop_J', 'Mass','Length', 'Node IDs...'
     DO i=1,p%NMembers
         !Calculate member mass here; this should really be done somewhere else, yet it is not used anywhere else
         !IT WILL HAVE TO BE MODIFIED FOR OTHER THAN CIRCULAR PIPE ELEMENTS
         propids=p%Elems(i,iMProp:iMProp+1) ! TODO use member
         mLength=MemberLength(Init%Members(i,1),Init,ErrStat,ErrMsg) ! TODO double check mass and length
         IF (ErrStat .EQ. ErrID_None) THEN
          mType =  Init%Members(I, iMType) ! 
          if (mType==idMemberBeam) then
             mMass= MemberMass(Init%PropSetsB(propids(1),4),Init%PropSetsB(propids(1),5),Init%PropSetsB(propids(1),6),   &
                               Init%PropSetsB(propids(2),4),Init%PropSetsB(propids(2),5),Init%PropSetsB(propids(2),6), mLength, .TRUE.)
  
             WRITE(UnSum, '(I9,I10,I10,I10,I10,E15.6,E15.6, A3,'//Num2LStr(Init%NDiv + 1 )//'(I6))') Init%Members(i,1:3),propids(1),propids(2),&
                   mMass,mLength,' ',(Init%MemberNodes(i, j), j = 1, Init%NDiv+1)
          else
             WRITE(UnSum, '(A)') 'TODO, member mass for non-beam elements'
          endif
  