SPLT_SplitTopology.F90

Go to the documentation of this file.

! *********************************************************************
! * Rocstar Simulation Suite                                          *
! * Copyright@2015, Illinois Rocstar LLC. All rights reserved.        *
! *                                                                   *
! * Illinois Rocstar LLC                                              *
! * Champaign, IL                                                     *
! * www.illinoisrocstar.com                                           *
! * sales@illinoisrocstar.com                                         *
! *                                                                   *
! * License: See LICENSE file in top level of distribution package or *
! * http://opensource.org/licenses/NCSA                               *
! *********************************************************************
! *********************************************************************
! * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,   *
! * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES   *
! * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND          *
! * NONINFRINGEMENT.  IN NO EVENT SHALL THE CONTRIBUTORS OR           *
! * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER       *
! * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,   *
! * Arising FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE    *
! * USE OR OTHER DEALINGS WITH THE SOFTWARE.                          *
! *********************************************************************
!******************************************************************************
!
! Purpose: generate new topology for the splitted grid.
!
! Description: none.
!
! Input: splitDirection = direction (i,j,k) in which the grid is to be splitted
!        regionsOld     = old topology.
!
! Output: regionsNew = new topology.
!
! Notes: none.
!
!******************************************************************************
!
! $Id: SPLT_SplitTopology.F90,v 1.4 2008/12/06 08:44:51 mtcampbe Exp $
!
! Copyright: (c) 2001 by the University of Illinois
!
!******************************************************************************

SUBROUTINE splittopology( splitDirection,regionsOld,regionsNew )

  USE moddatatypes
  USE modbndpatch, ONLY   : t_patch
  USE moddatastruct, ONLY : t_region
  USE modglobal, ONLY     : t_global
  USE modgrid, ONLY       : t_grid
  USE splt_modinterfaces, ONLY : copypatchdata
  USE moderror
  USE modparameters
  IMPLICIT NONE

! ... parameters
  INTEGER :: splitdirection

  TYPE(t_region), POINTER :: regionsold(:), regionsnew(:)

! ... loop variables
  INTEGER :: ireg, ipatchold

! ... local variables
  INTEGER :: lbound, l1beg, l1end, l2beg, l2end
  INTEGER :: ibeg, iend, jbeg, jend, kbeg, kend, ipatch

  TYPE(t_global), POINTER :: globalnew
  TYPE(t_grid) , POINTER  :: gridold, gridnew
  TYPE(t_patch), POINTER  :: patchold, patchnew

!******************************************************************************

  globalnew => regionsnew(1)%global

  CALL registerfunction( globalnew,'SplitTopology',&
  'SPLT_SplitTopology.F90' )

! copy old patches ------------------------------------------------------------

  gridold => regionsold(1)%levels(1)%grid

  DO ipatchold=1,regionsold(1)%nPatches
    patchold => regionsold(1)%levels(1)%patches(ipatchold)
    lbound   = patchold%lbound
    l1beg    = patchold%l1beg
    l1end    = patchold%l1end
    l2beg    = patchold%l2beg
    l2end    = patchold%l2end
    ibeg     = 1                ! offsets of new regions wrp. to old grid
    jbeg     = 1
    kbeg     = 1

! - loop over all new regions

    DO ireg=1,globalnew%nRegions

      gridnew => regionsnew(ireg)%levels(1)%grid
      iend    = ibeg + gridnew%ipc - 1
      jend    = jbeg + gridnew%jpc - 1
      kend    = kbeg + gridnew%kpc - 1

! --- i=const faces

      IF (lbound==1 .AND. ibeg==1) THEN
        IF ((l1beg<=jend .AND. l1end>=jbeg) .AND. &
            (l2beg<=kend .AND. l2end>=kbeg)) THEN
          regionsnew(ireg)%nPatches = regionsnew(ireg)%nPatches + 1
          CALL copypatchdata( regionsnew,patchold,splitdirection,ireg, &
                              jbeg,kbeg,gridnew%jpc,gridnew%kpc )
        ENDIF

      ELSE IF (lbound==2 .AND. iend==gridold%ipc) THEN
        IF ((l1beg<=jend .AND. l1end>=jbeg) .AND. &
            (l2beg<=kend .AND. l2end>=kbeg)) THEN
          regionsnew(ireg)%nPatches = regionsnew(ireg)%nPatches + 1
          CALL copypatchdata( regionsnew,patchold,splitdirection,ireg, &
                              jbeg,kbeg,gridnew%jpc,gridnew%kpc )
        ENDIF

