scale_mesh_cubedspheredom3d Module Reference

module FElib / Mesh / Cubed-sphere 3D domain More...

Data Types
type	meshcubedspheredom3d

Functions/Subroutines
subroutine	meshcubedspheredom3d_init (this, negx, negy, negz, rplanet, dom_zmin, dom_zmax, refelem, nlocalmeshperprc, nproc, myrank, fz, shallow_approx)

Detailed Description

module FElib / Mesh / Cubed-sphere 3D domain

Description

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

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

meshcubedspheredom3d_init()

subroutine scale_mesh_cubedspheredom3d::meshcubedspheredom3d_init	(	class(meshcubedspheredom3d), intent(inout)	this,
		integer, intent(in)	negx,
		integer, intent(in)	negy,
		integer, intent(in)	negz,
		real(rp), intent(in)	rplanet,
		real(rp), intent(in)	dom_zmin,
		real(rp), intent(in)	dom_zmax,
		type(hexahedralelement), intent(in), target	refelem,
		integer, intent(in)	nlocalmeshperprc,
		integer, intent(in), optional	nproc,
		integer, intent(in), optional	myrank,
		real(rp), dimension(negz+1), intent(in), optional	fz,
		logical, intent(in), optional	shallow_approx )

Definition at line 96 of file scale_mesh_cubedspheredom3d.F90.

    
    use scale_const, only: &
      pi => const_pi
    implicit none
 
    class(MeshCubedSphereDom3D), intent(inout) :: this
    integer, intent(in) :: NeGX
    integer, intent(in) :: NeGY
    integer, intent(in) :: NeGZ
    real(RP), intent(in) :: RPlanet
    real(RP), intent(in) :: dom_Zmin
    real(RP), intent(in) :: dom_zmax
    type(HexahedralElement), intent(in), target :: refElem
    integer, intent(in) :: NLocalMeshPerPrc
    integer, intent(in), optional :: nproc
    integer, intent(in), optional :: myrank
    real(RP), intent(in), optional :: FZ(NeGZ+1)
    logical, intent(in), optional :: shallow_approx
 
    integer :: k
    real(RP) :: dz
    !-----------------------------------------------------------------------------
 
    this%NeGX = negx
    this%NeGY = negy
    this%NeGZ = negz
 
    this%xmin_gl = - 0.25_rp * pi 
    this%xmax_gl = + 0.25_rp * pi 
    this%ymin_gl = - 0.25_rp * pi 
    this%ymax_gl = + 0.25_rp * pi 
    this%zmin_gl = dom_zmin
    this%zmax_gl = dom_zmax
    this%RPlanet = rplanet
    this%dom_vol = 4.0_rp / 3.0_rp * pi * ( ( dom_zmax + rplanet )**3 - ( dom_zmin + rplanet )**3 )
 
 
    !- 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
 
    !-
    if ( present(shallow_approx) ) then
      this%shallow_approx = shallow_approx
    else
      this%shallow_approx = .true.
    end if
 
    !--
    call meshbase3d_init( this, refelem, nlocalmeshperprc, 6, &
      nproc, myrank )
 
    !---
    call this%refElem2D%Init( this%refElem3D%PolyOrder_h, refelem%IsLumpedMatrix() )
    call this%mesh2D%Init( negx, negy, rplanet, this%refElem2D, nlocalmeshperprc, &
      nproc, myrank )
      
    !-- Modify the information of dimension for the cubed sphere mesh
    call this%SetDimInfo( meshbase3d_dimtypeid_x, "x", "1", "X-coordinate" )
    call this%SetDimInfo( meshbase3d_dimtypeid_y, "y", "1", "Y-coordinate" )
    call this%SetDimInfo( meshbase3d_dimtypeid_z, "z", "m", "Z-coordinate" )
    call this%SetDimInfo( meshbase3d_dimtypeid_xyz, "xyz", "1", "XYZ-coordinate" )
    call this%SetDimInfo( meshbase3d_dimtypeid_zt, "zt", "1", "XYZ-coordinate" )
    call this%SetDimInfo( meshbase3d_dimtypeid_xyzt, "xyzt", "1", "XYZ-coordinate" )
  
    return

References scale_localmesh_base::bctype_interior, scale_cubedsphere_coord_cnv::cubedspherecoordcnv_cs2lonlatpos(), scale_cubedsphere_coord_cnv::cubedspherecoordcnv_getmetric(), scale_mesh_base3d::meshbase3d_dimtypeid_x, scale_mesh_base3d::meshbase3d_dimtypeid_xyz, scale_mesh_base3d::meshbase3d_dimtypeid_xyzt, scale_mesh_base3d::meshbase3d_dimtypeid_y, scale_mesh_base3d::meshbase3d_dimtypeid_z, scale_mesh_base3d::meshbase3d_dimtypeid_zt, scale_mesh_base3d::meshbase3d_final(), scale_mesh_base3d::meshbase3d_init(), scale_mesh_base3d::meshbase3d_setgeometricinfo(), scale_mesh_cubedspheredom2d::meshcubedspheredom2d_check_division_params(), scale_mesh_cubedspheredom2d::meshcubedspheredom2d_setuplocaldom(), and scale_meshutil_cubedsphere3d::meshutilcubedsphere3d_buildglobalmap().