Mesh3D_t Derived Type

  1. SELF_Mesh_3D_t.f90
  2. SELF_Mesh_3D_t
  3. Mesh3D_t

type, public, extends(SEMMesh) :: Mesh3D_t


Inherits

type~~mesh3d_t~~InheritsGraph type~mesh3d_t Mesh3D_t type~semmesh SEMMesh type~mesh3d_t->type~semmesh type~domaindecomposition DomainDecomposition type~semmesh->type~domaindecomposition decomp type~domaindecomposition_t DomainDecomposition_t type~domaindecomposition->type~domaindecomposition_t

Inherited by

type~~mesh3d_t~~InheritedByGraph type~mesh3d_t Mesh3D_t type~mesh3d Mesh3D type~mesh3d->type~mesh3d_t type~mesh3d~2 Mesh3D type~mesh3d~2->type~mesh3d_t type~dgmodel3d_t DGModel3D_t type~dgmodel3d_t->type~mesh3d mesh type~amrcontroller3d AMRController3D type~amrcontroller3d->type~mesh3d baseMesh, activeMesh type~dgmodel3d~2 DGModel3D type~dgmodel3d~2->type~dgmodel3d_t type~dgmodel3d DGModel3D type~dgmodel3d->type~dgmodel3d_t type~ecdgmodel3d_t ECDGModel3D_t type~ecdgmodel3d_t->type~dgmodel3d~2 type~nulldgmodel3d_t NullDGModel3D_t type~nulldgmodel3d_t->type~dgmodel3d~2 type~lineareuler3d_t LinearEuler3D_t type~lineareuler3d_t->type~dgmodel3d~2 type~advection_diffusion_3d_t advection_diffusion_3d_t type~advection_diffusion_3d_t->type~dgmodel3d~2 type~ecdgmodel3d ECDGModel3D type~ecdgmodel3d->type~ecdgmodel3d_t type~nulldgmodel3d NullDGModel3D type~nulldgmodel3d->type~nulldgmodel3d_t type~ecdgmodel3d~2 ECDGModel3D type~ecdgmodel3d~2->type~ecdgmodel3d_t type~nulldgmodel3d~2 NullDGModel3D type~nulldgmodel3d~2->type~nulldgmodel3d_t type~lineareuler3d LinearEuler3D type~lineareuler3d->type~lineareuler3d_t type~lineareuler3d~2 LinearEuler3D type~lineareuler3d~2->type~lineareuler3d_t type~advection_diffusion_3d advection_diffusion_3d type~advection_diffusion_3d->type~advection_diffusion_3d_t type~advection_diffusion_3d~2 advection_diffusion_3d type~advection_diffusion_3d~2->type~advection_diffusion_3d_t type~esatmo3d_t ESAtmo3D_t type~esatmo3d_t->type~ecdgmodel3d type~ecadvection3d_t ECAdvection3D_t type~ecadvection3d_t->type~ecdgmodel3d type~esatmo3d ESAtmo3D type~esatmo3d->type~esatmo3d_t type~ecadvection3d~2 ECAdvection3D type~ecadvection3d~2->type~ecadvection3d_t type~ecadvection3d ECAdvection3D type~ecadvection3d->type~ecadvection3d_t type~esatmo3d~2 ESAtmo3D type~esatmo3d~2->type~esatmo3d_t

Contents

Variables

Type-Bound Procedures

Source Code

Mesh3D_t

Components

TypeVisibilityAttributesNameInitial
character(len=255), public, allocatable:: BCNames(:)
integer, public, pointer, dimension(:,:):: BCType
integer, public, pointer, dimension(:,:):: CGNSCornerMap
integer, public, pointer, dimension(:,:):: CGNSSideMap
type(DomainDecomposition), public :: decomp
integer, public, pointer, dimension(:,:):: elemInfo
integer, public, allocatable:: elemMaterial(:)
integer, public, pointer, dimension(:,:,:,:):: globalNodeIDs
character(len=SELF_MESH_MATNAME_LENGTH), public, allocatable:: materialNames(:)
integer, public, pointer, dimension(:,:):: mortarInfo=> null()
integer, public :: nBCs
integer, public :: nCornerNodes
integer, public :: nElem
integer, public :: nGeo
integer, public :: nGlobalElem
integer, public :: nMaterials =0
integer, public :: nMortars =0
integer, public :: nNodes
integer, public :: nSides
integer, public :: nUniqueNodes
integer, public :: nUniqueSides
real(kind=prec), public, pointer, dimension(:,:,:,:,:):: nodeCoords
integer, public :: quadrature
integer, public, pointer, dimension(:,:,:):: sideInfo
integer, public, pointer, dimension(:,:):: sideMap

