subroutine setboundarycondition_DGModel1D_t(this)
! Here, we use the pre-tendency method to calculate the
! derivative of the solution using a bassi-rebay method
! We then do a boundary interpolation and side exchange
! on the gradient field
implicit none
class(DGModel1D_t),intent(inout) :: this
! local
integer :: ivar
integer :: N,nelem
real(prec) :: x
nelem = this%geometry%nelem ! number of elements in the mesh
N = this%solution%interp%N ! polynomial degree
! left-most boundary
if(this%mesh%bcid(1) == SELF_BC_PRESCRIBED) then
x = this%geometry%x%boundary(1,1,1)
this%solution%extBoundary(1,1,1:this%nvar) = &
this%hbc1d_Prescribed(x,this%t)
elseif(this%mesh%bcid(1) == SELF_BC_RADIATION) then
this%solution%extBoundary(1,1,1:this%nvar) = &
this%hbc1d_Radiation(this%solution%boundary(1,1,1:this%nvar),-1.0_prec)
elseif(this%mesh%bcid(1) == SELF_BC_NONORMALFLOW) then
this%solution%extBoundary(1,1,1:this%nvar) = &
this%hbc1d_NoNormalFlow(this%solution%boundary(1,1,1:this%nvar),-1.0_prec)
else ! Periodic
this%solution%extBoundary(1,1,1:this%nvar) = this%solution%boundary(2,nelem,1:this%nvar)
endif
! right-most boundary
if(this%mesh%bcid(1) == SELF_BC_PRESCRIBED) then
x = this%geometry%x%boundary(2,nelem,1)
this%solution%extBoundary(2,nelem,1:this%nvar) = &
this%hbc1d_Prescribed(x,this%t)
elseif(this%mesh%bcid(1) == SELF_BC_RADIATION) then
this%solution%extBoundary(2,nelem,1:this%nvar) = &
this%hbc1d_Radiation(this%solution%boundary(2,nelem,1:this%nvar),-1.0_prec)
elseif(this%mesh%bcid(1) == SELF_BC_NONORMALFLOW) then
this%solution%extBoundary(2,nelem,1:this%nvar) = &
this%hbc1d_NoNormalFlow(this%solution%boundary(2,nelem,1:this%nvar),-1.0_prec)
else ! Periodic
this%solution%extBoundary(2,nelem,1:this%nvar) = this%solution%boundary(1,1,1:this%nvar)
endif
endsubroutine setboundarycondition_DGModel1D_t