anchor.diff:         Mesh%ElemTable(ELEMENT_POINT)%Elements(Mesh%ElemTable(ELEMENT_POINT)%nelem)%ElemNodes(1) = P1
          Mesh%ElemTable(ELEMENT_POINT)%Elements(Mesh%ElemTable(ELEMENT_POINT)%nelem)%Xelement = ELEMENT_POINT
        ELSE
  +       ErrStat = ErrID_Fatal
          ErrMess = 'MeshConstructElement_1PT called for invalid element type'
  -       ErrStat = -1 !bjj: invalid ErrStat code. Use ErrID_Fatal or similar
        ENDIF
        RETURN
   
  
positive.code:      IF ( mesh_debug ) print*,'Called MeshConstructElement_6PT'
       ErrStat = ErrID_None
       ErrMess = ''
       ErrStat = 1
       ErrMess = 'MeshConstructElement_6PT not supported'
     END SUBROUTINE MeshConstructElement_6PT
  
  
negative.code:      LOGICAL IgSib
       INTEGER i, j, k
  
       ErrStat = ErrID_None 
  
       IF ( .NOT. Mesh%Initialized ) RETURN 
  
  