block_obj_SubDivideBlock Subroutine

private pure subroutine block_obj_SubDivideBlock(coord, Nb, minlo, maxhi, sublo, subhi)

Computes the bounds of the sub-block form the bounds of the parent block. Each sub-block gets about the same number of grid points in each direction.

Type Bound

block_obj

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: coord(3)

Coordinates of MPI rank on MPI grid
integer, intent(in) :: Nb(3)

Number of blocks in each direction
integer, intent(in) :: minlo(3)

Global lo bounds
integer, intent(in) :: maxhi(3)

Global hi bounds
integer, intent(out) :: sublo(3)

lo bounds of sub-block
integer, intent(out) :: subhi(3)

hi bounds of sub-block

Called by

proc~~block_obj_subdivideblock~~CalledByGraph proc~block_obj_subdivideblock block_obj%block_obj_SubDivideBlock proc~block_obj_partition block_obj%block_obj_Partition proc~block_obj_partition->proc~block_obj_subdivideblock proc~cdifs_obj_preparesolverblock cdifs_obj_PrepareSolverBlock proc~cdifs_obj_preparesolverblock->proc~block_obj_partition proc~grans_obj_preparesolverblock grans_obj_PrepareSolverBlock proc~grans_obj_preparesolverblock->proc~block_obj_partition proc~cdifs_obj_preparesolver cdifs_obj_PrepareSolver proc~cdifs_obj_preparesolver->proc~cdifs_obj_preparesolverblock proc~grans_obj_preparesolver grans_obj_PrepareSolver proc~grans_obj_preparesolver->proc~grans_obj_preparesolverblock interface~cdifs_obj_preparesolver cdifs_obj%cdifs_obj_PrepareSolver interface~cdifs_obj_preparesolver->proc~cdifs_obj_preparesolver interface~grans_obj_preparesolver grans_obj%grans_obj_PrepareSolver interface~grans_obj_preparesolver->proc~grans_obj_preparesolver

Source Code

    pure subroutine block_obj_SubDivideBlock(coord,Nb,minlo,maxhi,sublo,subhi)
      !> Computes the bounds of the sub-block form
      ! the bounds of the parent block. Each sub-block
      ! gets about the same number of grid points in each
      ! direction.
      implicit none
      integer, intent(in)  :: coord(3)                                         !! Coordinates of MPI rank on MPI grid
      integer, intent(in)  :: Nb(3)                                            !! Number of blocks in each direction
      integer, intent(in)  :: minlo(3)                                         !! Global lo bounds
      integer, intent(in)  :: maxhi(3)                                         !! Global hi bounds
      integer, intent(out) :: sublo(3)                                         !! lo bounds of sub-block
      integer, intent(out) :: subhi(3)                                         !! hi bounds of sub-block
      ! Work variables
      integer :: q(3)
      integer :: r(3)
      integer :: dir

      ! Number of grid points per sub-block in each direction
      q  = (maxhi - minlo + 1)/Nb
      r  = mod(maxhi - minlo + 1,Nb)

      ! Set the lo and hi bounds
      do dir=1,3
        if (coord(dir).le.r(dir)) then
          sublo(dir) = minlo(dir)   + (coord(dir)-1)*(q(dir)+1)
          subhi(dir) = sublo(dir)-1 + q(dir)+1
        else
          sublo(dir) = minlo(dir)   + (coord(dir)-1)*q(dir)+r(dir)
          subhi(dir) = sublo(dir)-1 + q(dir)
        end if
      end do

      return
    end subroutine block_obj_SubDivideBlock