This subroutine loads the information needed to read data from a file.
Example: Read time periodic velocity signals from file 'data.dat' as inlet condition without interpolation between the velocity signals.
boundary_condition =
{
{ label = 'inlet',
kind = 'inlet_ubb',
velocityX = { kind = 'combined',
temporal = { predefined = 'datafile',
filename = 'data.dat',
intp = 'none',
periodic = true
}
},
velocityY = 0.0
velocityZ = 0.0
}
}
The options for interpolating the values are \li \a 'none': No interpolation is used. The signals will be chosen by: [t1,t2[ = v1, [t2,t3[ = v2, ... \li \a 'linear': A linear interpolation is done between neighboring values according to the following equation: \f[ \frac{v_{i+1}-v_{i}}{t_{i+1}-t_{i}} \cdot (t_{sim} \% t_{period} - t_{i}) + v_{i} \f]
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(tem_from_file_temporal_type), | intent(inout) | :: | me |
temporal predefined from_file type |
||
| type(flu_State) | :: | conf |
lua state type |
|||
| integer, | intent(in) | :: | thandle |
aotus parent handle |
subroutine load_temporal_from_file( me, conf, thandle )
! ---------------------------------------------------------------------------
!> temporal predefined from_file type
type(tem_from_file_temporal_type),intent(inout) :: me
!> lua state type
type(flu_State) :: conf
!> aotus parent handle
integer, intent(in) :: thandle
! ---------------------------------------------------------------------------
! local variables
integer :: iError
! aotus handle for ramping table
integer :: rampHandle
! ---------------------------------------------------------------------------
! read the filename of the datafile
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'filename', &
& val = me%datafile, &
& ErrCode = iError )
if (btest(iError, aoterr_Fatal)) then
write(logUnit(1),*)'FATAL Error occured, while retrieving datafile:'
if (btest(iError, aoterr_NonExistent)) &
& write(logUnit(1),*)'Variable not existent!'
if (btest(iError, aoterr_WrongType)) &
& write(logUnit(1),*)'Variable has wrong type!'
write(logUnit(1),*)'STOPPING'
call tem_abort()
end if
! read the interpolation kind for the data
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'intp', &
& val = me%intp, &
& ErrCode = iError, &
& default = 'none' )
if (btest(iError, aoterr_Fatal)) then
write(logUnit(1),*)'FATAL Error occured, while retrieving interpolation:'
if (btest(iError, aoterr_NonExistent)) &
& write(logUnit(1),*)'Variable not existent!'
if (btest(iError, aoterr_WrongType)) &
& write(logUnit(1),*)'Variable has wrong type!'
write(logUnit(1),*)'STOPPING'
call tem_abort()
end if
! open the ramping table
call aot_table_open( L = conf, &
& thandle = rampHandle, &
& parent = thandle, &
& key = 'ramping')
if (rampHandle /= 0) then
me%ramp = .true.
! read the ramping time
call aot_get_val( L = conf, &
& thandle = rampHandle, &
& key = 'rampVal', &
& val = me%rampVal, &
& ErrCode = iError)
if (btest(iError, aoterr_Fatal)) then
write(logUnit(1),*)'FATAL Error occured, while retrieving ramping '// &
& 'value:'
if (btest(iError, aoterr_NonExistent)) &
& write(logUnit(1),*)'Variable not existent!'
if (btest(iError, aoterr_WrongType)) &
& write(logUnit(1),*)'Variable has wrong type!'
write(logUnit(1),*)'STOPPING'
call tem_abort()
end if
! read the ramping time
call aot_get_val( L = conf, &
& thandle = rampHandle, &
& key = 'rampT', &
& val = me%rampT, &
& ErrCode = iError)
if (btest(iError, aoterr_Fatal)) then
write(logUnit(1),*)'FATAL Error occured, while retrieving ramping time:'
if (btest(iError, aoterr_NonExistent)) &
& write(logUnit(1),*)'Variable not existent!'
if (btest(iError, aoterr_WrongType)) &
& write(logUnit(1),*)'Variable has wrong type!'
write(logUnit(1),*)'STOPPING'
call tem_abort()
end if
end if
! read the factor
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'fac', &
& val = me%fac, &
& ErrCode = iError, &
& default = 1.0_rk )
if (btest(iError, aoterr_Fatal)) then
write(logUnit(1),*)'FATAL Error occured, while retrieving factor: '
if (btest(iError, aoterr_NonExistent)) &
& write(logUnit(1),*)'Variable not existent!'
if (btest(iError, aoterr_WrongType)) &
& write(logUnit(1),*)'Variable has wrong type!'
write(logUnit(1),*)'STOPPING'
call tem_abort()
end if
! read in wether the data should be treated periodic
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'periodic', &
& val = me%periodic, &
& ErrCode = iError, &
& default = .false. )
if (btest(iError, aoterr_Fatal)) then
write(logUnit(1),*)'FATAL Error occured, while retrieving periodic:'
if (btest(iError, aoterr_NonExistent)) &
& write(logUnit(1),*)'Variable not existent!'
if (btest(iError, aoterr_WrongType)) &
& write(logUnit(1),*)'Variable has wrong type!'
write(logUnit(1),*)'STOPPING'
call tem_abort()
end if
write(logUnit(3),*)"Using data from file: "//trim(me%datafile)
write(logUnit(3),*)" interpolation: "//trim(me%intp)
if (rampHandle /= 0) then
write(logUnit(3),*)" ramping time: ",me%rampT
write(logUnit(3),*)" ramping value: ",me%rampVal
end if
write(logUnit(3),*)" periodic: ",me%periodic
! load the data from file
call load_datafile( me )
end subroutine load_temporal_from_file