SetUpCaseBCS Subroutine

subroutine SetUpCaseBCS()

Uses

  • proc~~setupcasebcs~13~~UsesGraph proc~setupcasebcs~13 SetUpCaseBCS module~leapbc leapBC proc~setupcasebcs~13->module~leapbc iso_fortran_env iso_fortran_env module~leapbc->iso_fortran_env module~leapblock leapBlock module~leapbc->module~leapblock module~leapeulerian leapEulerian module~leapbc->module~leapeulerian module~leapio leapIO module~leapbc->module~leapio module~leapkinds leapKinds module~leapbc->module~leapkinds module~leapparallel leapParallel module~leapbc->module~leapparallel module~leaputils leapUtils module~leapbc->module~leaputils mpi_f08 mpi_f08 module~leapbc->mpi_f08 module~leapblock->iso_fortran_env module~leapblock->module~leapkinds module~leapblock->module~leapparallel module~leapblock->mpi_f08 module~leapio_hdf5 leapIO_hdf5 module~leapblock->module~leapio_hdf5 module~leapeulerian->iso_fortran_env module~leapeulerian->module~leapblock module~leapeulerian->module~leapio module~leapeulerian->module~leapkinds module~leapeulerian->module~leapparallel module~leapeulerian->module~leaputils module~leapeulerian->mpi_f08 module~leapio_h5hut leapIO_h5hut module~leapio->module~leapio_h5hut module~leapio->module~leapio_hdf5 module~leapio_silo leapIO_silo module~leapio->module~leapio_silo module~leapio_xdmf leapIO_xdmf module~leapio->module~leapio_xdmf module~leapkinds->iso_fortran_env module~leapparallel->iso_fortran_env module~leapparallel->module~leapkinds module~leapparallel->mpi_f08 module~leaputils->module~leapkinds module~leapio_h5hut->module~leapkinds module~leapio_h5hut->module~leapparallel module~leapio_h5hut->module~leapio_hdf5 module~leapio_hdf5->module~leapkinds module~leapio_hdf5->module~leapparallel module~leapio_hdf5->module~leaputils hdf5 hdf5 module~leapio_hdf5->hdf5 module~leapio_silo->module~leapkinds module~leapio_silo->module~leapparallel module~leapio_silo->module~leaputils module~leapio_silo->mpi_f08 module~leapio_xdmf->module~leapkinds module~leapio_xdmf->module~leaputils

Defines boundary conditions.

Arguments

None

Calls

