scale_mesh_cubedspheredom2d Module Reference

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

Data Types
type	meshcubedspheredom2d
	Derived type to manage a cubed-sphere 2D computational domain. More...

Functions/Subroutines
subroutine	meshcubedspheredom2d_init (this, negx, negy, rplanet, refelem, nlocalmeshperprc, nproc, myrank)
subroutine, public	meshcubedspheredom2d_check_division_params (nprcx_lc, nprcy_lc, prc_num, local_mesh_num_global, call_prc_abort)
subroutine, public	meshcubedspheredom2d_setuplocaldom (lcmesh, tileid, panelid, i, j, nprcx, nprcy, dom_xmin, dom_xmax, dom_ymin, dom_ymax, planet_radius, nex, ney)

Detailed Description

module FElib / Mesh / Cubed-sphere 2D domain

Description

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

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

meshcubedspheredom2d_init()

subroutine scale_mesh_cubedspheredom2d::meshcubedspheredom2d_init	(	class(meshcubedspheredom2d), intent(inout)	this,
		integer, intent(in)	negx,
		integer, intent(in)	negy,
		real(rp), intent(in)	rplanet,
		type(quadrilateralelement), intent(in), target	refelem,
		integer, intent(in)	nlocalmeshperprc,
		integer, intent(in), optional	nproc,
		integer, intent(in), optional	myrank )

Definition at line 84 of file scale_mesh_cubedspheredom2d.F90.

    
    use scale_const, only: &
      pi => const_pi
    implicit none
 
    class(MeshCubedSphereDom2D), intent(inout) :: this
    integer, intent(in) :: NeGX
    integer, intent(in) :: NeGY
    real(RP), intent(in) :: RPlanet
    type(QuadrilateralElement), intent(in), target :: refElem
    integer, intent(in) :: NLocalMeshPerPrc
    integer, intent(in), optional :: nproc
    integer, intent(in), optional :: myrank
    !-----------------------------------------------------------------------------
 
    this%NeGX = negx
    this%NeGY = negy
 
    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%RPlanet = rplanet
    this%dom_vol = 4.0_rp * pi * rplanet**2
 
 
    call meshbase2d_init( this, refelem, nlocalmeshperprc, &
      nproc, myrank )
 
    call this%SetDimInfo( meshbase2d_dimtypeid_x, "x", "1", "X-coordinate" )
    call this%SetDimInfo( meshbase2d_dimtypeid_y, "y", "1", "Y-coordinate" )
    call this%SetDimInfo( meshbase2d_dimtypeid_xy, "xy", "1", "XY-coordinate" )
    call this%SetDimInfo( meshbase2d_dimtypeid_xyt, "xyt", "1", "XY-coordinate" )
  
    return

References scale_mesh_base2d::meshbase2d_dimtypeid_x, scale_mesh_base2d::meshbase2d_dimtypeid_xy, scale_mesh_base2d::meshbase2d_dimtypeid_xyt, scale_mesh_base2d::meshbase2d_dimtypeid_y, scale_mesh_base2d::meshbase2d_final(), scale_mesh_base2d::meshbase2d_init(), meshcubedspheredom2d_check_division_params(), and meshcubedspheredom2d_setuplocaldom().

meshcubedspheredom2d_check_division_params()

subroutine, public scale_mesh_cubedspheredom2d::meshcubedspheredom2d_check_division_params	(	integer, intent(out)	nprcx_lc,
		integer, intent(out)	nprcy_lc,
		integer, intent(in)	prc_num,
		integer, intent(in)	local_mesh_num_global,
		logical, intent(in), optional	call_prc_abort )

Definition at line 200 of file scale_mesh_cubedspheredom2d.F90.

    
    use scale_prc, only: prc_abort
    implicit none
 
    integer, intent(in) :: PRC_NUM
    integer, intent(in) :: LOCAL_MESH_NUM_global
    integer, intent(out) :: NprcX_lc
    integer, intent(out) :: NprcY_lc
    logical, intent(in), optional :: call_prc_abort
 
    integer :: TILE_NUM_PER_PANEL
    logical :: call_prc_abort_
    !-----------------------------------------------------------------------------
 
    if (present(call_prc_abort)) then
      call_prc_abort_ = call_prc_abort
    else
      call_prc_abort_ = .false.
    end if
 
    if ( mod(local_mesh_num_global, 6) /= 0 ) then
      log_error("MeshCubedSphereDom2D_division_params",*) "The total number of local mesh must be a multiple of 6. Check!"
      if (call_prc_abort_) call prc_abort
    end if
 
    tile_num_per_panel = local_mesh_num_global / 6
 
    if ( prc_num <= 6 ) then
      if ( ( prc_num == 1 ) .or. &
          ( prc_num <= 6 .and. (mod(prc_num,2)==0 .or. mod(prc_num,3)==0)) ) then
        nprcx_lc = 1; nprcy_lc = 1
      else
        log_error("MeshCubedSphereDom2D_division_params",*) "The number of proceses is inappropriate. Check!"
        if (call_prc_abort_) call prc_abort
      end if
    else
      if ( mod(prc_num,6) == 0 ) then
        nprcx_lc = int(sqrt(dble(tile_num_per_panel)))
        nprcy_lc = tile_num_per_panel / nprcx_lc
        if ( nprcx_lc /= nprcy_lc ) then
          log_error("MeshCubedSphereDom2D_division_params",*) "The number of proceses is inappropriate. Check!"
          if (call_prc_abort_) call prc_abort  
        end if
      else
        log_error("MeshCubedSphereDom2D_division_params",*) "The number of proceses is inappropriate. Check!"
        if (call_prc_abort_) call prc_abort        
      end if
    end if
 
    return

