SELF_Mesh_3D_t Module

  1. SELF_Mesh_3D_t.f90
  2. SELF_Mesh_3D_t

Uses

  • module~~self_mesh_3d_t~~UsesGraph module~self_mesh_3d_t SELF_Mesh_3D_t module~self_mesh SELF_Mesh module~self_mesh_3d_t->module~self_mesh module~self_constants SELF_Constants module~self_mesh_3d_t->module~self_constants iso_c_binding iso_c_binding module~self_mesh_3d_t->iso_c_binding module~self_domaindecomposition~2 SELF_DomainDecomposition module~self_mesh_3d_t->module~self_domaindecomposition~2 HDF5 HDF5 module~self_mesh_3d_t->HDF5 module~self_supportroutines SELF_SupportRoutines module~self_mesh_3d_t->module~self_supportroutines module~self_lagrange~3 SELF_Lagrange module~self_mesh_3d_t->module~self_lagrange~3 module~self_quadrature SELF_Quadrature module~self_mesh_3d_t->module~self_quadrature module~self_hdf5 SELF_HDF5 module~self_mesh_3d_t->module~self_hdf5 module~self_mesh->module~self_constants module~self_mesh->iso_c_binding module~self_mesh->module~self_domaindecomposition~2 module~self_constants->iso_c_binding iso_fortran_env iso_fortran_env module~self_constants->iso_fortran_env module~self_domaindecomposition_t SELF_DomainDecomposition_t module~self_domaindecomposition~2->module~self_domaindecomposition_t module~self_supportroutines->module~self_constants module~self_supportroutines->iso_fortran_env module~self_lagrange~3->module~self_constants module~self_lagrange~3->iso_c_binding module~self_lagrange~3->iso_fortran_env module~self_lagrange_t SELF_Lagrange_t module~self_lagrange~3->module~self_lagrange_t module~self_quadrature->module~self_constants module~self_quadrature->iso_fortran_env module~self_hdf5->module~self_constants module~self_hdf5->HDF5 module~self_hdf5->iso_fortran_env mpi mpi module~self_hdf5->mpi module~self_lagrange_t->module~self_constants module~self_lagrange_t->iso_c_binding module~self_lagrange_t->HDF5 module~self_lagrange_t->module~self_supportroutines module~self_lagrange_t->module~self_quadrature module~self_lagrange_t->module~self_hdf5 module~self_lagrange_t->iso_fortran_env module~self_domaindecomposition_t->module~self_constants module~self_domaindecomposition_t->iso_c_binding module~self_domaindecomposition_t->module~self_supportroutines module~self_domaindecomposition_t->module~self_lagrange~3 module~self_domaindecomposition_t->mpi