proc~~setupcasebcs~13~~CallsGraph proc~setupcasebcs~13 SetUpCaseBCS none~get parser_obj%Get proc~setupcasebcs~13->none~get proc~bc_set_add bc_set%bc_set_Add proc~setupcasebcs~13->proc~bc_set_add proc~bc_set_final bc_set%bc_set_Final proc~setupcasebcs~13->proc~bc_set_final proc~bc_set_getbcpointer bc_set%bc_set_GetBCPointer proc~setupcasebcs~13->proc~bc_set_getbcpointer proc~bc_set_getextents bc_set%bc_set_GetExtents proc~setupcasebcs~13->proc~bc_set_getextents proc~bc_set_init bc_set%bc_set_Init proc~setupcasebcs~13->proc~bc_set_init proc~bc_set_setbc bc_set%bc_set_SetBC proc~setupcasebcs~13->proc~bc_set_setbc proc~parser_obj_read0d parser_obj%parser_obj_read0D none~get->proc~parser_obj_read0d proc~parser_obj_read1d parser_obj%parser_obj_read1D none~get->proc~parser_obj_read1d proc~bc_set_checkbounds bc_set%bc_set_CheckBounds proc~bc_set_add->proc~bc_set_checkbounds proc~bc_set_expand bc_set%bc_set_Expand proc~bc_set_add->proc~bc_set_expand proc~bc_set_getsidedirbynormal bc_set%bc_set_GetSideDirByNormal proc~bc_set_add->proc~bc_set_getsidedirbynormal proc~bc_set_updateextents bc_set%bc_set_UpdateExtents proc~bc_set_add->proc~bc_set_updateextents proc~hashtbl_obj_hashstring hashtbl_obj%hashtbl_obj_HashString proc~bc_set_add->proc~hashtbl_obj_hashstring proc~hashtbl_obj_put hashtbl_obj%hashtbl_obj_Put proc~bc_set_add->proc~hashtbl_obj_put proc~region_obj_init region_obj%region_obj_Init proc~bc_set_add->proc~region_obj_init proc~region_obj_final region_obj%region_obj_Final proc~bc_set_final->proc~region_obj_final proc~bc_set_getregionindex bc_set%bc_set_GetRegionIndex proc~bc_set_getbcpointer->proc~bc_set_getregionindex proc~region_obj_getbcindex region_obj%region_obj_GetBCIndex proc~bc_set_getbcpointer->proc~region_obj_getbcindex proc~bc_set_getextents->proc~bc_set_getregionindex proc~hashtbl_obj_init hashtbl_obj%hashtbl_obj_Init proc~bc_set_init->proc~hashtbl_obj_init proc~bc_set_setbc->proc~bc_set_getextents proc~bc_set_setbc->proc~bc_set_getregionindex proc~region_obj_add region_obj%region_obj_Add proc~bc_set_setbc->proc~region_obj_add proc~bc_set_getregionindex->proc~hashtbl_obj_hashstring none~get~4 hashtbl_obj%Get proc~bc_set_getregionindex->none~get~4 proc~bc_set_updateextents->proc~bc_set_getregionindex proc~sllist_obj_put sllist_obj%sllist_obj_Put proc~hashtbl_obj_put->proc~sllist_obj_put none~assigndefault parser_obj%AssignDefault proc~parser_obj_read0d->none~assigndefault proc~parser_obj_fetchlabelid parser_obj%parser_obj_FetchLabelID proc~parser_obj_read0d->proc~parser_obj_fetchlabelid proc~parser_obj_read1d->none~assigndefault proc~parser_obj_read1d->proc~parser_obj_fetchlabelid proc~region_obj_add->proc~hashtbl_obj_hashstring proc~region_obj_add->proc~hashtbl_obj_put proc~region_obj_expand region_obj%region_obj_Expand proc~region_obj_add->proc~region_obj_expand proc~hashtbl_obj_final hashtbl_obj%hashtbl_obj_Final proc~region_obj_final->proc~hashtbl_obj_final proc~region_obj_getbcindex->proc~hashtbl_obj_hashstring proc~region_obj_getbcindex->none~get~4 proc~region_obj_init->proc~hashtbl_obj_init proc~parser_obj_assigndefault0d parser_obj%parser_obj_AssignDefault0D none~assigndefault->proc~parser_obj_assigndefault0d proc~parser_obj_assigndefault1d parser_obj%parser_obj_AssignDefault1D none~assigndefault->proc~parser_obj_assigndefault1d proc~hashtbl_obj_get_int4 hashtbl_obj%hashtbl_obj_Get_int4 none~get~4->proc~hashtbl_obj_get_int4 proc~hashtbl_obj_get_int8 hashtbl_obj%hashtbl_obj_Get_int8 none~get~4->proc~hashtbl_obj_get_int8 proc~hashtbl_obj_get_real_dp hashtbl_obj%hashtbl_obj_Get_real_dp none~get~4->proc~hashtbl_obj_get_real_dp proc~hashtbl_obj_get_real_sp hashtbl_obj%hashtbl_obj_Get_real_sp none~get~4->proc~hashtbl_obj_get_real_sp proc~sllist_obj_put->proc~sllist_obj_put none~get~3 sllist_obj%Get proc~hashtbl_obj_get_int4->none~get~3 proc~hashtbl_obj_get_int8->none~get~3 proc~hashtbl_obj_get_real_dp->none~get~3 proc~hashtbl_obj_get_real_sp->none~get~3 proc~sllist_obj_get_int4 sllist_obj%sllist_obj_Get_int4 none~get~3->proc~sllist_obj_get_int4 proc~sllist_obj_get_int8 sllist_obj%sllist_obj_Get_int8 none~get~3->proc~sllist_obj_get_int8 proc~sllist_obj_get_real_dp sllist_obj%sllist_obj_Get_real_dp none~get~3->proc~sllist_obj_get_real_dp proc~sllist_obj_get_real_sp sllist_obj%sllist_obj_Get_real_sp none~get~3->proc~sllist_obj_get_real_sp proc~sllist_obj_get_int4->proc~sllist_obj_get_int4 proc~sllist_obj_get_int8->proc~sllist_obj_get_int8 proc~sllist_obj_get_real_dp->proc~sllist_obj_get_real_dp proc~sllist_obj_get_real_sp->proc~sllist_obj_get_real_sp

