sourcemethod_DGModel1D Subroutine

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

    call gpuCheck(hipMemcpy(c_loc(this%solution%interior), &
                            this%solution%interior_gpu,sizeof(this%solution%interior), &
                            hipMemcpyDeviceToHost))

    call gpuCheck(hipMemcpy(c_loc(this%solutiongradient%interior), &
                            this%solutiongradient%interior_gpu,sizeof(this%solutiongradient%interior), &
                            hipMemcpyDeviceToHost))

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

        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
    enddo

    call gpuCheck(hipMemcpy(this%source%interior_gpu, &
                            c_loc(this%source%interior), &
                            sizeof(this%source%interior), &
                            hipMemcpyHostToDevice))

  endsubroutine sourcemethod_DGModel1D