scale_mesh_cubedom3d Module Reference

module FElib / Mesh / Cubic 3D domain More...

Data Types
type	meshcubedom3d

Functions/Subroutines
subroutine	meshcubedom3d_init (this, negx, negy, negz, dom_xmin, dom_xmax, dom_ymin, dom_ymax, dom_zmin, dom_zmax, isperiodicx, isperiodicy, isperiodicz, refelem, nlocalmeshperprc, nprcx, nprcy, nproc, myrank, fz)
subroutine, public	meshcubedom3d_coord_conv (x, y, z, xx, xy, xz, yx, yy, yz, zx, zy, zz, vx, vy, vz, elem)

Detailed Description

module FElib / Mesh / Cubic 3D domain

Description

Manage mesh data of cubic 3D domain for element-based methods

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

meshcubedom3d_init()

subroutine scale_mesh_cubedom3d::meshcubedom3d_init	(	class(meshcubedom3d), intent(inout)	this,
		integer, intent(in)	negx,
		integer, intent(in)	negy,
		integer, intent(in)	negz,
		real(rp), intent(in)	dom_xmin,
		real(rp), intent(in)	dom_xmax,
		real(rp), intent(in)	dom_ymin,
		real(rp), intent(in)	dom_ymax,
		real(rp), intent(in)	dom_zmin,
		real(rp), intent(in)	dom_zmax,
		logical, intent(in)	isperiodicx,
		logical, intent(in)	isperiodicy,
		logical, intent(in)	isperiodicz,
		type(hexahedralelement), intent(in), target	refelem,
		integer, intent(in)	nlocalmeshperprc,
		integer, intent(in)	nprcx,
		integer, intent(in)	nprcy,
		integer, intent(in), optional	nproc,
		integer, intent(in), optional	myrank,
		real(rp), dimension(negz+1), intent(in), optional	fz )

Definition at line 95 of file scale_mesh_cubedom3d.F90.

    
    implicit none
 
    class(MeshCubeDom3D), intent(inout) :: this
    integer, intent(in) :: NeGX
    integer, intent(in) :: NeGY
    integer, intent(in) :: NeGZ
    real(RP), intent(in) :: dom_xmin
    real(RP), intent(in) :: dom_xmax
    real(RP), intent(in) :: dom_ymin
    real(RP), intent(in) :: dom_ymax
    real(RP), intent(in) :: dom_Zmin
    real(RP), intent(in) :: dom_zmax
    logical, intent(in) :: isPeriodicX
    logical, intent(in) :: isPeriodicY
    logical, intent(in) :: isPeriodicZ
    type(HexahedralElement), intent(in), target :: refElem
    integer, intent(in) :: NLocalMeshPerPrc
    integer, intent(in) :: NprcX
    integer, intent(in) :: NprcY
    integer, intent(in), optional :: nproc
    integer, intent(in), optional :: myrank
    real(RP), intent(in), optional :: FZ(NeGZ+1)
 
    integer :: k
    real(RP) :: dz
    !-----------------------------------------------------------------------------
    
    this%NeGX = negx
    this%NeGY = negy
    this%NeGZ = negz
 
    this%xmin_gl       = dom_xmin
    this%xmax_gl       = dom_xmax
    this%ymin_gl       = dom_ymin
    this%ymax_gl       = dom_ymax
    this%zmin_gl       = dom_zmin
    this%zmax_gl       = dom_zmax
    this%dom_vol = (this%xmax_gl - this%xmin_gl) * (this%ymax_gl - this%ymin_gl) * (this%zmax_gl - this%zmin_gl)  
 
    this%isPeriodicX = isperiodicx
    this%isPeriodicY = isperiodicy
    this%isPeriodicZ = isperiodicz
 
    this%NprcX = nprcx
    this%NprcY = nprcy
    this%NprcZ = 1
 
    !- Fz
    allocate( this%FZ(this%NeGZ+1) )
    if ( present(fz) ) then
      this%FZ(:) = fz(:)
    else
      this%FZ(1          ) = dom_zmin
      this%FZ(this%NeGZ+1) = dom_zmax
      dz = (dom_zmax - dom_zmin) / dble(this%NeGZ)
      do k=2, this%NeGZ
        this%FZ(k) = this%FZ(k-1) + dz
      end do
    end if
 
    !--
    call meshbase3d_init( this, refelem, nlocalmeshperprc, 6, &
                          nproc, myrank                       )
 
    !--- 2D mesh
 
    call this%refElem2D%Init( this%refElem3D%PolyOrder_h, refelem%IsLumpedMatrix() )
    
    this%mesh2D%isPeriodicX = isperiodicx
    this%mesh2D%isPeriodicY = isperiodicy
    this%mesh2D%NprcX = nprcx
    this%mesh2D%NprcY = nprcy
 
    call meshbase2d_init( this%mesh2D, this%refElem2D, nlocalmeshperprc,           &
                          nproc, myrank                                            )
 
    return