Used by

  • module~~self_mesh_3d_t~~UsedByGraph module~self_mesh_3d_t SELF_Mesh_3D_t module~self_mesh_3d SELF_Mesh_3D module~self_mesh_3d->module~self_mesh_3d_t module~self_mesh_3d~2 SELF_Mesh_3D module~self_mesh_3d~2->module~self_mesh_3d_t module~self_mappedscalar_3d_t SELF_MappedScalar_3D_t module~self_mappedscalar_3d_t->module~self_mesh_3d module~self_geometry_3d SELF_Geometry_3D module~self_mappedscalar_3d_t->module~self_geometry_3d module~self_dgmodel3d SELF_DGModel3D module~self_dgmodel3d->module~self_mesh_3d module~self_dgmodel3d->module~self_geometry_3d module~self_dgmodel3d_t SELF_DGModel3D_t module~self_dgmodel3d->module~self_dgmodel3d_t module~self_geometry_3d->module~self_mesh_3d module~self_dgmodel3d_t->module~self_mesh_3d module~self_dgmodel3d_t->module~self_geometry_3d module~self_mappedscalar_3d SELF_MappedScalar_3D module~self_dgmodel3d_t->module~self_mappedscalar_3d module~self_mappedvector_3d~2 SELF_MappedVector_3D module~self_dgmodel3d_t->module~self_mappedvector_3d~2 module~self_mappedvector_3d_t SELF_MappedVector_3D_t module~self_mappedvector_3d_t->module~self_mesh_3d module~self_mappedvector_3d_t->module~self_geometry_3d module~self_esatmo3d~2 SELF_ESAtmo3D module~self_esatmo3d~2->module~self_mesh_3d module~self_esatmo3d~2->module~self_geometry_3d module~self_esatmo3d_t SELF_ESAtmo3D_t module~self_esatmo3d~2->module~self_esatmo3d_t module~self_ecdgmodel3d_t SELF_ECDGModel3D_t module~self_esatmo3d~2->module~self_ecdgmodel3d_t module~self_ecadvection3d~2 SELF_ECAdvection3D module~self_ecadvection3d~2->module~self_mesh_3d module~self_ecadvection3d~2->module~self_geometry_3d module~self_ecadvection3d~2->module~self_ecdgmodel3d_t module~self_ecadvection3d_t SELF_ECAdvection3D_t module~self_ecadvection3d~2->module~self_ecadvection3d_t module~self_mappedscalar_3d->module~self_mappedscalar_3d_t module~self_mappedscalar_3d~2 SELF_MappedScalar_3D module~self_mappedscalar_3d~2->module~self_mappedscalar_3d_t module~self_mappedtwopointvector_3d_t SELF_MappedTwoPointVector_3D_t module~self_mappedtwopointvector_3d_t->module~self_geometry_3d module~self_points SELF_Points module~self_points->module~self_geometry_3d module~self_points->module~self_mappedscalar_3d module~self_points_t SELF_Points_t module~self_points->module~self_points_t module~self_points_t->module~self_geometry_3d module~self_points_t->module~self_mappedscalar_3d module~self_dgmodel3d~2 SELF_DGModel3D module~self_dgmodel3d~2->module~self_dgmodel3d_t module~self_mappedvector_3d~2->module~self_mappedvector_3d_t module~self_mappedvector_3d SELF_MappedVector_3D module~self_mappedvector_3d->module~self_mappedvector_3d_t module~self_esatmo3d_t->module~self_mappedscalar_3d module~self_ecdgmodel3d SELF_ECDGModel3D module~self_esatmo3d_t->module~self_ecdgmodel3d module~self_mappedtwopointvector_3d~2 SELF_MappedTwoPointVector_3D module~self_mappedtwopointvector_3d~2->module~self_mappedtwopointvector_3d_t module~self_ecdgmodel3d_t->module~self_dgmodel3d~2 module~self_ecdgmodel3d_t->module~self_mappedtwopointvector_3d~2 module~self_points~2 SELF_Points module~self_points~2->module~self_points_t module~self_mappedtwopointvector_3d SELF_MappedTwoPointVector_3D module~self_mappedtwopointvector_3d->module~self_mappedtwopointvector_3d_t module~self_advection_diffusion_3d_t self_advection_diffusion_3d_t module~self_advection_diffusion_3d_t->module~self_dgmodel3d~2 module~self_lineareuler3d_t self_LinearEuler3D_t module~self_lineareuler3d_t->module~self_dgmodel3d~2 module~self_nulldgmodel3d_t self_NullDGModel3D_t module~self_nulldgmodel3d_t->module~self_dgmodel3d~2 module~self_esatmo3d SELF_ESAtmo3D module~self_esatmo3d->module~self_esatmo3d_t module~self_ecdgmodel3d~2 SELF_ECDGModel3D module~self_ecdgmodel3d~2->module~self_ecdgmodel3d_t module~self_advection_diffusion_3d self_advection_diffusion_3d module~self_advection_diffusion_3d->module~self_advection_diffusion_3d_t module~self_advection_diffusion_3d~2 self_advection_diffusion_3d module~self_advection_diffusion_3d~2->module~self_advection_diffusion_3d_t module~self_ecdgmodel3d->module~self_ecdgmodel3d_t module~self_lineareuler3d self_LinearEuler3D module~self_lineareuler3d->module~self_lineareuler3d_t module~self_lineareuler3d~2 self_LinearEuler3D module~self_lineareuler3d~2->module~self_lineareuler3d_t module~self_nulldgmodel3d~2 self_NullDGModel3D module~self_nulldgmodel3d~2->module~self_nulldgmodel3d_t module~self_nulldgmodel3d self_NullDGModel3D module~self_nulldgmodel3d->module~self_nulldgmodel3d_t module~self_ecadvection3d_t->module~self_ecdgmodel3d module~self_ecadvection3d SELF_ECAdvection3D module~self_ecadvection3d->module~self_ecadvection3d_t

Contents

Variables

Derived Types

Functions

Subroutines


Variables

