anchor.diff:      N(:) = N(:) + NFIRST
       NLAST  = N(9)
       
  -    call HexTwnSeven (IA(NP(1)),DA(NP(2)),DA(NP(3)),DA(NP(4)),DA(NP(4)),IA(NP(5)),   &
  +    call HexTwn (IA(NP(1)),DA(NP(2)),DA(NP(3)),DA(NP(4)),DA(NP(4)),IA(NP(5)),   &
                     A(N(1)),A(N(2)),A(N(3)),A(N(4)),A(N(5)),A(N(6)),A(N(7)),A(N(8)))
       
       !Reuse DA(NP(4)) at Solution Phase 3 as displacement U
       return
   end subroutine HexTriQuad
   
  -subroutine HexTwnSeven (ID,X,Y,Z,U,MHT,E, PoissonRatio, Density, Gravity, LM, PositionData, MaterialData, Node)
  +subroutine HexTwn (ID,X,Y,Z,U,MHT,E, PoissonRatio, Density, Gravity, LM, PositionData, MaterialData, Node)
       use globals
       use MemAllocate
       use MathKernel
  
positive.code:     logical ::  OK_Flag
      
      select case (ElementShapeNodes)
      case (27)
      !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>WORKING PROGRESS
          xi      = Transformed(1)
          eta     = Transformed(2)
  
negative.code:     integer :: NumberOfElements, NumberOfMaterials, ElementGroupSize, QuadratureOrder=2
      integer :: N(11) !Pointers
      
      NPAR(5) = 20
      NPAR(4) = 0
      
      NumberOfElements = NPAR(2)
  