SetUpCaseIB Subroutine

subroutine SetUpCaseIB()

Uses

  • proc~~setupcaseib~4~~UsesGraph proc~setupcaseib~4 SetUpCaseIB module~immersed_boundaries immersed_boundaries proc~setupcaseib~4->module~immersed_boundaries module~immersed_boundaries_markers immersed_boundaries_markers module~immersed_boundaries->module~immersed_boundaries_markers module~immersed_boundaries_solids immersed_boundaries_solids module~immersed_boundaries->module~immersed_boundaries_solids iso_fortran_env iso_fortran_env module~immersed_boundaries_markers->iso_fortran_env module~leapbc leapBC module~immersed_boundaries_markers->module~leapbc module~leapblock leapBlock module~immersed_boundaries_markers->module~leapblock module~leapdiffop leapDiffOp module~immersed_boundaries_markers->module~leapdiffop module~leapeulerian leapEulerian module~immersed_boundaries_markers->module~leapeulerian module~leaphypre leapHypre module~immersed_boundaries_markers->module~leaphypre module~leapio leapIO module~immersed_boundaries_markers->module~leapio module~leapkinds leapKinds module~immersed_boundaries_markers->module~leapkinds module~leaplagrangian leapLagrangian module~immersed_boundaries_markers->module~leaplagrangian module~leapmonitor leapMonitor module~immersed_boundaries_markers->module~leapmonitor module~leapparallel leapParallel module~immersed_boundaries_markers->module~leapparallel module~leapparser leapParser module~immersed_boundaries_markers->module~leapparser module~leaptimer leapTimer module~immersed_boundaries_markers->module~leaptimer module~immersed_boundaries_solids->module~immersed_boundaries_markers module~immersed_boundaries_solids->module~leapblock module~immersed_boundaries_solids->module~leapeulerian module~immersed_boundaries_solids->module~leapio module~immersed_boundaries_solids->module~leapkinds module~immersed_boundaries_solids->module~leaplagrangian module~immersed_boundaries_solids->module~leapparallel module~leapbc->iso_fortran_env module~leapbc->module~leapblock module~leapbc->module~leapeulerian module~leapbc->module~leapio module~leapbc->module~leapkinds 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~leapio_hdf5 leapIO_hdf5 module~leapblock->module~leapio_hdf5 module~leapblock->mpi_f08 module~leapdiffop->module~leapbc module~leapdiffop->module~leapblock module~leapdiffop->module~leapeulerian module~leapdiffop->module~leapkinds module~leapdiffop->module~leapparallel 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~leaphypre->module~leapblock module~leaphypre->module~leapeulerian module~leaphypre->module~leapkinds module~leaphypre->module~leapparallel iso_c_binding iso_c_binding module~leaphypre->iso_c_binding module~leaphypre->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~leaplagrangian->module~leapblock module~leaplagrangian->module~leapio module~leaplagrangian->module~leapkinds module~leaplagrangian->module~leapparallel module~leapfilters leapfilters module~leaplagrangian->module~leapfilters module~leaplagrangian->module~leaputils module~leaplagrangian->mpi_f08 module~leapmonitor->iso_fortran_env module~leapmonitor->module~leapkinds module~leapmonitor->module~leapparallel module~leapmonitor->module~leaputils module~leapparallel->iso_fortran_env module~leapparallel->module~leapkinds module~leapparallel->mpi_f08 module~leapparser->iso_fortran_env module~leapparser->module~leapkinds module~leapcli leapCli module~leapparser->module~leapcli module~leaptimer->module~leapkinds module~leaptimer->module~leapparallel module~leaptimer->module~leapparser module~leaptimer->module~leaputils module~leapcli->module~leapkinds module~leapfilters->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 module~leaputils->module~leapkinds

Builds and writes immersed boundary.

Arguments

None

Calls

