sourcemethod_DGModel2D_t Subroutine

  subroutine sourcemethod_DGModel2D_t(this)
    implicit none
    class(DGModel2D_t),intent(inout) :: this
    ! Local
    integer :: iel
    integer :: i
    integer :: j
    real(prec) :: s(1:this%nvar),dsdx(1:this%nvar,1:2)

    do iel = 1,this%mesh%nelem
      do j = 1,this%solution%interp%N+1
        do i = 1,this%solution%interp%N+1

          s = this%solution%interior(i,j,iel,1:this%nvar)
          dsdx = this%solutionGradient%interior(i,j,iel,1:this%nvar,1:2)
          this%source%interior(i,j,iel,1:this%nvar) = this%source2d(s,dsdx)

        enddo
      enddo
    enddo

  endsubroutine sourcemethod_DGModel2D_t