Type-Bound Procedures

procedure, public :: Free => Free_Mesh3D_t

  • public subroutine Free_Mesh3D_t(this)

    Arguments

    TypeIntentOptionalAttributesName
    class(Mesh3D_t), intent(inout) :: this

procedure, public :: Init => Init_Mesh3D_t

  • public subroutine Init_Mesh3D_t(this, nGeo, nElem, nSides, nNodes, nBCs)

    Arguments

    TypeIntentOptionalAttributesName
    class(Mesh3D_t), intent(inout) :: this
    integer, intent(in) :: nGeo
    integer, intent(in) :: nElem
    integer, intent(in) :: nSides
    integer, intent(in) :: nNodes
    integer, intent(in) :: nBCs

generic, public :: PeriodicStructuredMesh => UniformPeriodicMesh_Mesh3D_t

  • public subroutine UniformPeriodicMesh_Mesh3D_t(this, nxPerTile, nyPerTile, nzPerTile, nTileX, nTileY, nTileZ, dx, dy, dz, comm)

    Create a fully triply-periodic structured hexahedral mesh and store it in SELF's unstructured mesh format. Element geometry and ordering are identical to UniformStructuredMesh; the only difference is connectivity. The faces on the six domain boundaries are wired as interior faces whose neighbor is the element on the opposite side of the domain. This realises the triply periodic box T^3 = [0,Lx] x [0,Ly] x [0,Lz] required by, e.g., the Arnold-Beltrami-Childress (ABC) flow benchmark.

    Read more…

    Arguments

    TypeIntentOptionalAttributesName
    class(Mesh3D_t), intent(out) :: this
    integer, intent(in) :: nxPerTile
    integer, intent(in) :: nyPerTile
    integer, intent(in) :: nzPerTile
    integer, intent(in) :: nTileX
    integer, intent(in) :: nTileY
    integer, intent(in) :: nTileZ
    real(kind=prec), intent(in) :: dx
    real(kind=prec), intent(in) :: dy
    real(kind=prec), intent(in) :: dz
    integer, intent(in), optional :: comm

procedure, public :: Read_HOHQMesh => Read_HOHQMesh_Mesh3D_t

  • public subroutine Read_HOHQMesh_Mesh3D_t(this, meshFile, comm)

    Reader for HOHQMesh 3-D (hexahedral) text mesh files in the ISM and ISM-MM formats, as written by HOHQMesh's WriteISMHexMeshFile (Source/3DSource/Mesh3DOutputMethods.f90). Unlike the 2-D writer, the 3-D writer emits NO format header line: the first line is always the count line "nNodes nElems polyOrder". The two variants differ only in the per-element corner-node line: * ISM : 8 corner-node ids * ISM-MM : 8 corner-node ids followed by a material-name string The variant is auto-detected from the presence of the 9th token.

    Read more…

    Arguments

    TypeIntentOptionalAttributesName
    class(Mesh3D_t), intent(out) :: this
    character, intent(in) :: meshFile
    integer, intent(in), optional :: comm

procedure, public :: Read_HOPr => Read_HOPr_Mesh3D_t

  • public subroutine Read_HOPr_Mesh3D_t(this, meshFile, comm)

    Arguments

    TypeIntentOptionalAttributesName
    class(Mesh3D_t), intent(out) :: this
    character, intent(in) :: meshFile
    integer, intent(in), optional :: comm

procedure, public :: RecalculateFlip => RecalculateFlip_Mesh3D_t

procedure, public :: ResetBoundaryConditionType => ResetBoundaryConditionType_Mesh3D_t

  • public subroutine ResetBoundaryConditionType_Mesh3D_t(this, bcid)

    This method can be used to reset all of the boundary elements boundary condition type to the desired value.

    Read more…

    Arguments

    TypeIntentOptionalAttributesName
    class(Mesh3D_t), intent(inout) :: this
    integer, intent(in) :: bcid

