tem_load_spongeLayer_radial Subroutine

public subroutine tem_load_spongeLayer_radial(me, conf, thandle, nDim, nComp, stateName)

This subroutine load data for radial sponge layer Example:

 spatial = {
   --supported options: 'spongelayer_radial','sponge_radial_2d',
   --                   'viscous_spongelayer_radial',
   --                   'viscous_spongelayer_radial_2d'
   predefined = 'viscous_spongelayer_radial', 
   origin = {0.0,0.0,0.0},
   radius = 1.0, -- Sponge start
   thickness = 0.3,
   damp_profile = 'linear', --'exponential', 'polynomial_n5', 'polynomial_n6'
   damp_factor = 0.5,
   damp_exponent = 1.0,
   target_state = {
     Default: density, velocityX, velocityY, velocityZ and pressure
     viscosity = 1e-3
 }

Arguments

Type IntentOptional Attributes Name
type(tem_spongeLayer_radial_type), intent(out) :: me

Radial spongeLayer data type
type(flu_State) :: conf

lua state type
integer, intent(in) :: thandle

aotus parent handle
integer, intent(in) :: nDim

number of Dimension for nonViscous sponges
integer, intent(in) :: nComp

Number of component of St-Fun variable under which this spatial function is defined
character(len=*), intent(in), optional :: stateName

Load stateName from target_state table

Calls

proc~~tem_load_spongelayer_radial~~CallsGraph proc~tem_load_spongelayer_radial tem_load_spongeLayer_radial aot_get_val aot_get_val proc~tem_load_spongelayer_radial->aot_get_val proc~load_spongelayer load_spongeLayer proc~tem_load_spongelayer_radial->proc~load_spongelayer proc~tem_abort tem_abort proc~tem_load_spongelayer_radial->proc~tem_abort proc~load_spongelayer->aot_get_val proc~load_spongelayer->proc~tem_abort aot_table_close aot_table_close proc~load_spongelayer->aot_table_close aot_table_open aot_table_open proc~load_spongelayer->aot_table_open proc~load_defaulttargetstate load_defaultTargetState proc~load_spongelayer->proc~load_defaulttargetstate mpi_abort mpi_abort proc~tem_abort->mpi_abort proc~load_defaulttargetstate->aot_get_val proc~load_defaulttargetstate->proc~tem_abort proc~load_defaulttargetstate->aot_table_close proc~load_defaulttargetstate->aot_table_open

Called by

proc~~tem_load_spongelayer_radial~~CalledByGraph proc~tem_load_spongelayer_radial tem_load_spongeLayer_radial proc~load_spatial_predefined load_spatial_predefined proc~load_spatial_predefined->proc~tem_load_spongelayer_radial proc~tem_load_spatial tem_load_spatial proc~tem_load_spatial->proc~load_spatial_predefined proc~load_spacetime_predefined load_spacetime_predefined proc~load_spacetime_predefined->proc~tem_load_spatial proc~tem_load_ic tem_load_ic proc~tem_load_ic->proc~tem_load_spatial proc~tem_load_spacetime_single tem_load_spacetime_single proc~tem_load_spacetime_single->proc~load_spacetime_predefined

Source Code

  subroutine tem_load_spongeLayer_radial(me, conf, thandle, nDim, nComp, &
    &                                    stateName)
    ! --------------------------------------------------------------------------
    !> Radial spongeLayer data type
    type(tem_spongeLayer_radial_type), intent(out) :: me
    !> lua state type
    type(flu_State) :: conf
    !> aotus parent handle
    integer, intent(in) :: thandle
    !> number of Dimension for nonViscous sponges
    integer, intent(in) :: nDim
    !> Number of component of St-Fun variable under which this spatial function
    !! is defined
    integer, intent(in) :: nComp
    !> Load stateName from target_state table
    character(len=*), intent(in), optional :: stateName
    ! --------------------------------------------------------------------------
    integer :: iError
    integer :: vError(3), errfatal(3)
    ! --------------------------------------------------------------------------

    errfatal = aotErr_Fatal
    ! Sponge origin
    call aot_get_val( L       = conf,      &
      &               thandle = thandle,   &
      &               key     = 'origin',  &
      &               val     = me%origin, &
      &               ErrCode = vError     )
    if (any(btest( vError, errFatal )) ) then
      write(logUnit(1),*) 'ERROR reading the sponge origin, ' &
        &              // 'origin is not well defined. '   &
        &              // 'It should have 3 entries for each coordinate.'
      call tem_abort()
    end if
    write(logUnit(1),*) ' * Origin =', me%origin

    ! Sponge inner radius
    call aot_get_val( L       = conf,      &
      &               thandle = thandle,   &
      &               key     = 'radius',  &
      &               val     = me%radius, &
      &               ErrCode = iError     )
    if (btest(iError, aotErr_Fatal)) then
      write(logUnit(1),*) 'ERROR reading the sponge inner radius. '
      call tem_abort()
    end if
    write(logUnit(1),*) ' * Inner radius =', me%radius

    ! Load base information required for sponge layer definition like
    ! damp_factor, damp_exponent and target_state
    call load_spongeLayer( conf      = conf,                     &
      &                    thandle   = thandle,                  &
      &                    me        = me%spongeLayer_base_type, &
      &                    nDim      = nDim,                     &
      &                    nComp     = nComp,                    &
      &                    stateName = stateName                 )

  end subroutine tem_load_spongeLayer_radial