proc~~setupcaseib~4~~CallsGraph proc~setupcaseib~4 SetUpCaseIB markers markers proc~setupcaseib~4->markers n n proc~setupcaseib~4->n none~finalize~27 marker_set%Finalize proc~setupcaseib~4->none~finalize~27 none~get parser_obj%Get proc~setupcaseib~4->none~get none~initialize~27 marker_set%Initialize proc~setupcaseib~4->none~initialize~27 p p proc~setupcaseib~4->p proc~lagrangian_set_applyperiodicity lagrangian_set%lagrangian_set_ApplyPeriodicity proc~setupcaseib~4->proc~lagrangian_set_applyperiodicity proc~lagrangian_set_communicate lagrangian_set%lagrangian_set_Communicate proc~setupcaseib~4->proc~lagrangian_set_communicate proc~lagrangian_set_localize lagrangian_set%lagrangian_set_Localize proc~setupcaseib~4->proc~lagrangian_set_localize proc~lagrangian_set_resize lagrangian_set%lagrangian_set_Resize proc~setupcaseib~4->proc~lagrangian_set_resize proc~lagrangian_set_setwritefilename lagrangian_set%lagrangian_set_SetWriteFileName proc~setupcaseib~4->proc~lagrangian_set_setwritefilename proc~parallel_obj_rankisroot parallel_obj%parallel_obj_RankIsRoot proc~setupcaseib~4->proc~parallel_obj_rankisroot proc~respart_set_final ResPart_set%ResPart_set_Final none~finalize~27->proc~respart_set_final 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~respart_set_init ResPart_set%ResPart_set_Init none~initialize~27->proc~respart_set_init proc~lagrangian_set_communicate->proc~lagrangian_set_resize mpi_gather mpi_gather proc~lagrangian_set_communicate->mpi_gather mpi_recv mpi_recv proc~lagrangian_set_communicate->mpi_recv mpi_send mpi_send proc~lagrangian_set_communicate->mpi_send proc~lagrangian_set_recycle lagrangian_set%lagrangian_set_Recycle proc~lagrangian_set_communicate->proc~lagrangian_set_recycle proc~block_obj_getowningcell block_obj%block_obj_GetOwningCell proc~lagrangian_set_localize->proc~block_obj_getowningcell proc~lagrangian_set_recycle->proc~lagrangian_set_resize 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~respart_set_final->none~finalize~27 proc~lagrangian_set_freempitype lagrangian_set%lagrangian_set_FreeMPIType proc~respart_set_final->proc~lagrangian_set_freempitype proc~respart_set_init->none~initialize~27 proc~respart_set_init->proc~lagrangian_set_resize proc~lagrangian_set_creatempitype lagrangian_set%lagrangian_set_CreateMPIType proc~respart_set_init->proc~lagrangian_set_creatempitype proc~respart_set_setobjecttype ResPart_set%ResPart_set_SetObjectType proc~respart_set_init->proc~respart_set_setobjecttype 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 SetMPIDataTypeParams SetMPIDataTypeParams proc~lagrangian_set_creatempitype->SetMPIDataTypeParams mpi_type_commit mpi_type_commit proc~lagrangian_set_creatempitype->mpi_type_commit mpi_type_create_resized mpi_type_create_resized proc~lagrangian_set_creatempitype->mpi_type_create_resized mpi_type_create_struct mpi_type_create_struct proc~lagrangian_set_creatempitype->mpi_type_create_struct mpi_type_get_extent mpi_type_get_extent proc~lagrangian_set_creatempitype->mpi_type_get_extent mpi_type_size mpi_type_size proc~lagrangian_set_creatempitype->mpi_type_size mpi_type_free mpi_type_free proc~lagrangian_set_freempitype->mpi_type_free

Called by

proc~~setupcaseib~4~~CalledByGraph proc~setupcaseib~4 SetUpCaseIB program~main~14 main program~main~14->proc~setupcaseib~4

Source Code

    subroutine SetUpCaseIB
      !> Builds and writes immersed boundary.
      use immersed_boundaries
      implicit none
      type(marker_set)    :: IB
      character(str64)    :: filename
      real(wp)            :: diam
      real(wp)            :: xc(3)
      real(wp)            :: d_axis
      integer             :: n_axis
      real(wp)            :: d_theta
      integer             :: n_theta
      integer             :: i,j,m
      real(wp), parameter :: Pi=4.0*atan(1.0_wp)

      ! Get info from parser
      call parser%Get("IB IC file",        filename)
      call parser%Get("Cylinder diameter", diam    )
      call parser%Get("Cylinder center",   xc      )

      ! Initialize Immersed Boundaries
      call IB%Initialize('IB',block,parallel)

      ! Create a cylinder
      if (parallel%RankIsRoot()) then

        ! Count the number of markers
        d_axis = 0.5_wp*block%dx(3)
        n_axis = max(1,floor((block%pmax(3)-block%pmin(3))/d_axis))
        d_axis = (block%pmax(3)-block%pmin(3))/real(n_axis,wp)

        n_theta=ceiling(Pi*diam/d_axis)
        d_theta=2.0_wp*Pi/real(n_theta,wp)

        ! Create markers
        call IB%Resize(n_theta*n_axis)

        select type (markers=>IB%p)
        type is (marker_obj)
          do i=1,n_theta
            do j=1,n_axis
              ! ID
              m = n_axis*(i-1)+j
              markers(m)%id  = int(m,kind=8)
              ! Position
              markers(m)%p(1)= xc(1)+0.5_wp*diam*cos((i-1)*d_theta + d_theta/2.0_wp)
              markers(m)%p(2)= xc(2)+0.5_wp*diam*sin((i-1)*d_theta + d_theta/2.0_wp)
              markers(m)%p(3)= xc(3)+(j-1)*d_axis + d_axis/2.0_wp + block%pmin(3)
              ! Surface Area
              markers(m)%SA  = 0.5_wp*diam*d_theta*d_axis
              ! Normal
              markers(m)%n(1)= cos((i-1)*d_theta + d_theta/2.0_wp)
              markers(m)%n(2)= sin((i-1)*d_theta + d_theta/2.0_wp)
              markers(m)%n(3)= 0.0_wp
              ! velocity
              markers(m)%v   = 0.0_wp
            end do
          end do
        end select
      end if

      ! Treatment for periodicity
      call IB%ApplyPeriodicity()

      ! Send to the right rank
      call IB%Communicate()

      ! Localize marker points
      call IB%Localize()

      ! Write data to disk
      call IB%SetWriteFileName(filename)
      call IB%Write(0,0.0_WP)

      ! Finalize
      call IB%Finalize()
      return
    end subroutine SetUpCaseIB