procedure, public :: SimpleMortarMesh => SimpleMortarMesh_Mesh3D_t

  • public subroutine SimpleMortarMesh_Mesh3D_t(this, dx, bcids, comm, flips)

    Create the smallest 3D 2:1 nonconforming (mortar) mesh: one "big" element of size 2dx x 2dx x 2*dx whose east face is shared with the faces of four "small" dx x dx x dx elements stacked 2x2 in (y,z). The mesh is conforming everywhere except at the single mortar interface.

    Read more…

    Arguments

    TypeIntentOptionalAttributesName
    class(Mesh3D_t), intent(out) :: this
    real(kind=prec), intent(in) :: dx
    integer, intent(in) :: bcids(1:6)
    integer, intent(in), optional :: comm
    integer, intent(in), optional :: flips(1:4)

generic, public :: StructuredMesh => UniformStructuredMesh_Mesh3D_t

  • public subroutine UniformStructuredMesh_Mesh3D_t(this, nxPerTile, nyPerTile, nzPerTile, nTileX, nTileY, nTileZ, dx, dy, dz, bcids, comm)

    Create a structured mesh and store it in SELF's unstructured mesh format. The mesh is created in tiles of size (tnx,tny,tnz). Tiling is used to determine the element ordering.

    Read more…

    Arguments

    TypeIntentOptionalAttributesName
    class(Mesh3D_t), intent(out) :: this
    integer, intent(in) :: nxPerTile
    integer, intent(in) :: nyPerTile
    integer, intent(in) :: nzPerTile
    integer, intent(in) :: nTileX
    integer, intent(in) :: nTileY
    integer, intent(in) :: nTileZ
    real(kind=prec), intent(in) :: dx
    real(kind=prec), intent(in) :: dy
    real(kind=prec), intent(in) :: dz
    integer, intent(in) :: bcids(1:6)
    integer, intent(in), optional :: comm

procedure, private :: UniformPeriodicMesh_Mesh3D_t

  • public subroutine UniformPeriodicMesh_Mesh3D_t(this, nxPerTile, nyPerTile, nzPerTile, nTileX, nTileY, nTileZ, dx, dy, dz, comm)

    Create a fully triply-periodic structured hexahedral mesh and store it in SELF's unstructured mesh format. Element geometry and ordering are identical to UniformStructuredMesh; the only difference is connectivity. The faces on the six domain boundaries are wired as interior faces whose neighbor is the element on the opposite side of the domain. This realises the triply periodic box T^3 = [0,Lx] x [0,Ly] x [0,Lz] required by, e.g., the Arnold-Beltrami-Childress (ABC) flow benchmark.

    Read more…

    Arguments

    TypeIntentOptionalAttributesName
    class(Mesh3D_t), intent(out) :: this
    integer, intent(in) :: nxPerTile
    integer, intent(in) :: nyPerTile
    integer, intent(in) :: nzPerTile
    integer, intent(in) :: nTileX
    integer, intent(in) :: nTileY
    integer, intent(in) :: nTileZ
    real(kind=prec), intent(in) :: dx
    real(kind=prec), intent(in) :: dy
    real(kind=prec), intent(in) :: dz
    integer, intent(in), optional :: comm

procedure, private :: UniformStructuredMesh_Mesh3D_t

  • public subroutine UniformStructuredMesh_Mesh3D_t(this, nxPerTile, nyPerTile, nzPerTile, nTileX, nTileY, nTileZ, dx, dy, dz, bcids, comm)

    Create a structured mesh and store it in SELF's unstructured mesh format. The mesh is created in tiles of size (tnx,tny,tnz). Tiling is used to determine the element ordering.

    Read more…

    Arguments

    TypeIntentOptionalAttributesName
    class(Mesh3D_t), intent(out) :: this
    integer, intent(in) :: nxPerTile
    integer, intent(in) :: nyPerTile
    integer, intent(in) :: nzPerTile
    integer, intent(in) :: nTileX
    integer, intent(in) :: nTileY
    integer, intent(in) :: nTileZ
    real(kind=prec), intent(in) :: dx
    real(kind=prec), intent(in) :: dy
    real(kind=prec), intent(in) :: dz
    integer, intent(in) :: bcids(1:6)
    integer, intent(in), optional :: comm