Referenced by scale_mesh_cubedspheredom2d::meshcubedspheredom2d::assigndomid(), meshcubedspheredom2d_init(), and scale_mesh_cubedspheredom3d::meshcubedspheredom3d_init().

meshcubedspheredom2d_setuplocaldom()

subroutine, public scale_mesh_cubedspheredom2d::meshcubedspheredom2d_setuplocaldom	(	type(localmesh2d), intent(inout)	lcmesh,
		integer, intent(in)	tileid,
		integer, intent(in)	panelid,
		integer, intent(in)	i,
		integer, intent(in)	j,
		integer, intent(in)	nprcx,
		integer, intent(in)	nprcy,
		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)	planet_radius,
		integer, intent(in)	nex,
		integer, intent(in)	ney )

Definition at line 256 of file scale_mesh_cubedspheredom2d.F90.

 
    use scale_meshutil_cubedsphere2d, only: &
      meshutilcubedsphere2d_genconnectivity,   &
      meshutilcubedsphere2d_genrectdomain,     &
      meshutilcubedsphere2d_buildinteriormap,  &
      meshutilcubedsphere2d_genpatchboundarymap
 
    use scale_cubedsphere_coord_cnv, only: &
      cubedspherecoordcnv_getmetric,      &
      cubedspherecoordcnv_cs2lonlatpos
    
    use scale_localmesh_base, only: bctype_interior
 
    implicit none
      
    type(LocalMesh2D), intent(inout) :: lcmesh
    integer, intent(in) :: tileID
    integer, intent(in) :: panelID
    integer, intent(in) :: i, j
    integer, intent(in) :: NprcX, NprcY
    real(RP), intent(in) :: dom_xmin, dom_xmax
    real(RP), intent(in) :: dom_ymin, dom_ymax
    real(RP), intent(in) :: planet_radius
    integer, intent(in) :: NeX, NeY
    
    class(ElementBase2D), pointer :: elem
    real(RP) :: delx, dely
    integer :: ke
 
    real(RP), allocatable :: gam(:,:)
    !-----------------------------------------------------------------------------
 
    elem => lcmesh%refElem2D
    
    lcmesh%tileID = tileid
    lcmesh%panelID = panelid
    !$acc update device(lcmesh%tileID, lcmesh%panelID)
        
    !--
    lcmesh%Ne  = nex * ney
    lcmesh%Nv  = (nex + 1)*(ney + 1)
    lcmesh%NeS = 1
    lcmesh%NeE = lcmesh%Ne
    lcmesh%NeA = lcmesh%Ne + 2*(nex + ney)
    !$acc update device(lcmesh%Ne, lcmesh%Nv, lcmesh%NeS, lcmesh%NeE, lcmesh%NeA)
 
    lcmesh%NeX = nex
    lcmesh%NeY = ney
    !$acc update device(lcmesh%NeX, lcmesh%NeY)
 
    !--
    delx = ( dom_xmax - dom_xmin ) / dble(nprcx)
    dely = ( dom_ymax - dom_ymin ) / dble(nprcy)
 
    lcmesh%xmin = dom_xmin + (i-1)*delx
    lcmesh%xmax = dom_xmin +  i   *delx
    lcmesh%ymin = dom_ymin + (j-1)*dely
    lcmesh%ymax = dom_ymin +  j   *dely
    !$acc update device(lcmesh%xmin, lcmesh%xmax, lcmesh%ymin, lcmesh%ymax)
 
    !--
    allocate( lcmesh%pos_ev(lcmesh%Nv,2) )
    allocate( lcmesh%EToV(lcmesh%Ne,elem%Nv) )
    allocate( lcmesh%EToE(lcmesh%Ne,elem%Nfaces) )
    allocate( lcmesh%EToF(lcmesh%Ne,elem%Nfaces) )
    allocate( lcmesh%BCType(lcmesh%refElem%Nfaces,lcmesh%Ne) )
    allocate( lcmesh%VMapM(elem%NfpTot, lcmesh%Ne) )
    allocate( lcmesh%VMapP(elem%NfpTot, lcmesh%Ne) )
    allocate( lcmesh%MapM(elem%NfpTot, lcmesh%Ne) )
    allocate( lcmesh%MapP(elem%NfpTot, lcmesh%Ne) )
    !$acc enter data create(lcmesh%pos_ev, lcmesh%EToV, lcmesh%EToE, lcmesh%EToF, lcmesh%BCType, &
    !$acc   lcmesh%VMapM, lcmesh%VMapP, lcmesh%MapM, lcmesh%MapP)
 
    lcmesh%BCType(:,:) = bctype_interior
    !$acc update device(lcmesh%BCType)
 
    !----
 
    call meshutilcubedsphere2d_genrectdomain( lcmesh%pos_ev, lcmesh%EToV, & ! (out)
      lcmesh%NeX, lcmesh%xmin, lcmesh%xmax,                               & ! (in)
      lcmesh%NeY, lcmesh%ymin, lcmesh%ymax          )                       ! (in)
    !$acc update device(lcmesh%pos_ev, lcmesh%EToV)
    
    !---
    call meshbase2d_setgeometricinfo(lcmesh, meshcubedspheredom2d_coord_conv, meshcubedspheredom2d_calc_normal )
    
    call cubedspherecoordcnv_getmetric( &
      lcmesh%pos_en(:,:,1), lcmesh%pos_en(:,:,2), elem%Np * lcmesh%Ne, planet_radius, & ! (in)
      lcmesh%G_ij, lcmesh%GIJ, lcmesh%Gsqrt(:,lcmesh%NeS:lcmesh%NeE)                  ) ! (out)
    !$acc update host(lcmesh%G_ij, lcmesh%GIJ)
 
    allocate( gam(elem%Np,lcmesh%Ne) )
    gam(:,:) = 1.0_rp
    !$acc data copyin(gam)
 
    call cubedspherecoordcnv_cs2lonlatpos( &
      lcmesh%panelID, lcmesh%pos_en(:,:,1), lcmesh%pos_en(:,:,2), gam(:,:), & ! (in)
      lcmesh%Ne * lcmesh%refElem2D%Np,                                      & ! (in)
      lcmesh%lon(:,:), lcmesh%lat(:,:)                                      ) ! (out)
    !$acc update host(lcmesh%lon, lcmesh%lat)
 
    !$acc end data
    
    !---
 
    call meshutilcubedsphere2d_genconnectivity( lcmesh%EToE, lcmesh%EToF, & ! (out)
      lcmesh%EToV, lcmesh%Ne, elem%Nfaces )                                 ! (in)
    !$acc update device(lcmesh%EToE, lcmesh%EToF)
 
    !---
    call meshutilcubedsphere2d_buildinteriormap( &
      lcmesh%VmapM, lcmesh%VMapP, lcmesh%MapM, lcmesh%MapP,                &
      lcmesh%pos_en, lcmesh%pos_ev, lcmesh%EToE, lcmesh%EtoF, lcmesh%EtoV, &
      elem%Fmask, lcmesh%Ne, elem%Np, elem%Nfp, elem%Nfaces, lcmesh%Nv     )
 
    call meshutilcubedsphere2d_genpatchboundarymap( &
      lcmesh%VMapB, lcmesh%MapB, lcmesh%VMapP,                           &
      lcmesh%pos_en, lcmesh%xmin, lcmesh%xmax, lcmesh%ymin, lcmesh%ymax, &
      elem%Fmask, lcmesh%Ne, elem%Np, elem%Nfp, elem%Nfaces, lcmesh%Nv   )
    !$acc update device(lcmesh%VMapM, lcmesh%VMapP, lcmesh%MapM, lcmesh%MapP)
    !$acc enter data copyin(lcmesh%VMapB, lcmesh%MapB)
    
    !--
    call fill_halo_metric( lcmesh%Gsqrt, &
      lcmesh%VMapM, lcmesh%VMapP, lcmesh, elem )
    !$acc update host(lcmesh%Gsqrt)
 
    return

References scale_localmesh_base::bctype_interior, scale_cubedsphere_coord_cnv::cubedspherecoordcnv_cs2lonlatpos(), scale_cubedsphere_coord_cnv::cubedspherecoordcnv_getmetric(), scale_mesh_base2d::meshbase2d_setgeometricinfo(), and scale_meshutil_cubedsphere2d::meshutilcubedsphere2d_buildglobalmap().

Referenced by scale_mesh_cubedspheredom2d::meshcubedspheredom2d::assigndomid(), meshcubedspheredom2d_init(), and scale_mesh_cubedspheredom3d::meshcubedspheredom3d_init().