anchor.diff:      CALL Eye( outputMesh%Orientation, ErrStat, ErrMsg )         
           
   END SUBROUTINE CreateY2Meshes
  -!------------------------------------------------------------------------------------------------------
  -!> Set the index array that maps SD internal nodes to the Y2Mesh nodes.
  -!! NOTE: SDtoMesh is not checked for size, nor are the index array values checked for validity, 
  -!!       so this routine could easily have segmentation faults if any errors exist.
  -SUBROUTINE SD_Y2Mesh_Mapping(p, SDtoMesh )
  -   TYPE(SD_ParameterType),       INTENT(IN   )  :: p           !< Parameters
  -   INTEGER(IntKi),               INTENT(  OUT)  :: SDtoMesh(:) !< index/mapping of mesh nodes with SD mesh
  -   ! locals
  -   INTEGER(IntKi)                               :: i
  -   INTEGER(IntKi)                               :: SDnode
  -   INTEGER(IntKi)                               :: y2Node
  -
  -   y2Node = 0
  -   ! Interface nodes (IDI)
  -   DO I = 1,SIZE(p%IDI,1)/6
  -      y2Node = y2Node + 1      
  -      SDnode = p%IDI(I*6) / 6     ! TODO TODO TODO integer division gives me the actual node index; it is not the nodeID
  -      SDtoMesh( SDnode ) = y2Node ! TODO add safety check
  -   END DO
  -   
  -   ! Interior nodes (IDL)
  -   DO I = 1,SIZE(p%IDL,1)/6 
  -      y2Node = y2Node + 1      
  -      SDnode = p%IDL(I*6) / 6     ! TODO TODO TODO integer division gives me the actual node index; it is not the nodeID
  -      SDtoMesh( SDnode ) = y2Node ! TODO add safety check
  -   END DO
  -
  -   ! Base Reaction nodes (IDC)
  -   DO I = 1,SIZE(p%IDC,1)/6 
  -      y2Node = y2Node + 1      
  -      SDnode = p%IDC(I*6) / 6     ! TODO TODO TODO integer division gives me the actual node index; it is not the nodeID
  -      SDtoMesh( SDnode ) = y2Node ! TODO add safety check
  -   END DO
  -
  -END SUBROUTINE SD_Y2Mesh_Mapping
   
   
   !---------------------------------------------------------------------------
  
positive.code:                 XYZ2=Init%Nodes(Init%Members(i,3),2:4)
             ENDIF
         ENDDO    
         CALL GetDirCos(XYZ1(1:3), XYZ2(1:3), DirCos, mlength, ErrStat, ErrMsg)
         DirCos=TRANSPOSE(DirCos) !This is now global to local
         WRITE(UnSum, '(I9,9(E15.6))') Init%Members(i,1), ((DirCos(k,j),j=1,3),k=1,3)
     ENDDO
  
negative.code:    CHARACTER(1024)        :: SummaryName    ! name of the SubDyn summary file
     INTEGER(IntKi)         :: i, j, k, propids(2)  !counter and temporary holders
     INTEGER(IntKi)         :: mType ! Member Type
     Real(ReKi)             :: mMass, mLength ! Member mass and length
     INTEGER(IntKi)         :: SDtoMeshIndx(Init%NNode)
     REAL(ReKi)             :: MRB(6,6)    !REDUCED SYSTEM Kmatrix, equivalent mass matrix
     REAL(ReKi)             :: XYZ1(3),XYZ2(3), DirCos(3,3) !temporary arrays, member i-th direction cosine matrix (global to local) and member length
     CHARACTER(*),PARAMETER                 :: SectionDivide = '____________________________________________________________________________________________________'
     CHARACTER(*),PARAMETER                 :: SubSectionDivide = '__________'
     CHARACTER(2),  DIMENSION(6), PARAMETER :: MatHds= (/'X ', 'Y ', 'Z ', 'XX', 'YY', 'ZZ'/)  !Headers for the columns and rows of 6x6 matrices
  