sourcemethod_DGModel1D_t Subroutine

public subroutine sourcemethod_DGModel1D_t(this)

Arguments

TypeIntentOptionalAttributesName
class(DGModel1D_t), intent(inout) :: this

Contents


Source Code

  subroutine sourcemethod_DGModel1D_t(this)
    implicit none
    class(DGModel1D_t),intent(inout) :: this
    ! Local
    integer :: iel
    integer :: i
    real(prec) :: f(1:this%solution%nvar),dfdx(1:this%solution%nvar)

    do concurrent(i=1:this%solution%N+1,iel=1:this%mesh%nElem)

      f = this%solution%interior(i,iel,1:this%solution%nvar)
      dfdx = this%solutionGradient%interior(i,iel,1:this%solution%nvar)

      this%source%interior(i,iel,1:this%solution%nvar) = &
        this%source1d(f,dfdx)

    enddo

  endsubroutine sourcemethod_DGModel1D_t