Boundary condition for the block face with zero area; turning into a pole
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
character(len=*), | intent(in) | :: | face | Name of the face at which boundary condition is called |
||
type(boundarytype), | intent(in) | :: | bc | |||
type(extent), | intent(in) | :: | dims |
subroutine pole(face, bc, dims)
!< Boundary condition for the block face
!< with zero area; turning into a pole
implicit none
type(extent), intent(in) :: dims
type(boundarytype), intent(in) :: bc
character(len=*), intent(in) :: face
!< Name of the face at which boundary condition is called
call copy3(density , "flat", face, bc, dims)
call copy3(x_speed , "flat", face, bc, dims)
call copy3(y_speed , "flat", face, bc, dims)
call copy3(z_speed , "flat", face, bc, dims)
call copy3(pressure, "flat", face, bc, dims)
select case (turbulence)
case('none')
!do nothing
continue
case('sa', 'saBC')
call copy3(tv, "flat", face, bc, dims)
case('sst', 'sst2003')
call copy3(tk, "flat", face, bc, dims)
call copy3(tw, "flat", face, bc, dims)
case('kkl')
call copy3(tk, "flat", face, bc, dims)
call copy3(tkl, "flat", face, bc, dims)
case DEFAULT
Fatal_error
end select
select case(trim(transition))
case('lctm2015')
call copy3(tgm, "flat", face, bc, dims)
case DEFAULT
continue
end select
end subroutine pole