get the z component
read list of vertices
how many lists provided, here we just have one here we need to read out the kind and than point to the specific movement of the polygon
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(tem_polygon_material_type), | intent(out) | :: | me |
Polygon data structure to fill with information provided by the user in config. |
||
| type(flu_State) | :: | conf |
Handle to the Lua script containing the polygon definition |
|||
| integer, | intent(in), | optional | :: | thandle |
Handle for the table containing the polygon definition. |
subroutine tem_polygon_material_multi_load(me, conf, thandle)
! ----------------------------------------------------------------------
!> Polygon data structure to fill with information provided
!! by the user in config.
type(tem_polygon_material_type), intent(out) :: me
!> Handle to the Lua script containing the polygon definition
type(flu_state) :: conf
!> Handle for the table containing the polygon definition.
integer, intent(in), optional :: thandle
! ----------------------------------------------------------------------
real(kind=rk), allocatable :: defout(:)
integer :: vertex_table, valtable, vertices_table
integer :: movement_table, parameter_table
integer :: iVertex, ipoly
integer :: iError
integer, allocatable :: vError(:)
integer :: iError_v(2)
! ----------------------------------------------------------------------
write(logUnit(1),*) 'Loading predefined function for multi body polygon:'
valtable =0
!> get the z component
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'zmin', &
& val = me%zmin, &
& default = 0.0_rk, &
& ErrCode = iError )
if (btest(iError, aoterr_fatal)) then
write(logunit(1),*) 'ERROR in tem_polygon_material_load: Not able to'
write(logunit(1),*) ' get a value for zmin!'
write(logunit(1),*) ' This also applies if no inval is provided.'
call tem_abort()
end if
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'zmax', &
& val = me%zmax, &
& default = 0.0_rk, &
& ErrCode = iError )
if (btest(iError, aoterr_fatal)) then
write(logunit(1),*) 'ERROR in tem_polygon_material_load: Not able to'
write(logunit(1),*) ' get a value for zmax!'
write(logunit(1),*) ' This also applies if no inval is provided.'
call tem_abort()
end if
!> read list of vertices
call aot_table_open( L = conf, &
& parent = thandle, &
& key = 'vertices', &
& thandle = vertices_table )
if (vertices_table == 0) then
write(logunit(1),*) 'ERROR in tem_polygon_material_load: No vertices'
write(logunit(1),*) ' defined, unable to set up a polynomial!'
write(logunit(1),*) ' Please define vertices tables with lists of'
write(logunit(1),*) ' 2D points.'
call tem_abort()
end if
me%nPoly = aot_table_length( L = conf, &
& thandle = vertices_table )
!> how many lists provided, here we just have one
allocate(me%poly_list(me%npoly))
do iPoly = 1, me%npoly
call aot_table_open( L = conf, &
& parent = vertices_table, &
& pos = ipoly, &
& thandle = vertex_table )
if (vertex_table == 0) then
write(logunit(1),*) 'ERROR in tem_polygon_material_load: No vertices'
write(logunit(1),*) ' Please define vertex tables with lists of'
write(logunit(1),*) ' 2D points.'
call tem_abort()
end if
me%poly_list(ipoly)%nVertices = aot_table_length( L = conf, &
& thandle = vertex_table )
allocate(me%poly_list(ipoly)%vertex(me%poly_list(ipoly)%nVertices, 2))
do iVertex=1,me%poly_list(ipoly)%nVertices
call aot_get_val( val = me%poly_list(ipoly)%vertex(iVertex, :), &
& ErrCode = iError_v, &
& L = conf, &
& thandle = vertex_table, &
& pos = iVertex )
if (any(iError_v /= 0)) then
write(logunit(1),*) 'ERROR in tem_polygon_material_load: Not able to'
write(logunit(1),*) ' obtain vertex ', iVertex, '!'
write(logunit(1),*) ' Vertices have to be vectors of length 2,'
write(logunit(1),*) ' with real numbers as entries.'
call tem_abort()
end if
end do
call aot_table_close( L = conf, &
& thandle = vertex_table )
end do
call aot_table_close( L = conf, &
& thandle = vertices_table )
call aot_table_open( L = conf, &
& parent = thandle, &
& key = 'movement', &
& thandle = movement_table )
if (movement_table == 0) then
write(logunit(1),*) 'ERROR in tem_polygon_material_load: No movement'
write(logunit(1),*) ' defined, unable to continue!'
write(logunit(1),*) &
& ' Please define movement table with a movement kind'
call tem_abort()
end if
call aot_get_val( L = conf, &
& thandle = movement_table, &
& key = 'movement_kind', &
& val = me%moving%movement_kind, &
& default = 'No movement', &
& ErrCode = iError )
call aot_table_close( L = conf, &
& thandle = movement_table )
select case(me%moving%movement_kind)
case( 'lin_multi_body_2d', 'lin_multi_body_3d')
call aot_table_open( L = conf, &
& parent = thandle, &
& key = 'lin_parameter', &
& thandle = parameter_table )
if (parameter_table == 0) then
write(logunit(1),*) &
& 'ERROR in tem_polygon_material_load: lin_parameter'
write(logunit(1),*) ' not defined, unable to set movement!'
write(logunit(1),*) ' Define lin_parameter table for multi body'
call tem_abort()
end if
me%nComponents = aot_table_length( L = conf, &
& thandle = parameter_table )
call aot_table_close( L = conf, &
& thandle = parameter_table )
allocate(me%moving%lin_parameter(me%nComponents))
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'lin_parameter', &
& val = me%moving%lin_parameter, &
& maxlength = me%nComponents, &
& ErrCode = VError )
if (any(btest(vError, aoterr_Fatal))) then
write(logunit(1),*) &
& 'ERROR in tem_polygon_material_load: No movement'
write(logunit(1),*) &
& ' values defined, unable to set up the movement!'
write(logunit(1),*) &
& ' Please define lin_parameter!'
write(logunit(1),*) &
& ' For linear movement set velocity in X,Y for 2D!'
call tem_abort()
end if
case( 'sin_multi_body_2d','sin_multi_body_3d')
call aot_table_open( L = conf, &
& parent = thandle, &
& key = 'sin_parameter', &
& thandle = parameter_table )
if (parameter_table == 0) then
write(logunit(1),*) &
& 'ERROR in tem_polygon_material_load: sin_parameter'
write(logunit(1),*) ' not defined, unable to set movement!'
write(logunit(1),*) &
& ' Please define sin_parameter table for multi body'
call tem_abort()
end if
me%nComponents = aot_table_length( L = conf, &
& thandle = parameter_table )
call aot_table_close( L = conf, &
& thandle = parameter_table )
allocate(me%moving%sin_parameter(me%nComponents))
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'sin_parameter', &
& val = me%moving%sin_parameter, &
& maxlength = me%nComponents, &
& ErrCode = VError )
if (any(btest(vError, aoterr_Fatal))) then
write(logunit(1),*) 'ERROR in tem_polygon_material_load: No movement'
write(logunit(1),*) &
& ' values defined, unable to set up the movement!'
write(logunit(1),*) ' Please define sin_parameter!'
write(logunit(1),*) &
& ' For sinefuction set amplitude and frequency in X,'
write(logunit(1),*) ' than the information for the Y direction!'
end if
case( 'rot_multi_body_2d','rot_multi_body_3d')
call aot_table_open( L = conf, &
& parent = thandle, &
& key = 'rot_parameter', &
& thandle = parameter_table )
if (parameter_table == 0) then
write(logunit(1),*) &
& 'ERROR in tem_polygon_material_load: rot_parameter'
write(logunit(1),*) ' not defined, unable to set movement!'
write(logunit(1),*) ' Please define rot_parameter table'
call tem_abort()
end if
me%nComponents = aot_table_length( L = conf, &
& thandle = parameter_table )
call aot_table_close( L = conf, &
& thandle = parameter_table )
allocate(me%moving%rot_parameter(me%nComponents))
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'rot_parameter', &
& val = me%moving%rot_parameter, &
& maxlength = me%nComponents, &
& ErrCode = VError )
if (any(btest(vError, aoterr_Fatal))) then
write(logunit(1),*) 'ERROR in tem_polygon_material_load: No movement'
write(logunit(1),*) &
& ' values defined, unable to set up the movement!'
write(logunit(1),*) ' Please define rot_parameter!'
write(logunit(1),*) &
& ' For rotation provide the X and Y direction of the'
write(logunit(1),*) ' rotation_point and the rotation speed! '
end if
case default
write(logUnit(1),*) 'NO movement prescribed for the polygon'
end select
deallocate(vError)
write(logUnit(1),*) 'The multi geometry movement is defined by: ', &
& me%moving%movement_kind
!! here we need to read out the kind and than point to the specific
!! movement of the polygon
call aot_table_open( L = conf, &
& parent = thandle, &
& key = 'inval', &
& thandle = valtable )
if (valtable == 0) then
! inval not provided as a table, try to read it as a scalar.
allocate(me%inval(1))
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'inval', &
& val = me%inval(1), &
& default = 1.0_rk, &
& ErrCode = iError )
if (btest(iError, aoterr_Fatal)) then
write(logunit(1),*) 'ERROR in tem_polygon_material_load: Not able to'
write(logunit(1),*) ' get value for inval!'
call tem_abort()
end if
me%nComponents = 1
allocate(me%outval(me%nComponents))
! Outval needs to be consistent with the inval definition, if inval was
! defined as a scalar, outval also has to be a scalar!
! We do not check for tables with single entries in this case.
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'outval', &
& val = me%outval(1), &
& default = 0.0_rk, &
& ErrCode = iError )
if (btest(iError, aoterr_fatal)) then
write(logunit(1),*) 'ERROR in tem_polygon_material_load: Not able to'
write(logunit(1),*) ' get a value for outval!'
write(logunit(1),*) ' Note, that outval needs be a scalar, as'
write(logunit(1),*) ' inval is provided as a scalar.'
write(logunit(1),*) ' This also applies if no inval is provided.'
call tem_abort()
end if
else
! Intable is a table, close it an read it into an array.
call aot_table_close(L = conf, thandle = valtable)
! Value to use inside the polygon
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'inval', &
& val = me%inval, &
& maxlength = 20, &
& default = [1.0_rk], &
& ErrCode = vError )
if (any(btest(vError, aoterr_Fatal))) then
write(logunit(1),*) 'ERROR in tem_polygon_material_load: Not able to'
write(logunit(1),*) ' get values for inval!'
call tem_abort()
end if
me%nComponents = size(me%inval)
deallocate(vError)
! Definition of outval needs to be consistent with inval, it has to have
! the same number of components, and also needs to be a vector.
! However, we define a default of all zeroes, so if outval is 0 for all
! components, this definition can be omitted in the user definition.
allocate(me%outval(me%nComponents))
allocate(vError(me%nComponents))
allocate(defout(me%nComponents))
defout = 0.0_rk
! Value to use outside the polygon
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'outval', &
& val = me%outval, &
& default = defout, &
& ErrCode = vError )
if (any(btest(vError,aoterr_Fatal))) then
write(logunit(1),*) 'ERROR in tem_polygon_material_load: Not able to'
write(logunit(1),*) ' get a value for outval!'
write(logunit(1),*) ' Note, that outval needs to have the same'
write(logunit(1),*) ' length as inval.'
call tem_abort()
end if
deallocate(vError)
deallocate(defout)
end if
end subroutine tem_polygon_material_multi_load