References scale_localmesh_base::bctype_interior, scale_mesh_base2d::meshbase2d_init(), scale_mesh_base3d::meshbase3d_final(), scale_mesh_base3d::meshbase3d_init(), scale_mesh_base3d::meshbase3d_setgeometricinfo(), meshcubedom3d_coord_conv(), scale_mesh_rectdom2d::meshrectdom2d_setuplocaldom(), scale_meshutil_3d::meshutil3d_buildglobalmap(), scale_meshutil_3d::meshutil3d_buildinteriormap(), scale_meshutil_3d::meshutil3d_genconnectivity(), scale_meshutil_3d::meshutil3d_gencubedomain(), and scale_meshutil_3d::meshutil3d_genpatchboundarymap().

meshcubedom3d_coord_conv()

subroutine, public scale_mesh_cubedom3d::meshcubedom3d_coord_conv	(	real(rp), dimension(elem%np), intent(out)	x,
		real(rp), dimension(elem%np), intent(out)	y,
		real(rp), dimension(elem%np), intent(out)	z,
		real(rp), dimension(elem%np), intent(out)	xx,
		real(rp), dimension(elem%np), intent(out)	xy,
		real(rp), dimension(elem%np), intent(out)	xz,
		real(rp), dimension(elem%np), intent(out)	yx,
		real(rp), dimension(elem%np), intent(out)	yy,
		real(rp), dimension(elem%np), intent(out)	yz,
		real(rp), dimension(elem%np), intent(out)	zx,
		real(rp), dimension(elem%np), intent(out)	zy,
		real(rp), dimension(elem%np), intent(out)	zz,
		real(rp), dimension(elem%nv), intent(in)	vx,
		real(rp), dimension(elem%nv), intent(in)	vy,
		real(rp), dimension(elem%nv), intent(in)	vz,
		type(elementbase3d), intent(in)	elem )

Definition at line 515 of file scale_mesh_cubedom3d.F90.

 
    implicit none
 
    type(ElementBase3D), intent(in) :: elem
    real(RP), intent(out) :: x(elem%Np), y(elem%Np), z(elem%Np)
    real(RP), intent(out) :: xX(elem%Np), xY(elem%Np), xZ(elem%Np)
    real(RP), intent(out) :: yX(elem%Np), yY(elem%Np), yZ(elem%Np)
    real(RP), intent(out) :: zX(elem%Np), zY(elem%Np), zZ(elem%Np)
    real(RP), intent(in) :: vx(elem%Nv), vy(elem%Nv), vz(elem%Nv)
 
    !-------------------------------------------------
    
    x(:) = vx(1) + 0.5_rp*(elem%x1(:) + 1.0_rp)*(vx(2) - vx(1))
    y(:) = vy(1) + 0.5_rp*(elem%x2(:) + 1.0_rp)*(vy(3) - vy(1))
    z(:) = vz(1) + 0.5_rp*(elem%x3(:) + 1.0_rp)*(vz(5) - vz(1))
 
    xx(:) = 0.5_rp*(vx(2) - vx(1)) !matmul(refElem%Dx1,mesh%x1(:,n))
    xy(:) = 0.0_rp                 !matmul(refElem%Dx2,mesh%x1(:,n))
    xz(:) = 0.0_rp                 !matmul(refElem%Dx3,mesh%x1(:,n))
    yx(:) = 0.0_rp                 !matmul(refElem%Dx1,mesh%x2(:,n))
    yy(:) = 0.5_rp*(vy(3) - vy(1)) !matmul(refElem%Dx2,mesh%x2(:,n))
    yz(:) = 0.0_rp                 !matmul(refElem%Dx3,mesh%x2(:,n))
    zx(:) = 0.0_rp                 !matmul(refElem%Dx1,mesh%x3(:,n))
    zy(:) = 0.0_rp                 !matmul(refElem%Dx2,mesh%x3(:,n))
    zz(:) = 0.5_rp*(vz(5) - vz(1)) !matmul(refElem%Dx3,mesh%x3(:,n))
 
    return

Referenced by meshcubedom3d_init(), and scale_mesh_cubedom3d::meshcubedom3d::set_geometric_with_vcoord().