Called by

proc~~setupcasebcs~13~~CalledByGraph proc~setupcasebcs~13 SetUpCaseBCS program~main~13 main program~main~13->proc~setupcasebcs~13

Source Code

    subroutine SetUpCaseBCS()
      !> Defines boundary conditions.
      use leapBC
      implicit none
      ! Work variables
      type(bc_set)      :: bcs
      type(extent_obj)  :: extents
      integer           :: i,j,k
      real(wp)          :: xhi(3)
      real(wp)          :: xlo(3)
      real(wp)          :: L(3)
      real(wp)          :: Vlid
      real(wp), pointer :: V1bc(:,:,:)
      real(wp), pointer :: V2bc(:,:,:)
      real(wp), pointer :: V3bc(:,:,:)

      ! Get lid velocity and domain size
      call parser%Get("Lid velocity", Vlid, default = 1.0_wp  )
      call parser%Get("Domain size", L )

      ! Initialize utility that handles boundary conditions
      call bcs%Initialize(block,parallel)

      ! Setup regions where boundary conditions apply
      ! - Left boundary
      xlo = [-0.5_wp*L(1),-0.5_wp*L(2),-0.5_wp*L(3)]
      xhi = [-0.5_wp*L(1), 0.5_wp*L(2), 0.5_wp*L(3)]
      call bcs%Add('Left', xlo, xhi, normal = '+x1')

      ! - Right boundary
      xlo = [ 0.5_wp*L(1),-0.5_wp*L(2),-0.5_wp*L(3)]
      xhi = [ 0.5_wp*L(1), 0.5_wp*L(2), 0.5_wp*L(3)]
      call bcs%Add('Right', xlo, xhi, normal = '-x1')

      ! - Bottom boundary
      xlo = [-0.5_wp*L(1),-0.5_wp*L(2),-0.5_wp*L(3)]
      xhi = [ 0.5_wp*L(1),-0.5_wp*L(2), 0.5_wp*L(3)]
      call bcs%Add('Bottom', xlo, xhi, normal = '+x2')

      ! - Top boundary
      xlo = [-0.5_wp*L(1), 0.5_wp*L(2),-0.5_wp*L(3)]
      xhi = [ 0.5_wp*L(1), 0.5_wp*L(2), 0.5_wp*L(3)]
      call bcs%Add('Top', xlo, xhi, normal = '-x2')

      call bcs%SetBC('Top',    BC_INFLOW)
      call bcs%SetBC('Bottom', BC_WALL  )
      call bcs%SetBC('Left',   BC_WALL  )
      call bcs%SetBC('Right',  BC_WALL  )

      ! Set the inflow values
      extents = bcs%GetExtents('Top')
      call bcs%GetBCPointer('Top','V1',V1bc)
      call bcs%GetBCPointer('Top','V2',V2bc)
      call bcs%GetBCPointer('Top','V3',V3bc)
      do k=extents%lo(3),extents%hi(3)
        do j=extents%lo(2),extents%hi(2)
          do i=extents%lo(1),extents%hi(1)
            V1bc(i,j,k) = Vlid
            V2bc(i,j,k) = 0.0_wp
            V3bc(i,j,k) = 0.0_wp
          end do
        end do
      end do

      ! Write boundary conditions
      call bcs%Write(0,0.0_wp)

      ! Clear data
      V1bc => null()
      V2bc => null()
      V3bc => null()
      call bcs%Finalize()

      return
    end subroutine SetUpCaseBCS