! --- j=const. faces

      ELSE IF (lbound==3 .AND. jbeg==1) THEN
        IF ((l1beg<=kend .AND. l1end>=kbeg) .AND. &
            (l2beg<=iend .AND. l2end>=ibeg)) THEN
          regionsnew(ireg)%nPatches = regionsnew(ireg)%nPatches + 1
          CALL copypatchdata( regionsnew,patchold,splitdirection,ireg, &
                              kbeg,ibeg,gridnew%kpc,gridnew%ipc )
        ENDIF

      ELSE IF (lbound==4 .AND. jend==gridold%jpc) THEN
        IF ((l1beg<=kend .AND. l1end>=kbeg) .AND. &
            (l2beg<=iend .AND. l2end>=ibeg)) THEN
          regionsnew(ireg)%nPatches = regionsnew(ireg)%nPatches + 1
          CALL copypatchdata( regionsnew,patchold,splitdirection,ireg, &
                              kbeg,ibeg,gridnew%kpc,gridnew%ipc )
        ENDIF

! --- k=const. faces

      ELSE IF (lbound==5 .AND. kbeg==1) THEN
        IF ((l1beg<=iend .AND. l1end>=ibeg) .AND. &
            (l2beg<=jend .AND. l2end>=jbeg)) THEN
          regionsnew(ireg)%nPatches = regionsnew(ireg)%nPatches + 1
          CALL copypatchdata( regionsnew,patchold,splitdirection,ireg, &
                              ibeg,jbeg,gridnew%ipc,gridnew%jpc )
        ENDIF

      ELSE IF (lbound==6 .AND. kend==gridold%kpc) THEN
        IF ((l1beg<=iend .AND. l1end>=ibeg) .AND. &
            (l2beg<=jend .AND. l2end>=jbeg)) THEN
          regionsnew(ireg)%nPatches = regionsnew(ireg)%nPatches + 1
          CALL copypatchdata( regionsnew,patchold,splitdirection,ireg, &
                              ibeg,jbeg,gridnew%ipc,gridnew%jpc )
        ENDIF
      ENDIF    ! lbound

! --- increase offsets

      IF (splitdirection == 1) THEN
        ibeg = ibeg + gridnew%ipc
      ELSE IF (splitdirection == 2) THEN
        jbeg = jbeg + gridnew%jpc
      ELSE
        kbeg = kbeg + gridnew%kpc
      ENDIF

    ENDDO  ! iReg

  ENDDO    ! iPatchOld

! generate connections between regions ----------------------------------------
! faces 2, 4, 6

  DO ireg=1,globalnew%nRegions-1
    regionsnew(ireg)%nPatches = regionsnew(ireg)%nPatches + 1
    ipatch             = regionsnew(ireg)%nPatches
    patchnew           => regionsnew(ireg)%levels(1)%patches(ipatch)
    gridnew            => regionsnew(ireg)%levels(1)%grid
    patchnew%bcType    = bc_regionconf
    patchnew%bcCoupled = 0
    patchnew%align     = .true.
    patchnew%srcRegion = ireg + 1

    IF (splitdirection == 1) THEN
      patchnew%lbound    = 2
      patchnew%l1beg     = 1
      patchnew%l1end     = gridnew%jpc
      patchnew%l2beg     = 1
      patchnew%l2end     = gridnew%kpc
      patchnew%srcLbound = 1
      patchnew%srcL1beg  = 1
      patchnew%srcL1end  = gridnew%jpc
      patchnew%srcL2beg  = 1
      patchnew%srcL2end  = gridnew%kpc
    ELSE IF (splitdirection == 2) THEN
      patchnew%lbound    = 4
      patchnew%l1beg     = 1
      patchnew%l1end     = gridnew%kpc
      patchnew%l2beg     = 1
      patchnew%l2end     = gridnew%ipc
      patchnew%srcLbound = 3
      patchnew%srcL1beg  = 1
      patchnew%srcL1end  = gridnew%kpc
      patchnew%srcL2beg  = 1
      patchnew%srcL2end  = gridnew%ipc
    ELSE
      patchnew%lbound    = 6
      patchnew%l1beg     = 1
      patchnew%l1end     = gridnew%ipc
      patchnew%l2beg     = 1
      patchnew%l2end     = gridnew%jpc
      patchnew%srcLbound = 5
      patchnew%srcL1beg  = 1
      patchnew%srcL1end  = gridnew%ipc
      patchnew%srcL2beg  = 1
      patchnew%srcL2end  = gridnew%jpc
    ENDIF
  ENDDO  ! iReg