TypeVisibilityAttributesNameInitial
integer, private :: CGNStoSELFflip(1:6,1:6,1:4)
integer, public, parameter:: selfSide3D_Bottom =1
integer, public, parameter:: selfSide3D_East =3
integer, public, parameter:: selfSide3D_North =4
integer, public, parameter:: selfSide3D_South =2
integer, public, parameter:: selfSide3D_Top =6
integer, public, parameter:: selfSide3D_West =5

Derived Types

type, public, extends(SEMMesh) :: 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 :: nBCs
integer, public :: nCornerNodes
integer, public :: nElem
integer, public :: nGeo
integer, public :: nGlobalElem
integer, public :: nMaterials =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
procedure, public :: Init => Init_Mesh3D_t
generic, public :: PeriodicStructuredMesh => UniformPeriodicMesh_Mesh3D_t
procedure, public :: Read_HOHQMesh => Read_HOHQMesh_Mesh3D_t
procedure, public :: Read_HOPr => Read_HOPr_Mesh3D_t
procedure, public :: RecalculateFlip => RecalculateFlip_Mesh3D_t
procedure, public :: ResetBoundaryConditionType => ResetBoundaryConditionType_Mesh3D_t
generic, public :: StructuredMesh => UniformStructuredMesh_Mesh3D_t
procedure, private :: UniformPeriodicMesh_Mesh3D_t
procedure, private :: UniformStructuredMesh_Mesh3D_t
procedure, public :: UpdateDevice => UpdateDevice_Mesh3D_t
procedure, public :: Write_Mesh => Write_Mesh3D_t

Functions

public pure function elementid(i, j, k, ti, tj, tk, nxpertile, nypertile, nzpertile, ntilex, ntiley, ntilez) result(eid)

Arguments

TypeIntentOptionalAttributesName
integer, intent(in) :: i
integer, intent(in) :: j
integer, intent(in) :: k
integer, intent(in) :: ti
integer, intent(in) :: tj
integer, intent(in) :: tk
integer, intent(in) :: nxpertile
integer, intent(in) :: nypertile
integer, intent(in) :: nzpertile
integer, intent(in) :: ntilex
integer, intent(in) :: ntiley
integer, intent(in) :: ntilez

Return Value integer

public pure function gid2eid(gx, gy, gz, nxpertile, nypertile, nzpertile, ntilex, ntiley, ntilez) result(eid)

Map a global element position (gx,gy,gz), with gx in [1,nX] etc., to the tile-ordered element id used throughout the structured mesh. This is the inverse of the (i,ti) decomposition: gx = i + nxpertile*(ti-1).

Arguments

TypeIntentOptionalAttributesName
integer, intent(in) :: gx
integer, intent(in) :: gy
integer, intent(in) :: gz
integer, intent(in) :: nxpertile
integer, intent(in) :: nypertile
integer, intent(in) :: nzpertile
integer, intent(in) :: ntilex
integer, intent(in) :: ntiley
integer, intent(in) :: ntilez

Return Value integer


Subroutines

public subroutine Free_Mesh3D_t(this)

Arguments

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

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

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

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

public subroutine RecalculateFlip_Mesh3D_t(this)

Arguments

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

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

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

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

public subroutine UpdateDevice_Mesh3D_t(this)

Arguments

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

public subroutine Write_Mesh3D_t(this, meshFile)

Arguments

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

public subroutine build_nodeCoords_for_hex(mesh, e, nGeo, allCorners, allFlags, allFaces, nodeXYZ)

Fill mesh%nodeCoords(:,:,:,:,e) for one hexahedral element by transfinite (Coons) interpolation of its six face grids. Faces flagged in the mesh file use the file's face-point grids; the remaining faces are bilinear patches of their four corner nodes evaluated at Chebyshev-Gauss-Lobatto parametric coordinates. Edge curves are extracted from the face grids (preferring a flagged face when an edge borders one flagged and one bilinear face) and the standard Boolean-sum formula x = Pxi + Peta + Pzeta - PxiPeta - PxiPzeta - PetaPzeta + PxiPeta*Pzeta combines face, edge, and corner contributions. For an element with all-straight faces this reduces to trilinear interpolation of the eight corners.

Arguments

TypeIntentOptionalAttributesName
class(Mesh3D_t), intent(inout) :: mesh
integer, intent(in) :: e
integer, intent(in) :: nGeo
integer, intent(in) :: allCorners(:,:)
integer, intent(in) :: allFlags(:,:)
real(kind=prec), intent(in) :: allFaces(:,:,:,:,:)
real(kind=prec), intent(in) :: nodeXYZ(:,:)