anchor.diff:     DirCos = 0
      
      IF ( EqualRealNos(xz, 0.0) ) THEN
  -      IF( y2 < y1) THEN
  +      IF( y2 < y1) THEN  !x is kept along global x
            DirCos(1, 1) = 1.0
            DirCos(2, 3) = -1.0
            DirCos(3, 2) = 1.0
  
positive.code:    
     ! Initialize TempMembers and Elems
     DO I = 1, p%NMembers
        
        !IF ( SwapElementNodes( Init%Nodes(Init%Members(I, 2),:), Init%Nodes(Init%Members(I, 3),:) ) ) THEN
        !   TempMembers(I, 1) = I
        !   TempMembers(I, 4) = Init%Members(I, 5)
        !   TempMembers(I, 5) = Init%Members(I, 4)
        !
        !   p%Elems(I,     1) = I
        !   p%Elems(I, NPE+2) = Init%Members(I, 5)
        !   p%Elems(I, NPE+3) = Init%Members(I, 4)
        !
        !   Node1 = Init%Members(I, 3)
        !   Node2 = Init%Members(I, 2)
        !ELSE
           TempMembers(I, 1) = I
           TempMembers(I, 4) = Init%Members(I, 4)
           TempMembers(I, 5) = Init%Members(I, 5)
        
           p%Elems(I,     1) = I
           p%Elems(I, NPE+2) = Init%Members(I, 4)
           p%Elems(I, NPE+3) = Init%Members(I, 5)
        
           Node1 = Init%Members(I, 2)
           Node2 = Init%Members(I, 3)
        
        !END IF
        
        flg1 = 0
        flg2 = 0
  
negative.code:    INTEGER,                INTENT(IN   )   :: p1
     INTEGER,                INTENT(IN   )   :: p2
     TYPE(SD_ParameterType), INTENT(INOUT)   :: p
    
     
     ! Local Variables
     
     p%Elems(k, 1) = k
     p%Elems(k, 2) = n1
     p%Elems(k, 3) = n2
     p%Elems(k, 4) = p1
     p%Elems(k, 5) = p2
  
  END SUBROUTINE GetNewElem
  !------------------------------------------------------------------------------------------------------
  