! faces 1, 3, 5

  DO ireg=2,globalnew%nRegions
    regionsnew(ireg)%nPatches = regionsnew(ireg)%nPatches + 1
    ipatch             = regionsnew(ireg)%nPatches
    patchnew           => regionsnew(ireg)%levels(1)%patches(ipatch)
    gridnew            => regionsnew(ireg)%levels(1)%grid
    patchnew%bcType    = bc_regionconf
    patchnew%bcCoupled = 0
    patchnew%align     = .true.
    patchnew%srcRegion = ireg - 1

    IF (splitdirection == 1) THEN
      patchnew%lbound    = 1
      patchnew%l1beg     = 1
      patchnew%l1end     = gridnew%jpc
      patchnew%l2beg     = 1
      patchnew%l2end     = gridnew%kpc
      patchnew%srcLbound = 2
      patchnew%srcL1beg  = 1
      patchnew%srcL1end  = gridnew%jpc
      patchnew%srcL2beg  = 1
      patchnew%srcL2end  = gridnew%kpc
    ELSE IF (splitdirection == 2) THEN
      patchnew%lbound    = 3
      patchnew%l1beg     = 1
      patchnew%l1end     = gridnew%kpc
      patchnew%l2beg     = 1
      patchnew%l2end     = gridnew%ipc
      patchnew%srcLbound = 4
      patchnew%srcL1beg  = 1
      patchnew%srcL1end  = gridnew%kpc
      patchnew%srcL2beg  = 1
      patchnew%srcL2end  = gridnew%ipc
    ELSE
      patchnew%lbound    = 5
      patchnew%l1beg     = 1
      patchnew%l1end     = gridnew%ipc
      patchnew%l2beg     = 1
      patchnew%l2end     = gridnew%jpc
      patchnew%srcLbound = 6
      patchnew%srcL1beg  = 1
      patchnew%srcL1end  = gridnew%ipc
      patchnew%srcL2beg  = 1
      patchnew%srcL2end  = gridnew%jpc
    ENDIF
  ENDDO  ! iReg

! finalize --------------------------------------------------------------------

  CALL deregisterfunction( globalnew )

END SUBROUTINE splittopology

!******************************************************************************
!
! RCS Revision history:
!
! $Log: SPLT_SplitTopology.F90,v $
! Revision 1.4  2008/12/06 08:44:51  mtcampbe
! Updated license.
!
! Revision 1.3  2008/11/19 22:18:01  mtcampbe
! Added Illinois Open Source License/Copyright
!
! Revision 1.2  2004/12/03 03:33:38  wasistho
! rflo_modinterfacessplit to splt_modinterfaces
!
! Revision 1.1  2004/12/03 02:41:15  wasistho
! added prefix
!
! Revision 1.1  2004/12/03 00:45:30  wasistho
! lower to upper case
!
! Revision 1.6  2003/05/15 02:57:07  jblazek
! Inlined index function.
!
! Revision 1.5  2003/03/20 22:32:37  haselbac
! Renamed ModInterfaces
!
! Revision 1.4  2003/03/20 19:45:54  haselbac
! Corrected mistake in phased check-in
!
! Revision 1.3  2002/09/27 00:57:11  jblazek
! Changed makefiles - no makelinks needed.
!
! Revision 1.2  2002/09/05 17:40:22  jblazek
! Variable global moved into regions().
!
! Revision 1.1  2002/07/12 21:50:08  jblazek
! Added tool to split single grid into multiple regions.
!
!******************************************************************************