procedure, public :: UpdateDevice => UpdateDevice_Mesh3D_t

procedure, public :: Write_Mesh => Write_Mesh3D_t

  • public subroutine Write_Mesh3D_t(this, meshFile)

    Arguments

    TypeIntentOptionalAttributesName
    class(Mesh3D_t), intent(inout) :: this
    character, intent(in) :: meshFile

Source Code

  type,extends(SEMMesh) :: Mesh3D_t
    integer,pointer,dimension(:,:,:) :: sideInfo
    real(prec),pointer,dimension(:,:,:,:,:) :: nodeCoords
    integer,pointer,dimension(:,:) :: elemInfo
    integer,pointer,dimension(:,:,:,:) :: globalNodeIDs
    integer,pointer,dimension(:,:) :: CGNSCornerMap
    integer,pointer,dimension(:,:) :: sideMap
    integer,pointer,dimension(:,:) :: CGNSSideMap
    integer,pointer,dimension(:,:) :: BCType
    character(LEN=255),allocatable :: BCNames(:)
    ! 2:1 nonconforming (mortar) interface bookkeeping.
    !
    ! A mortar interface joins one "big" element face to the four quarter-faces of its
    ! 2:1 refined neighbors ("small" faces). Faces that participate in a mortar carry
    ! sideInfo(3) = 0 (no conforming neighbor) and sideInfo(5) = 0 (no boundary
    ! condition), so all conforming-face machinery (SideExchange, RecalculateFlip,
    ! boundary-condition mapping) skips them; sideInfo(1) is set to the mortar index
    ! for reference.
    !
    ! The big face is split into four sub-faces ("quadrants") indexed in the big
    ! face's trace coordinates: q = kx + 2*(ky-1) with kx, ky in {1,2} the
    ! half-interval index along the face's first and second coordinate
    ! (1 = [-1,0], 2 = [0,1]).
    !
    ! mortarInfo(1:14,1:nMortars) is replicated on all ranks (element ids are global):
    !   1  - big element id ; 2 - big local face id
    !   3  - small element id on the sub-face covering big-face quadrant 1
    !   4  - 10*(small local face) + flip for that sub-face
    !   5,6   - same for quadrant 2
    !   7,8   - same for quadrant 3
    !   9,10  - same for quadrant 4
    !   11:14 - global side ids of sub-faces 1..4 (used as MPI message tags)
    ! The flip maps big-face quadrant coordinates to the small face's own coordinates
    ! (the same receiver-to-donor convention as sideInfo(4); see MortarFaceMap).
    integer :: nMortars = 0
    integer,pointer,dimension(:,:) :: mortarInfo => null()
    ! Material tracking: every element has an integer material id
    ! indexing into materialNames. Single-material constructors and
    ! readers (HOPr, structured, periodic) leave nMaterials = 1 with
    ! the name "default". The HOHQMesh ISM-MM reader populates the
    ! table with the material strings from the .mesh file.
    integer :: nMaterials = 0
    integer,allocatable :: elemMaterial(:)
    character(LEN=SELF_MESH_MATNAME_LENGTH),allocatable :: materialNames(:)

  contains

    procedure,public :: Init => Init_Mesh3D_t
    procedure,public :: Free => Free_Mesh3D_t
    procedure,public :: UpdateDevice => UpdateDevice_Mesh3D_t

    generic,public :: StructuredMesh => UniformStructuredMesh_Mesh3D_t
    procedure,private :: UniformStructuredMesh_Mesh3D_t

    procedure,public :: SimpleMortarMesh => SimpleMortarMesh_Mesh3D_t

    generic,public :: PeriodicStructuredMesh => UniformPeriodicMesh_Mesh3D_t
    procedure,private :: UniformPeriodicMesh_Mesh3D_t

    procedure,public :: Read_HOPr => Read_HOPr_Mesh3D_t
    procedure,public :: Read_HOHQMesh => Read_HOHQMesh_Mesh3D_t

    procedure,public :: ResetBoundaryConditionType => ResetBoundaryConditionType_Mesh3D_t

    procedure,public :: Write_Mesh => Write_Mesh3D_t

    procedure,public :: RecalculateFlip => RecalculateFlip_Mesh3D_t

  endtype Mesh3D_t