scale_meshutil_cubedsphere3d Module Reference

module FElib / Mesh / utility for 3D cubed-sphere mesh More...

Functions/Subroutines
subroutine, public	meshutilcubedsphere3d_buildglobalmap (panelid_table, pi_table, pj_table, pk_table, tileid_map, tilefaceid_map, tilepanelid_map, ntile, nez)

Detailed Description

module FElib / Mesh / utility for 3D cubed-sphere mesh

Description

A module useful for generating 3D cubed-sphere mesh

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

meshutilcubedsphere3d_buildglobalmap()

subroutine, public scale_meshutil_cubedsphere3d::meshutilcubedsphere3d_buildglobalmap	(	integer, dimension(ntile), intent(out)	panelid_table,
		integer, dimension(ntile), intent(out)	pi_table,
		integer, dimension(ntile), intent(out)	pj_table,
		integer, dimension(ntile), intent(out)	pk_table,
		integer, dimension(6,ntile), intent(out)	tileid_map,
		integer, dimension(6,ntile), intent(out)	tilefaceid_map,
		integer, dimension(6,ntile), intent(out)	tilepanelid_map,
		integer, intent(in)	ntile,
		integer, intent(in)	nez )

Definition at line 46 of file scale_meshutil_cubedsphere3d.F90.

    
    ! use scale_prc, only: PRC_isMaster
    use scale_meshutil_3d, only: &
      meshutil3d_genconnectivity
    use scale_meshutil_cubedsphere2d, only: &
      meshutilcubedsphere2d_modifyconnectivity
    implicit none
 
    integer, intent(in) :: Ntile
    integer, intent(out) :: panelID_table(Ntile)
    integer, intent(out) :: pi_table(Ntile)
    integer, intent(out) :: pj_table(Ntile)
    integer, intent(out) :: pk_table(Ntile)
    integer, intent(out) :: tileID_map(6,Ntile)
    integer, intent(out) :: tileFaceID_map(6,Ntile)
    integer, intent(out) :: tilePanelID_map(6,Ntile)
    integer, intent(in) :: NeZ
 
    integer :: NtilePerPanel
    integer :: NeX, NeY, NvX, NvY, NvZ
    integer, allocatable :: nodesID_3d(:,:,:,:)
    integer, allocatable :: EToV(:,:)
    integer, allocatable :: EToE(:,:)
    integer, allocatable :: EToF(:,:)
    integer :: i, j, k, f
    integer :: panelID
    integer :: tileID, tileID_R
    integer :: counter
    
    integer :: pi_, pj_
    !-----------------------------------------------------------------------------
 
    ntileperpanel = ntile / 6
    ney = int( sqrt(dble(ntileperpanel)) )
    nex = ntileperpanel/ney
    nvx = nex + 1
    nvy = ney + 1
    nvz = nez + 1
    allocate( nodesid_3d(nvx,nvy,nvz,6) )
    allocate( etov(ntile,8), etoe(ntile,6), etof(ntile,6) )
 
    counter = 0
    do panelid = 1, 6
    do k = 1, nvz
    do j = 1, nvy
    do i = 1, nvx
      counter = counter + 1
      nodesid_3d(i,j,k,panelid) = counter
    end do
    end do
    end do
    end do
 
    !----
 
    tileid = 0
    do panelid = 1, 6
    do k = 1, nez
    do j = 1, ney
    do i = 1, nex
      tileid = tileid + 1
      panelid_table(tileid) = panelid
      pi_table(tileid) = i; pj_table(tileid) = j; pk_table(tileid) = k
      etov(tileid,:) = (/ nodesid_3d(i,j  ,k  ,panelid), nodesid_3d(i+1,j  ,k  ,panelid), &
                          nodesid_3d(i,j+1,k  ,panelid), nodesid_3d(i+1,j+1,k  ,panelid), &
                          nodesid_3d(i,j  ,k+1,panelid), nodesid_3d(i+1,j  ,k+1,panelid), &
                          nodesid_3d(i,j+1,k+1,panelid), nodesid_3d(i+1,j+1,k+1,panelid)  /)
    end do
    end do
    end do
    end do
 
    call meshutil3d_genconnectivity( etoe, etof, &
      etov, ntile, 6 )
    tileid_map(:,:) = transpose(etoe)
    tilefaceid_map(:,:) = transpose(etof)
    
    do tileid=1, ntile
    do f=1, 6
      tileid_r = tileid_map(f,tileid)
      tilepanelid_map(f,tileid) = panelid_table(tileid_r)
    end do
    end do
 
    call meshutilcubedsphere2d_modifyconnectivity( &
      tilepanelid_map, tileid_map, tilefaceid_map,          & ! (inout)
      panelid_table, pi_table, pj_table, nex, ney, ntile, 6 ) ! (in)
        
    return

References scale_meshutil_3d::meshutil3d_genconnectivity(), and scale_meshutil_cubedsphere2d::meshutilcubedsphere2d_modifyconnectivity().

Referenced by scale_mesh_cubedspheredom3d::meshcubedspheredom3d_init().