load crvp variables to set initial condition
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(flu_State) | :: | conf |
lua state type |
|||
| integer, | intent(in) | :: | thandle |
aotus parent handle |
||
| type(ic_2dcrvp_type), | intent(out) | :: | me |
Global gauss pulse data type |
subroutine load_ic_2dcrvp( conf, thandle, me )
! ---------------------------------------------------------------------------
!> lua state type
type(flu_State) :: conf
!> aotus parent handle
integer, intent(in) :: thandle
!> Global gauss pulse data type
type(ic_2dcrvp_type), intent(out) :: me
! ---------------------------------------------------------------------------
integer :: cent_handle
integer :: i
integer :: iError
! ---------------------------------------------------------------------------
!center
call aot_table_open( L = conf, &
& parent = thandle, &
& thandle = cent_handle, &
& key = 'center' )
if (aot_table_length(L=conf, thandle=cent_handle) == 3) then
do i=1,3
call aot_get_val( L = conf, &
& thandle = cent_handle, &
& pos = i, &
& val = me%center(i), &
& ErrCode = iError)
end do
else
write(*,*) 'ERROR while reading the center of a vortices,'
write(*,*) 'should have 3 entries for each coordinate!'
STOP
end if
call aot_table_close(conf, cent_handle)
write(logUnit(1),*) ' * center =', me%center
! get physical time steps dx and dt
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 't', &
& val = me%t, &
& ErrCode = iError, &
& default = 0._rk )
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'kappa', &
& val = me%kappa, &
& ErrCode = iError, &
& default = 1._rk )
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'p0', &
& val = me%p0, &
& ErrCode = iError, &
& default = 0._rk )
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'rho0', &
& val = me%rho0, &
& ErrCode = iError, &
& default = 1._rk )
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'cs', &
& val = me%cs, &
& ErrCode = iError )
if(( iError .gt. 0 )) then
if( me%p0 .le. 0._rk ) then
write(logUnit(1),*)' No speed of sound is given and the reference '
write(logUnit(1),*)' pressure for CRVP is 0. Choose > 0! '
write(logUnit(1),*)' Either give the speed of sound cs = ... or use '
write(logUnit(1),*)' a reference pressure > 0'
call tem_abort
end if
me%cs = sqrt( me%kappa*me%p0/me%rho0 )
end if
if( me%cs .le. 0._rk ) then
write(logUnit(1),*)' Error: Speed of sound <= 0 '
call tem_abort
end if
! Gaussian vortex core model for pressure
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'gpmodel', &
& val = me%pressGaussModel, &
& ErrCode = iError, &
& default = .true. )
me%matchFactor = 2.2_rk
if( me%pressGaussModel ) then
write(logUnit(1),*) ' * gaussian model for pressure in rC'
else
write(logUnit(1),*) ' * no model for pressure in rC'
end if
! Rankine vortex core model for velocity
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'rankine', &
& val = me%rankineModel, &
& ErrCode = iError, &
& default = .true. )
if( me%rankineModel ) then
write(logUnit(1),*) ' * Rankine vortex core model'
end if
!halfwidth
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'radius_rot', &
& val = me%radius_rot, &
& ErrCode = iError )
write(logUnit(1),*) ' * radius_rot =', me%radius_rot
me%radius_C = me%radius_rot / 3._rk
!amplitude
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'circulation', &
& val = me%circulation, &
& ErrCode = iError )
write(logUnit(1),*) ' * circulation =', me%circulation
me%Ma = me%circulation/(4._rk*pi*me%radius_rot*me%cs)
! some more info to the crvp initialization
write(logUnit(1),*) ' * Ma_rot =', me%Ma
! Set the cutoff definition
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'cutoff_length', &
& val = me%cutoff%length, &
& ErrCode = iError )
if ( btest( iError, aotErr_NonExistent ) ) then
me%cutoff%active = .false.
else
me%cutoff%active = .true.
endif
!me%cutoff%r_min = 0.6_rk
!me%cutoff%r_max = 0.95_rk
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'cutoff_rmin', &
& val = me%cutoff%r_min, &
& ErrCode = iError, &
& default = huge(me%cutoff%r_min))
call aot_get_val( L = conf, &
& thandle = thandle, &
& key = 'cutoff_rmax', &
& val = me%cutoff%r_max, &
& ErrCode = iError, &
& default = huge(me%cutoff%r_max))
if( me%cutoff%active ) then
me%cutoff%linear = .false.
me%cutoff%quadratic = .true.
write(logUnit(1),*)' * Cutoff length = ', me%cutoff%length
if( me%cutoff%linear ) then
write(logUnit(1),*)' * Cutoff progress linear'
end if
if( me%cutoff%quadratic ) then
write(logUnit(1),*)' * Cutoff progress quadratic'
end if
endif
end subroutine load_ic_2dcrvp