JZ: Currently, we read the number of necessary coefficients from the lua script. However, there are also some automatic ways to estimate the number of coefficients necessary to obtain a certain numerical accuaracy, e.g. in [Wiscombe1980] .
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(flu_State) | :: | conf |
lua state type |
|||
integer, | intent(in) | :: | thandle |
aotus parent handle |
The number of coefficients necessary for the Fourier series
function tem_get_ncoeffs_miescat(conf, thandle ) result(nCoeffs)
! --------------------------------------------------------------------------
!> lua state type
type(flu_State) :: conf
!> aotus parent handle
integer, intent(in) :: thandle
!> The number of coefficients necessary for the Fourier series
integer :: nCoeffs
! --------------------------------------------------------------------------
integer :: iError
! --------------------------------------------------------------------------
!> @todo JZ: Currently, we read the number of necessary coefficients
!! from the lua script. However, there are also some automatic ways
!! to estimate the number of coefficients necessary to obtain a certain
!! numerical accuaracy, e.g. in [Wiscombe1980] .
! Read the frequency
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'nCoeffs', &
& val = nCoeffs, &
& ErrCode = iError )
if ( iError.ne.0 ) then
call tem_abort( 'ERROR in tem_get_ncoeffs_miescat: Not able to read' &
& // ' nCoeffs for Mie-Scatter solution from lua' )
end if
end function tem_get_ncoeffs_miescat