GetGlobalFaceCon method in ElemData_Class

[vickysharma0812]

Implement following method.

SUBROUTINE Elemdata_GetGlobalFaceCon(obj, globalFaceCon, localFaceCon)
  TYPE(Elemdata_), INTENT(INOUT) :: obj
  INTEGER(I4B), INTENT(INOUT) :: globalFaceCon(:, :)
  INTEGER(I4B), OPTIONAL, INTENT(INOUT) :: localFaceCon(:, :)

  INTEGER(I4B) :: tFaces, &
                  tNodes, localFaces0(4_I4B, TypeRefelemOpt%maxFaces), &
                  faceElemType(TypeRefelemOpt%maxFaces), &
                  tFaceNodes(TypeRefelemOpt%maxFaces), &
                  iface, face_temp(4), aint

#ifdef DEBUG_VER
  CHARACTER(*), PARAMETER :: myName = "Elemdata_GetGlobalFaceCon()"
#endif

  CALL RefElemGetGeoParam(elemType=obj%name, &
                          tFaces=tFaces, tNodes=tNodes, faceCon=localFaces0, &
                          faceOpt=1_I4B, faceElemType=faceElemType, &
                          tFaceNodes=tFaceNodes)

#ifdef DEBUG_VER
  CALL e%RaiseError(modName//'::'//myName//' - '// &
                    '[WIP ERROR] :: This routine is under development')
#endif

  DO iface = 1, tFaces
    aint = tFaceNodes(iface)
    face_temp(1:aint) = obj%globalNodes(localFaces0(1:aint, iface))

    CALL FaceShapeMetaData_Quadrangle(face=face_temp(1:aint), &
                                   sorted_face=globalFaceCon(1:aint, iface), &
                                      localFaces=localFaceCon(1:aint, iface))
  END DO

END SUBROUTINE Elemdata_GetGlobalFaceCon