anchor.diff:      logical ::  OK_Flag
       
       select case (ElementShapeNodes)
  -    case (27)
  +    case (20)
       !>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>WORKING PROGRESS
           xi      = Transformed(1)
           eta     = Transformed(2)
  
positive.code:     integer :: NumberOfElements, NumberOfMaterials, ElementGroupSize, QuadratureOrder=2
      integer :: N(11) !Pointers
      
      NPAR(5) = 27
      NPAR(4) = 0
      
      NumberOfElements = NPAR(2)
  
negative.code:     N(:) = N(:) + NFIRST
      NLAST  = N(9)
      
      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 HexTwn (ID,X,Y,Z,U,MHT,E, PoissonRatio, Density, Gravity, LM, PositionData, MaterialData, Node)
      use globals
      use MemAllocate
      use MathKernel
  